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1.  INTRODUCTION 


Couplod  Iran;  mission  lines  have  continually  received  much  attention  in 
many  diverse  areas  of  application.  Multiconductor  transmission  lines  have 
been  investigated  in  early  power  system  studies  and  continue  to  receive 
attention  in  this  area  with  regard  to  the  transient  behavior  of  power  lines 
under  fault  and  lightning  induced  conditions  [l-13].  Modern  emphasis  on 
multilayer  distributed  circuits,  strip  lines  and  microstrip  associated  with 
integrated-circuit  technology  has  produced  c renewal  of  interest  [14-19,68] 
as  has  the  interest  in  predicting  transients  induced  on  cables  by  external 
electromagnetic  field  sources  such  as  high  power  radars  or  an  electromag- 
netic pulse  (EMP)  from  nuclear  detonations  [20-27].  Determining  cross- 
talk in  communication  circuits  [28-30]  and  digital  computer  wiring  inter- 
ference [30-32]  are  examples  of  other  areas  in  which  the  subject  of  multi- 
conductor transmission  lines  consistently  arise. 

Of  particular  interest  within  the  electromagnetic  compatibility  (EMU) 
community  is  the  prediction  of  coupling  between  wires  and  their  associated 
termination-networks  in  closely  coupled,  high  density  cable  bundles  and 
flat  pack  (ribbon)  cables  on  modern  electronic  systems.  Control  of  intra- 
system  electromagnetic  compatibility  for  systems  within  the  Department  of 
Defense  is  generally  governed  by  MID-STD-461  and  462.  These  are  general 
documents  which  prescribe  limits  on  emissions  and  susceptibilities  of  the 
individual  subsystems  and  equipments  with  regard  to  undesired  signals 
(interference)  and  do  not  in  themselves  consider  the  coupling  paths  between 


the  equipments  and  subsystems  within  systems.  The  undesired  signals  as 
used  in  this  context  are  with  respect  to  the  particular  equipment  or  sub- 
system, net  all  . t which  are  unuesired  from  the  overall  system  standpoint. 
For  example,  the  undesired  signals  may  be  truly  undesired  ones,  such  as 
transmitter  harmonics,  or  may  be  the  result  of  an  essential  signal,  such  as 
the  fundamental  frequr-  (,(.  of  a transmitter,  coupling  to  a receptor  for  which 
such  coupling  is  not  intended. 

Even  if  all  the  equipments  and  subsystems  within  a system  conform  to 
the  limits  in  MIL-STD-461,  it  is,  of  course,  not  necessarily  true  that  over- 
all system  compatibility  v/ill  be  achieved.  Since  these  limits  do  not  take 
into  account  the  various  coupling  mechanisms  and  proximities  of  the  equip- 
ments, a system  whose  equipments  and  subsystems  meet  MIL-STD-461  may 
prove  to  be  incompatible  and  numerous  instances  of  required  retrofit  and 
interference  suppression  measures  on  systems  meeting  these  limits  illus- 
trate this  fact.  Thus  overall  s ystem  compatibility  may  not  be  achieved 
unless  all  signals  (desired  ar»d  undesired)  and  actual  coupling  paths  within 
the  system  are  considered,  analytically.  This  deficiency  has  led  to  the 
development  of  various  computer-aided  intra- system  (as  ;<pposed  to  intei- 
system)  compatibility  prediction  programs  which  mathe/r.atic  lly  model  the 
systems  and  take  into  account  the  various  coupling  path.'  for  unintentional 
energy  transfer  (interference)  as  well  as  intentional  energy  transfer  [33.  37 ]„ 

The  various  coupling  paths  can  generally  be  classified  into  combinations 
of  wire,  antenna  and  metallic  box  coupling,  e.g.,  wire- to- wire,  antemia- 
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to-antenna,  antenna- to- wire,  box-U-box,  etc.  In  the  case  of  wire-to-wire 


coupled  interference  in  cable  bundles,  this  undesired  coupling  of  energy 
between  circuits  sharing  a common  bundle  may  be  more  severe  than  one 
ma^  realize.  For  example,  numerous  cases  (both  experiiixental  and  analyti- 
cal) may  be  shown  where,  for  certain  frequencies,  the  ratio  of  the  received 
interference  voltage  across  the  terminals  of  a device  to  the  voltage  ernitted 
by  another  device,  which  is  coupled  via  wire-to-wire  coupling  mechanisms, 
exceeds  unity.  The  two  devices  are  not  directly  connected  by  a common 
pair  of  wires;  the  wdres  connected  to  each  device  are  only  in  close  proximity 
in  a common  cable  bundle.  Rarely  does  one  encounter  voltage  transfer 
functions  with  magnitudes  greater  than  unity  in  antenna- to-antenna  inter- 
ference coupling  problems  and  this  illustrates  the  importance  of  considering 
the  mechanism  of  wire-coupled  interference  transfer. 

It  is  the  purpose  of  this  report  to  provide  a cc.nplete  and  unified  discus- 
sion of  multiconductor  transmission  line  theory  as  it  applies  to  the  predic- 
tion of  wdre-coupled  interference.  The  common  approaches  and  assumptions 
which  are  either  explicitly  or  implicitly  used  in  the  problem  formulations 
which  appear  throughout  the  literature  are  discussed.  In  addition  to  provid- 
ing a discussion  of  the  limitations  and  advantages  of  each  of  these  techniques, 
some  numerically  stable  and  efficient  techniques  for  solving  the  multicon- 
ductor transmission  line  problem  for  large  numbers  of  closely  coupled, 
dielectric-insulated  wires  will  be  presented.  Methods  for  computing  the 
per- unit-length  parameters  will  also  be  given.  Some  of  the  results  . an  be 
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found  in  various  places  in  the  literature  although  the  treatments  of  the  sub- 
ject of  multiconductor  lines  generally  either  discuss  the  solution  of  the 
equations  describing  the  transmission  line  and  associated  termination*net- 
works  with  the  entries  in  the  transmission  line  equations  (the  per- unit-length 
parameters)  assumed  to  be  obtainable  or  they  di^^cuss  the  derivation  of  the 
per-unit-length  parameters  without  regard  to  the  solution  of  the  equations 
describing  the  line.  The  purpose  of  this  report  is  to  provide  a comprehen- 
sive discussion  of  the  complete  problem  solution  and  in  addition  present 
some  new  techniques  for  considering  large  numbers  of  closely  coupled, 
dielectric-insulated  wires. 

Throughout  this  report,  the  emphasis  will  be  on  the  frequency  respons-e 
of  the  transmission  lines  rather  than  the  transient  response  since  EMC  con- 
trol documents  currently  apply  predominantly  to  the  frequency  domain.  If 
one  assumes  linear  termination  networks  (no  hysteresis,  etc.)  and  assumes 
no  nonlinear  effects  associated  with  the  transmission  lines  such  as  corona 
discharge,  then  the  equations  describing  the  problem  (the  transmission  lines 
and  associated  terminations)  will  be  linear  and  thus  the  frequercy  response 
provides  a completely  general  characterization. 

Matrix  formulation  of  the  equations  and  other  results  of  matrix  analysis 
will  be  used  where  necessary  for  a logical  and  concise  development  and  the 
reader  is  referred  to  [38]  or  other  texts  on  linear  algebra  listed  in  the  refer- 
ences. 
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U.  THE  TEM  MODE  FORMULATION  FOR  MULTI  CONDUCTOR  LINES 


Consider  a Ax  length  section  of  an  (n+l)-conductor,  uniform  transmis- 
sion line  in  a homogeneous  medium  shown  in  Fig.  1 lying  parallel  to  the  x 
direction  in  a rectangular  coordinate  system.  The  line  is  said  to  be  uniform 
if  there  is  no  cross-sectional  variation  with  x either  in  the  conductors  or  the 
characteristics  of  the  medium,  i.e.  , "end-on"  or  cros£<- sectional  views  in 
planes  perpendicular  to  x are  identical  for  all  x.  The  medium  surrounding 
the  conductors  and  contained  within  the  zero-th  conductor  is  assumed  to  be 
linear  and  isotropic  and  therefore  is  describable  by  the  scalars  e (permit- 
tivity), gt  (permeability),  and  a (conductivity)  which  are  independent  of  the 
electric  and  magnetic  fields  in  the  medium  but  may  be  functions  of  frequency. 
If  e,  u and  a are  independent  of  position  in  the  medium,  i.e.,  independent  of 
X,  y and  z,  the  medium  is  said  to  be  homogeneous.  Thus  for  uniform  lines, 
all  (n+1)  conductors  have  uniform  cross  sections  along  their  lengths  and  are 
parallel  to  each  other  and  the  x direction  and  in  the  case  of  an  inhomo- 
geneous medium,  the  characteristics  of  the  medium  (c,  u.  O')  exhibit  no 
cross-sectional  variation  with  x and  are  therefore  independent  of  x. 

The  conventional  distributed-parameter,  transmission  line  model,  of 
course,  describes  only  the  TEM  (Transverse  Electro-Magnetic)  mode  of 
propagation  on  the  line  and  higher  order  modes  are  not  considered.  The 
electric  field  intensity  vector,  P(x,  y,  z,t),  and  the  magnetic  field  intensity 
vector,  W(x,  y,  z,  t),  for  the  TEM  mode  of  propagation  both  lie  in  planes  (y,  z) 
transverse  or  perpendicular  to  the  direction  of  propagation  (the  x direction) 
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Figure  1.  An  {n+l)-conductor  uniform 
transmission  line. 
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and  t is  the  time  variable.  Thus  it  has  been  shovn  a number  of  times  that, 
assuming  (n-fl)  perfect  conductors,  a homogeneous  medium  and  the  TEM 
mode  of  propagation,  the  nonzero  components  of  the  field  vectors  (the  trans.. 
verse  electric  field,  P,j,(x,  y,  z,  t),  and  the  transverse  magnetic  field, 

-4 

K.p(x,  y,  z,  t))  at  each  x along  the  line  satisfy  the  same  spatial  distributions  as 
static  fields  [40].  Therefore  one  can  meaningfully  define  voltages  between 
the  conductors  and  currents  flowing  on  the  conductors  [40],  For  further 
clarification,  see  Appendix  A. 

The  emphasis  in  this  report  will  be  upon  determining  the  frequency 
response  of  the  transmission  lines  and  associated  termination-networks. 
Therefore  sinusoidal  excitation  is  a ssumed  with  the  field  vectors  written  as 
y,  z,  t)  = E(x,  y,  z)eJ^  ^ and  il{x,  y,  z,  t)  = ll(x,  y,  z)ej“^  ^ where  E(x,  y,  z)  and 

-» 

H(x,  y,  z)  are  complex- valued  vectors  independent  of  time  t and  uu  is  the 
radian  frequency  of  excitation  (a'  = 2Trf).  To  characterize  lines  in  a homo- 
geneous medium  such  as  in  Fig.  1 under  the  TEM  mode  assumption,  the 
potential, (x,  t),  of  the  i-th  conductor  with  respect  to  the  reference  con- 
ductor (the  zero  conductor)  and  the  current,  «^^(x,  t),  associated  with  the 
i-th  conductor  are  defined  for  i=l, --,n  (see  Fig,  Ic).  The  currents  are 

directed  in  the  positive  x direction  and  t!.e  current  in  the  reference  conduc- 

n 

tor  satisfies  J^(x,  t)  = - v ^^(x,t)  [40].  Voltages  and  currents  for  sin- 
usoidal excitation  aie  written  asZ^(x,  t)  = Vj^(x)e-^'^^  and^j^(x,  t)  = L(x)ej'^^ 
where  V.(x)  and  I-(x)  are  the  phasor  voltages  and  currents  respectively  and 
are  complex- valued  scalars  independent  of  time,  t.  In  the  cross-sectional 
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view  of  Fig.  Ib,  the  voltage  of  the  i-th  conductor  with  respect  to  the  zero-th 

conductor  chosen  as  a reference  is  defined  as  the  lir,-*  integral  of  along 

contour  C-  in  the  y,  z plane  and  the  current  associated  with  the  i-th  conductor 

-*  * 

is  defined  as  the  line  integral  of  along  the  closed  contour  Cj  in  the  y,  2 
plane.  The  assumption  of  TEM  mode  propagation  precludes  the  existence  of 
^ component  of  the  magnetic  field  intensity  vector  in  the  longitudinal  direc- 
tion (the  K direction).  This  assumption  coupled  with  the  assumption  of  per- 
fect conductors  insures  that  the  definition  of  the  voltage:^  is  unique  [40],  The 
assumption  of  a TEM  fields  structure  also  precludes  the  existence  of  a longi- 
tudinal component  of  the  electric  field  intensity  vector.  Theraforej  no  longi- 
tudinal cond>!ction  or  displacement  current  in  the  dielectric  is  considered 
and  any  current  flow  in  the  dielectric  will  be  confined  to  the  transverse 
plane.  This  assumption  coupled  with  the  assumption  of  perfect  coj.ductors 
insures  that  the  definition  of  the  line  currents  is  unique  [40].  These  results, 
of  course,  provide  the  basis  for  representing  transmission  lines  for  the  TEM 
mode  of  propagation  over  "electrically  short"  Ax  lengths  with  lumped  equiva- 
lent circuits  whose  parameters,  which  are  per-unit-length  qviantities  and  are 

-♦  -» 

derived  under  the  condition  that  the  transverse  iield  vectors,  .and  at 
each  X along  the  line  satisfy  static  distributions,  represent  the  TEM  mode 
of  propagation  for  non-static  excitation  [40].  These  important  conclusion  5 
are  demonstrated  in  Appendix  A. 

Imperfect  conductors,  inhomogeneous  media  and  electrically  large 
cross-sectional  line  dimensions  preclude  the  existence  of  only  the  TEM 


mode  for  the  following  reasons.  With  lossy  conductors,  there  will  neces» 
sarily  be  a longitudinal  component  of  the  electric  field  in  the  x direction  due  > 

to  the  nonzero  surface  impedance  of  the  conductors  [40].  If  the  surrounding 
medium  is  inhomogeneous,  then  wave  propagation  can  no  longer  be  TEM  as 
a result  of  the  different  phase  velocities  in  the  different  honnogeneous  por- 
tions of  the  media.  Imperfect  conductors  and  inhomogeneous  media  are 
nevertheless  considered  with  th«t:  distributed-parameter,  transmission  line 
model  under  the  assumption  that  the  conductor  losses  and  the  inhomogenei- 
ties in  the  media  do  not  significantly  perturb  the  field  distribution  from  a 
TEM  structure.  The  inclusion  of  inhomogeneous  media  which  is  termed  the 
"quasi-TEM  mode"  assumption  is  particularly  important  in  microstrip 
problems  and  other  associated  integrated-circuit  structures  [14-18,68]. 

Electrically  large  cross-sectional  dimensions  of  the  line  (conductor  separa-  * 

tion,  wire  radius,  etc.)  evidently  are  also  capable  of  producing  higher  order 
modes  a_.u  this  can  be  surmised  from  the  Cact  that  the  infinite  parallel-plate 
transmission  Une,  which  is  rigorously  solvable  and  capable  of  supporting 
the  TEM  mode  of  propagation,  will  support  only  the  TEM  mode  for  frequen- 
cies such  that  the  plate  spacing  is  less  thv.n  onc-half  wavelength.  Also,  it 
can  be  shown  that  a two- conductor  coaxial  line  will  support  nig  her  order 
modes  when  the  mean  circumference  of  the  annular  space  between  the  two 
conductors  is  greafer  than  one  wavelength.  Thus  throughout  this  report, 
tlie  cross- sectional  dimensions  of  the  line  vrill  be  assumed  to  be  electri- 
cally small,  i.  e.  , much  less  than  a wavelength,  so  that  transmission  line 
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theory  applies,  i.e.,  the  TEM  mode  is  the  dominant  mode  of  propagation. 

The  specific  cases  of  interest  to  be  considered  in  this  report  are  shown 
as  cross-sectional  views  in  the  y,  z plane  in  Fig,  2 and  Fig.  3.  In  Fig.  2, 
n wiret’  (circular  conductors)  are  shown  with  another  conductor,  the  refer- 
ence conductor,  denoted  as  the  zero-th  conductor.  In  Fig.  2a,  the  reference 
conductor  is  also  a wire  whereas  in  Fig.  2b  and  Fig.  2c  the  reference  con- 
ductors are  an  infinite  ground  plane  and  an  overall  circular  shield  respec- 
tively. These  lines  are  uniform  and  the  surrounding  medium  is  homogen- 
eous. In  Fig.  2a  and  Fig.  2b,  the  surrounding  medium  is  free  space  with 
parameters  and  cty,.  In  Fig.  2c,  the  medium  within  the  circular  shield  is 
homogeneous  with  parameters  e»  o.  (The  permeability  of  all  dielec- 

trics in  this  report  will  be  considered  to  be  that  of  free  space,  u^*) 

In  Fig,  3,  similar  cases  are  shown  with  the  wires  having  circular 
dielectric  insulations  (an  obviously  very  common  situation).  Thus  the 
medium  in  each  of  these  cases  is  inhomogeneous  although  the  lines  are 
nevertheless  uniform.  The  permeabilities  of  the  dielectric  insulations  are 
considered  to  be  that  of  free  space,  as  is  lypi<-al  of  difr,-le':  .‘■rics.  Each 

dielectric  insulation  is  described  by  the  scalars  permittivity  , £ , emd  con- 
ductivity, o i=0, 1, -,n  and  the  space  surrounding  the  dielectr.v  ''■ifjula... 

tions  is  considered  to  be  free  space. 

The  corresponding  cases  for  the  more  familiar  two-conductor  lines 
(n=l)  are  shown  in  Fig.  4 and  Fig.  5.  Note  in  Fig.  4 that  the  lines  of  E and 
H are  shown  perpendicular  to  each  other.  This  is  a natural  consequence  of 


Figure  3.  Multiconductor  transmission  lines  in  an 
inhomogeneous  medium. 


Figure  5. 


Two-conductor  transmission  lines 
III  an  inhomogeneous  medium. 
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the  TEM  mode  assumption  [40]. 

If  the  medium  is  homogeneous  as  in  Fig.  1 and  Fig.  2 and  all  (n+1)  con- 
ductors are  perfect  conductors,  then  losses  in  the  medium  can  be  included 
without  violating  the  TEM  mode  assumption  or  the  uniqueness  of  the  voltage 
and  current  definitions  [40].  However,  in  the  case  of  a homogeneous  me- 


dium in  Fig,  2,  it  is  only  logical  to  consider  a lossy  medium  for  the  case  in 
Fig.  2c  since  the  surrounding  medium  in  Fig.  2a  and  Fig.  2b  is  considered 
to  be  free  space.  Dielectric  losses  can  be  introduced  through  a finite,  non- 


zero ohmic  conductivity,  a (which  generally  will  be  quite  small  for  typi- 
cal insulation  materials)  and  also  through  dipole  relaxation  effects  [30].  To 
include  both  of  these  effects,  we  may  consider  the  material  to  be  charac- 
terized by  a complex,  effective  permittivity  (which  is  frequency  dependent) 
instead  of  a real  permittivity.  To  include  dipole  relaxation  losses,  the  per- 
mittivity may  be  considered  to  be  complex  as  [30]  e = d -jc".  Ampere's 

law  in  a homogeneous  medium  possessing  both  of  these  loss  quantities  be- 

(C7  1 + tJU  e")  ■* 

comes  vXH  = CTj  E + joue  E = [(a^  + uje")  + Jcue'JE  = jujc'  [ 1-j  ; ]E. 

The  real  part  of  the  complex  permittivity  is  expressed  as  e'  = where 

is  the  permittivity  of  free  space  and  Ej.  is  the  relative  dielectric  constant. 
The  effective  conductivity  of  the  ho  ogeneous  medium  then,  becomes 
a = + U)e"»  Thus  the  losses  of  the  medium  may  be  accounted  for  by  using 

a complex  effective  permittivity  e^££  = 6)  instead  of  a real  per- 

mittivity and  tan  6 = a/(uu£^  E^.)  is  the  loss  tangent  of  the  material  [40]. 
Ordinarily,  the  loss  tangent  and  the  relative  dielectric  constant  e^.  i^-re 
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given  for  materials  as  a function  of  freque,.:\cy.  Therefore,  it  is  quite  clear 
that  for  (n+1)  perfect  conductors  in  a homogeneous  medium,  losses  in  the 
medium,  i.e.,  a ^0,  can  be  included  without  violating  the  TEM  mode 
assumption  or  the  uniqueness  of  the  voltage  and  current  definitions  since  the 
real  permittivity  for  the  lossless  case  (o  =0)  is  merely  replaced  by  a com- 
plex permittivity,  e^ff#  to  account  for  losses  in  the  medium.  Since  the  TEM 
mode  assumption  is  legitimate  for  the  lossless,  homogeneous  case,  there  is 
no  reason  why  the  use  of  a complex  permittivity  instead  of  a real  permit- 
tivity should  change  this. 

The  lumped-circuit  model  for  a ^x  length  section  of  the  two-conductor 
lines  in  a homogeneous  medium  in  Fig.  4 are  shown  in  Fig.  6.  The  lines 
have  a total  length  £ and  Thevenin  equivalents  of  the  linear  terminations  at 
the  ends  of  the  line  are  shown. 

The  lumped-circuit  model  describing  the  TEM  mode  of  propagation  for 
a ^x  length  section  of  any  of  the  multiconductor  lines  in  a homogeneous 
medium  in  Fig.  1 and  Fig.  2 is  shown  in  Fig.  7.  All  t.  x length  models  fox 
other  sections  of  the  line  will  be  identical  since  the  line  is  uniform.  Since 
the  cross-sectional  dimensions  of  the  line  (conductor  spacing,  wire  radius, 
etc.)  are  all  assumed  to  be  "electrically  small"  and  ^x  is  assumed  to  be 
"electrically  short",  then  it  is  valid  to  characterize  a ^x  section  of  the  line 
with  a lumped  equivalent  circuit. 

Resistance  elements  r.,  , r , r^  and  conductance  elements  g;n. 

Cq  Cj  ”l0 

g--  are  included  to  represent  losses  associated  with  the  conductors  and 

* J 
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Figure  6.  The  termination-networks  and  equivalent  cir'^uits  for 
two- conductor  transmission  lines. 

-22- 
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medium  respectively.  The  inclusion  of  a surrounding  medium  having  a 
finite,  nonzero  conductivity  and  dipole  relaxation  losses  with  these  shv.’nt 
conductances  is  consistent  with  the  assumption  of  TEM  mode  propagation 
whereby  no  longitudinal  conduction  or  dieplacement  current  can  flow  in  the 
dielectric  and  any  current  flow  in  the  medium  is  confined  to  the  transverse 
plane.  The  shunt  conductances  account  for  the  portions  of  the  transverse 
currents  a.ssociated  with  conductive  and  dipole  relaxation  losses  of  the 
medium,  i.e.,  the  transverse  displacement  and  corductic  currents  due  to 
the  imaginary  part  of  Similarly,  shunt  capacitances  account  for  the 

transverse  displacement  currents  associated  with  the  real  part  of 
Also  self  inductance  terms  for  the  conductors,  a , i.,  i mutual  inductances 

0 i j 

between  the  conductors,  mutual  capacitances  between  the 

conductors,  c^q,  Cj^,  c.^,  are  shown  [39].  Lcysy  conductors  also  produce  a 

portion  of  the  self  inductances  due  to  skin  effect  which  is  represented  by  the 

elements  f which  are  internal  self  inductances  produced  by  cur- 

^0 

rents  internal  to  the  lossy  conductors  [2,3,30].  The  infinite  ground  plane 
and  circular  shield  in  Fig.  2b  and  Fig.  2c  are  considered  to  be  perfect 
conductors  and  for  these  cases  r = l = 0.  A method  of  including  a lossy 
ground  plane  is  given  in  [29]  and  is  frequently  used  to  represent  the  earth 
return  path  in  power  systems  [13]. 

Some  care  must  be  exercised  in  interpreting  the  elements 
and  *^jo»  strictly  "self  inductances"  and  "mutual  inductances" 

respectively  in  the  conventional  sense.  This  interpretation  relies  on  the 
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property  that  the  sum  of  the  currents  (at  a particular  x)  associated  with  all 


n 


(n+1)  confluctors  is  zero,  i.  e.  , T I.(x)  = 0.  An  excellent  discussion  of  this 

i=0  '■ 

is  presented  in  reference  [3],  Chapter  1 and  the  reader  is  referred  to  this 
for  further  clarification.  Our  results  will  not  rely  on  this  interpretation 
since  we  will  not  determine  these  individual  external  inductance  parameters 
but  will  instead  obtain  the  per- unit-length  external  inductance  matrix,  L.,  of 
the  line  directly.  The  entries  in  L,  which  are  the  essential  items  in  our 
analysis,  will  be  linear  combinations  of  these  per-unit-length  "inductances" 
and  once  L is  determined,  there  is  no  need  to  separate  its  entries. 

All  of  the  terms  resulting  from  losses,  r^^,  r^  , r^^,g^Q,  g^^, 

I t are,  in  general,  functions  of  frequency.  The  external  parameters, 
Iq,  a.,  ly  m...  m.Q,m.Q,  c.^,  c.^,  c.q,  g.^,  g.^,  g..,  are  derived  assuming 
perfect  conductors  such  that  the  transverse  fields  satisfy  a static  distribution 
at  each  x along  the  line  [39].  These  external  parameters  will  also  be  func- 
tions of  frequency  if  the  permeability,  permittivity  or  conductivity  of  the 
surrounding  medium  is  a function  of  frequency.  In  this  case,  the  parameters 
are  recomputed  for  each  frequency  assuming  the  transverse  fields  satisfy  a 
static  distribution  at  each  x along  the  line.  All  parameters  are  per-vinit- 
length  quantities  and  therefore  the  total  value  of  each  parameter  for  a 
length  model  in  Fig.  7 is  the  per-unit-length  value  multiplied  by  the  section 
length,  hx. 

It  is  important  to  note  that  this  is  an  exact  representation  of  the  TEM 
mode  of  propagation  for  (n+1)  perfect  conductors  in  a homogeneous  medium 
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as  in  Fig.  1 and  Fig.  2,  Imperfect  conductors  are  considered  as  an  approxi- 


mation through  r , r , r , under  tli-,  assumption  that  the 

Cq  Cj  CQ  V.£  Uj 

conductivities  of  the  conductors  are  very  large  and  much  greater  that  the 
conductivity  of  the  dielectric  medium  so  that  the  fields  structure  is  essen- 
tially TEM.  Although  the  presence  of  an  inhomogeneous  medium  as  in  Fig. 

3 precludes  the  existence  of  the  TEM  mode  except  perhaps  in  the  limiting 
case  of  zero  frequency,  the  equivalent-circuit  representation  in  Fig.  7 will 
be  assumed  to  be  an  adequate  representation  for  the  quasi-TEM  mode  for  the 
lines  in  an  inhonaogeneous  medium  in  Fig.  3.  The  parameters  for  this  case 
will  also  be  computed  at  each  frequency  by  assuming  (as  a first-order 
approximation)  that  the  field  vectors  are  entirely  transverse  and  satisfy  a 
static  distribution  at  each  x along  the  line. 

For  the  two- conductor  cases  in  Fig.  4,  the  transmission  line  equations 
can  be  derived  from  the  L x equivalent  circuits  in  Fig.  6 for  the  sinusoidal, 
steady  state  in  the  limit  as  -♦  0 as  a pair  of  coupled,  first-order,  ordi- 
nary, complex  differential  equations  [2,3] 

(la) 


^'dx^  + (^c  + ® 


vj: 

dl(x) 

dx 


+ (g  + ja'c)  V(x)  = 0 


(lb) 


where  Z and  Y are  the  per-unit-length  impedances  and  admittances  of  the 

line  respectively.  For  each  of  these  cases,  r_=  r„  + r_  , , 

1 *^0  1 0 

f + Xq  “ ^^0*  ^ ” ^10  ® ~ ®10*  ^ incident  electromagnetic  field 

illuminates  the  line  of  Fig.  4a,  the  equations  in  (1)  are  modified  to  include 
the  effects  of  the  incident  field  and  become  [20] 
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(2a) 


<1  V(x) 

4 Z I(x)  = V^Cx) 


•2^  + Y V(x)  = I,(x) 


(2b) 


v^^ere  V_(x)  and  I (x)  are  distributed  sources  along  the  line  induced  by  the 

S 8 

Spectral  components  of  the  incident  field  and  are  given  by  [20] 

>d 


pa  (inc) 

V (x)=ju,u^  \ H (y,x)dy 
s ^ 0 


(3a) 


Ig(x)=-Y\J  Ey(y,x)dy 


(3b) 


The  two  wires  in  Fig.  4a  lie  in  the  x, y plane  with  wire  0 at  y = 0 and  wire  1 

at  y = d.  The  components  of  the  incident  magnetic  and  electric  field  intensi- 

(inc) 

ties  at  the  radian  frequency  (u  in  the  z and  y directions  are  denoted  by  H (y>^ 
(inc) 

and  Ey^y,x)^ respectively. 

Similarly  for  multiconductor  lines,  the  transmission  line  equations  can 
be  derived  from  the  equivalent  circuit  in  Fig.  7 for  the  sinusoidal,  steady 
state  in  the  limit  as  ^x  -»  0 as  a pair  of  n coupled,  first-order,  ordinary, 
complex  differential  equations  in  matrix  form  as  (see  Appendix  B) 


V(x)  + Z I(x)  = V„(x) 


Ux)  + Y V(x)  =Jtg(x) 


(4a) 

(4b) 


vdiich  may  be  written  in  an  alternate  form  as  a set  of  2n  coupled  equations 
in  partitioned  form  as 


V(x) 

J(x) 

n^ 


0 


V(x) 

I(x) 


V^(x) 

i^(x) 


(ic) 
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A matrix,  M,  with  m rows  and  n columns  is  said  to  be  m y n and  the  ele- 
ment  in  the  i-th  row  and  j-th  columiii  is  designated  [m]..  with  i=l, ---,m 

•<w  \ 1 


and  j=l, The  dot  (•)  denotes  first  derivative  with  respect  to  x,  i.e.. 


position,  i.e.,  [ 0 ]. . = 0 for  i=l,---,mand  j=l,---,n.  The  elements  of 

m<^n  *-J 


the  n ^ 1 complex  column  vectors  V(x),  I(x),  V (x),  I (x)  are  [V(x)].  = V.(x), 
[l(x)].  = I.(x),  [v  (x)].  = V (x),  [l  (x)].  = I„(x)  where  the  element  of  an  n v 1 

— LI  — S I —5  I 

column  vector  JV  with  n rows  in  the  i-th  row  is  denoted  by  [^]j.  for  i=l, , n. 

The  per-unit-length  series  voltage  sources,  V_  (x),  and  shunt  current 

i 

sources,  I (x),  are  induced  by  the  spectral  components  of  the  incident 
®i 

field  and  are  complex. valued  and  functions  of  frequency  and  position,  x, 
along  the  line.  For  (n+1)  wires  in  a homogeneous  medium  in  Fig.  2a,  these 
sources  are  shown  in  Appendix  C and  in  [27]  to  be 

, r-^iO  (inc) 

Vg,(x)  = ^ (§.,x)d?.  (5a) 

0 


=io'+ 


j=i 

j^i 


'(5..X)  d?.  (5b) 


n 

+ T 

j=l 


{< 


gij  + j(D  C 


x)d?. 


where  5^  is  a straight-line  contour  between  wire  0 and  wire  i and  perpendic- 
ular to  wire  0 and  wire  i.  K^^”^\§.,x)  and  , x)  are  the  components  of 

Hi  I i 

the  incident  field  vectors  normal  to  a plane  formed  by  the  two  wires  and 
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parallel  to  (transverse  field)  respectively.  Solutions  for  Vg  (x)  and  I (x) 

i i 

for  ‘he  other  configurations  are  discussed  in  [22,  24,  25]  and  Appendix  C. 

The  n y n complex- valued  noatrices  Z and  Y are  the  per-unit-length 

impedance  and  admittance  matrices  respectively  and  are  symmetric,  i,  e.  , 
t t 

Z = Z and  Y = Y where  the  transpose  of  an  n y n matrix  M is  denoted  as 
M*".  These  matrices  are  independent  of  x since  the  lines  are  uniform  and 


are  separable  as 


Z = R + J OD  L + j u'  L 
Y = G + j(i)  C 


where  and  are  the  per-unit-length  conductor  resistance  and  conductor 
internal  inductance  matrices  respectively  and  are  real,  symmetric.  The 
external  parameter  noatrices,  G,  L and  C,  are  real  and  can  also  be  shown 
to  be  symmetric  (for  linear,  isotropic  media)  regardless  of  whether  the 
medium  is  homogeneous  or  inhomogeneous  thus  permitting  the  equivalent 
circuit  representation  in  Fig.  7.  [39].  The  matrices  G,  L and  C are  the 


per-unit-length  external  conductance,  inductance  and  capacitance  matrices 
respectively.  The  entries  in  these  matrices  are  obtained  in  Appendix  B 


and  are  given  by 


[R  ]..  = r + 

'^C  ll  C: 


• [R]..  = r^ 

^0  ^0 


[L  ]..  = X + /.  [L  ]. . = £ 
‘■-c-'ii  c.  cj-  ‘•'^chj 


[L].,  = X.  + £ - 2m.^  [l].  . = + m.  . - m.  - -m_ 

~Ti  i *0  iO  ~ ij  0 ij  iO  jO 
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fSlii=«iO+  SiJ  ‘Sly  (7d) 

Ji 

+ 

1 

for  i,  j = 1,  ---,n.  C and  G are  said  to  be  hyperdominant  since  each  term  on 
the  main  diagonal  is  greatei  than  the  sum  of  the  elements  in  that  row  [39] 
and  they  can  therefore  be  shown  to  be  positive  definite  meanin);  that  all  n 
eigenvalues  of  C and  all  n eigenvalues  of  G ai"e  positive  and  nonzero  [4l], 

The  derivation  of  the  per-unit-length  parameters  will  be  discussed  in  Section 
V. 


III.  SOLUTION  OF  THE  TRANSMISSION  LINE  EQUATIONS 


The  set  of  2n  first-order,  complex- valued,  ordinary  differential  equa- 
tions in  (4c)  \rfiich  describe  the  transmission  line  for  the  TEM  mode  of  pro- 
pagation and  the  sinusoidal,  steady  state  are  in  the  form  of  state  variable 
equations  [38,42],  Systems  of  first-order  differential  equations  in  the  state 
variable  form  have  received  considerable  attention  in  recent  years  in  the 
general  area  of  linear  systems  and  the  solution  to  (4c)  is 

(8) 

AA^ere  4 (x,x  ) is  the  2n  y 2n  comple.c- valued  state  transition  matrix  which  is 

the  solution  to  (4c)  with  Vg(x)  =_Ig(x)  = ^0  ^ and  the  parameter  Xq  is  some 

arbitrary  fixed  point  along  the  line  [38,42], 

Obviously  the  difficult  portion  of  the  analysis  (aside  from  the  difficulty 

in  computing  the  per- unit-length  parameters  and  equivalent  field  excitation 

sources  in  V (x)  and  I (x))  is  the  determination  of  the  state  transition 
— «s  ••••§ 

matrix  or  chain  parameter  matrix,  $(x,Xq).  Fortunately,  for  uniform  lines 
where  Z and  Y are  not  functions  of  x,  the  solution  is  fairly  simple  as  will  be 
shown  (although  there  are  some  important  computational  problems  when 
losses  are  inclv’ded).  For  nonuniform  lines  where  Z and  Y are  functions  of 
X,  i.e.,  Z(x)  and  Y(x),  (4c)  becomes  a set  of  nonconstant-coefficient  differen- 
tial  equations  (Bessel's  equation  is  an  exa.nple  of  a nonconstant-coefficient 
differential  equation)  [43-46].  For  these  types  of  lines,  (8)  holds  but  the 
ultimate  difficulty  is  the  determination  of  the  state  transition  matrix  and 


V(x) 

= *(3t.XQ) 

A 

+ \ $(x,x) 

A 

i(x) 

i<*0> 

X 

^0 

i^(x) 
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except  for  some  very  special  structures  one  must  resort  to  numerical 


methods  and  approximations  to  obtain  £(x,  x^)  [42],  If  the  line  is  "abruptly 

nonimiform"  as  with  branched  cables,  i.  e. , consists  of  uniform  subsections 

iii  cascade,  then  the  overall  chain  parameter  rixatrix,  $(x, x*),  is  the  pro- 

duct  (in  the  appropriate  order)  of  the  chain  matrices  of  the  individual 

uniform  subsections  between  x and  Xq  and  thus  is  straightforward  to  obtain. 

As  an  example  of  this  application  to  an  "abruptly  nonuniform"  line,  consider 

the  line  as  a cascade  of  N uniform  (n+l)-conductor  transmission  lines  with 

each  section  between  x = x.  and  x.  , described  by 

i i-l  ^ 


— * — 

r , -1 

V(x^) 

= Ii(^i»  ^i.l) 

>1 

A 

1 

fa 

i<*i> 

+ \ » (X  X) 

C?X 

(9) 

1 

i 

*U1 

[isi(-U 

4 f 

;( 

i 

X.  , < X < x. 
I- 1 — - I 


for  i=l,---,N  where  $:(x.,  x.  ,)  is  the  chain  parameter  matrix  for  the  i-th 

1 j,- 1 

section  betv'een  x = and  x =Xj^  ^^i-1^  si^^^^si  equiva- 

lent^induced  source  vectors  for  the  i-  th  section.  By  sequential  substitution,  the 
overall  chain  parameter  matrix  for  the  cascade  of  N sections  between  x^  and 
Xj^  (which  are  not  required  to  be  identical)  becomes 


V(xn) 


i<*N-  *0* 


■I 


^N-1^  1n-1^*N-1»  *N-2^ 
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, X^)  $j(Xj, 


(iO) 


v(-o) 


■r 

m 

s 


+ 


c N r 

3 j£n^^n» 

"^N-i 


X) 


V^N^x) 


X. 

ii+i**i+i’  *■>  'I  *’ 


X. 


i-1 


" * 1 
V .(x) 

— 8l' 


I .(X) 
r«i  -• 


, A 

ldS£ 


The  overall  chain  parameter  matrix  for  this  cascade  of  nonidentical  line 
sections  between  x = Xq  and  x = x^^  is  identified  in  (10)  as  the  matrix  product 
$^(Xj^,  Xq).  Note  that  the  indicated  order  of  multiplication  of  the  individual 
chain  parameter  matrices  must  be  preserved  since  they  do  not  generally 
commute.  Lumped-element  networks  at  discrete  points  along  the  line  can 
also  be  incorporated  into  the  problem  by  writing  the  ^natrix  chain  parame- 
ters of  these  networks  and  including  them  appropriately  into  the  product  of 
the  chain  parameter  matrices  of  the  individual  uniform  sections  in  the  above 
manner. 

When  the  line  is  uniform  (as  is  being  considered  here)  where  Z and  Y 
are  independent  of  x,  the  state  transition  matrix,  ^(x,x^),  can  be  shown  to 
be  a function  of  only  one  variable;  the  difference  quantity  (x-Xq)  \4Z],  Thus 
for  uniform  lines,  we  may  denote  the  state  transition  or  chain  parameter 
matrix  as  $(x-x.).  The  state  transition  matrix  has  the  property  that 
$(Xq,  Xq)  =^2n  -izn  ^ identity  matrix  with  = 1 and 

[1,  ]. . = 0 for  i,  j =1,  2n  and  i^j  [38].  This  should  be  clear  from  (8)  by 

~tn  ij 

setting  x equal  to  Xq,  Additionally,  it  may  be  shown  that  ^"^(x, Xq)  = £(Xq,x) 
where  the  inverse  of  an  n \ n matrix  M is  denoted  by  M"^  and  therefore  the 
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inverse  of  the  chain  parameter  matrix  may  be  trivially  determined  [4Z], 

This  is  qixite  obvious  from  (8)  for  ^(x)  =_Ig(x)  = by  interchanging  the 
roles  of  X and 

For  two- conductor  lines  (n=l)  the  transmission  line  equations  become  a 
set  of  two  complex- valued,  ordinary  differential  equations  given  in  (1).  The 
solution  of  the  transmission  line  equations  for  two- conductor  lines  can  be 
obtained  quite  easily  by  differentiating  (lb)  with  respect  to  x and  substituting 


(la)  to  yield 


d»  I(x) 

=YZI(x) 

= I(x) 


(11) 


where  the  propagation  constant,  v*  is 

V = Jrz  (12) 

. 

The  solution  to  (11)  becomes 

I(x)  = e"^^  1+  - e^^  r (13) 

where  and  I”  are  complex,  undetermined  constants.  Substituting  (13)  into 
(lb)  yields 

V(x)  = + e'^^  I’}  (14) 

where  the  characteristic  impedance,  Z^,  is  given  by 

Z^  =V  /Y  =/z7y  . (15) 

To  find  the  solutions  in  the  time  domain,  multiply  (13)  and  (14)  by  to 

ob  tain 


V(x,t)  = rz^  1+)  + r Zr  e'j"'* " ''*’1-1 


(16a) 


t) 


(x,  t) 
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4(x.t)  = 1+]  , [o<j'"^  I-]  , (16b) 

' * « — ^ -» 

4>"^(x.t)  cl'(x,t) 

Therefore  the  total  solution  consists  of  waves  traveling  in  the  4-x  direction 
(forward-traveling  waves)  denoted  by  t)  and  ^^(x,  t)  and  waves  traveling 

in  the  -x  direction  (backward-traveling  wavec)  denoted  by^(x„  t)  and^~(x,  t). 
The  characteristic  impedance,  Zq,  is  the  ratio  of  the  voltage  and  current  in 
the  respective  -vaves. 

For  two-conductor  lines,  the  chain  par.  jter  iriatrix  is  2 y 2 and  can 
easily  be  shown  to  be  [2,  3] 


«(x,Xo)  = 


cosh 


(v(x-Xq)} 


L-^sinh{v(x.Xo)} 


-Zq  Sinh  /Y(x-Xo^h 

/ iF^ 

cosh  I v(x.Xq)£J 


vdiere  cosh  and  sinh  are  the  hyperbolic  cosine  and  sine  respectively. 

Note  that  the  determinant  of  the  chain  parameter  matrix  is  unity.  Knovring 
this  quantity,  the  solution  for  the  voltage  and  current  at  any  point,  x,  along 
the  line  can  be  found  from  (8)  in  terms  of  the  voltage  and  current  at  some 
reference  point,  Xq,  as 

V(x)  = cosh{v(x-XQ)J  V(xq)  - Zq  sinh|^Y(x-xo)jl(xQ)  (18a) 


_x 

+ [cosh^v(x-x)  J V^(x)  - sinh^v(x-x)j  Ig(x)]  dx 
*0 

I(x)  = - sink  [ Y(x-  *0>}  V(Xq)  + cosh|'Y(x-XQ  )} 

c 

+ ^ [-  Y—  sinh  j^Y(x-x)  J Vg(x)  + cosh  j^Y(x-x)J  Ig(x)|  dx 


08b) 
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For  multiconductor  lines,  the  equations  in  (4c)  may  be  thought  of  as 
"strongly-coupled"  state  variable  equations  since  the  block  off-diagonal 
terms,  Z and  Y,  are  nonzero  vdiereas  the  block  main-diagonal  terms  are 
zero,  0 . The  chain  parameter  matrix,  ^(x,x^),  however,  may  be  dater- 
mined  in  the  following  manner  which  is  similar  to  the  method  for  solving  the 
two-conductor  line  employed  above  [18,26,4?].  Assuming  for  the  moment 
that  V„(x)  = I (x)  = nS  1,  differentiating  the  second  set  of  equations  (4b)  again 
with  respect  to  x,  I(x)  = -Y_y(x),  and  substituting  the  first  set  (4a),  Mx)  = 

-Z  I(x),  one  obtains  the  set  of  n second-order  differential  equations 

T(x)  = Y Z _I(x)  . (19) 

Note  that  even  though  Y and  Z are  symmetric,  it  is  not  necessarily  true 
that  the  matrix  product  YZ  (or  Z Y)  will  be  symmetric. 

The  solution  of  (19)  is  usually  obtained  with  similarity  transformations 
[13,18,26,38,41,42,47,48],  which  is  referred  to  in  the  power  transmission 
literature  as  "modal  decomposition"  [13].  Define  a change  of  variables, 

I(x)  = T_I^(x)  where  T is  an  nyn  nonsingular,  complex- valued  matrix  and 
I^(x)  is  an  n\l  comp  lex- valued  vector  of  "mode  currents".  Substituting 
this  change  of  variables  into  (19)  yields 

*I  = T"W  Z T . (20) 

Suppose  there  exists  an  n-yn  similarity  transformation,  T,  which  dia- 
gonalizes  Y Z,  i.e., 

T"^  Y Z T = Y*  (21) 

where  Y®  is  an  nyn  diagonal  m-'^trix  with 
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(22a) 


= V,» 


f;!:’)  ij  = 0 

i/j 


(22b) 


and  the  terms,  y®,  i = 1.  » a*..  » i 

p r^,  I X, ,nare  complex* valued 


becomes  a set  of  n uncoupled  diffe 
[18,26.47] 


scalars.  Then  (20) 


rential  .quationa  wieh  th.  aimple  aolution 


i(x)  = T I (X) 
m 


Vx 


where  U an  nXn  diagonal  matrix  with 


(23) 


^ u 


(24a) 

(24b) 


- 0 

co4(aa,™a.e„.U.aa  co...„.a, 

' i - fl  Ji  and  I ; = a general,  be  factions  of  frequency 
f47j.  These  undeCer,nlned  consfanf.  will  be  eralunfed  by  considering  ihe 
>»undary  condUions  or  fernrinaiion-nei^rhs  ai  fhe  ends  of  fbe  line.  Since 
^ (4b)_I(x)  ..  - yyix),  one  may  obtain  from  (23) 


Mx)  --T  .y'^ 

~ dx 


- v-I 


(25) 


= r‘Tv(e-j;*[neX*i-, 

= r'TyT-J(T(e-X-j;neV--,} 
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where  y diagonal  complex- valued  matrix  with 

i^iii  = n 


Mij=“  • 


(26a) 

(26b) 


One  can  easily  show  from  (21)  the  identity  Y"^  Ty  = ZTy“^  which  is  used  in 
(25). 

The  solution  in  the  time  domain  can  be  found  since  V(x,  t)  = V(x)e'^‘"‘' 
and^(x,  t)  = I(x)e'^'‘'*'  by  multiplying  (23)  and  (25)  by  e . It  should  then  be 
clear  that  the  total  solutions  consist  of  forward-traveling  waves,  |£^(x,  t), 
J)^(x,  t),  and  backward-traveling  waves,  (x,  t),  (x,  t),  on  the  line  with 

[18] 


where 


^(x,  t)  = 2^^{x,  t)  + 3^'  (x,  t) 

(27a) 

^(x,  t)  = ^^(x,  t)  - (x,  t) 

(27b) 

^+(x,t)  = T e'^^  I 

/S/  «— 

(28a) 

^'(x,  t)  = T eX*  I - 

ew  /s/ 

(28b) 

y-^(x.t)  = z^3  jV,t) 

(28c) 

•^“(x,  t)  = Z^^-(x,  t) 

(28d) 

characteristic-impedance  matrix"  relating  the 

voltages  and 

currents  in  the  waves  with  Z/-  defined  from  (23),  (25)  and  (28)  as 

Zc  = Y’l  T y T-1  = Z T Y*^ 

<-o  es/  /s/  ^ /Si/  rw»  /s/ 

Z^  = Y"WyZ  = Z(/yZ)-^ 

V /s/  V /-./  /s/  ^ W ^ /S/  • 


(29a) 

(29b) 
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The  symbolic  notation  in  (29b)  conforms  to  the  scalar  characteristic  imped- 
ance for  two- conductor  lines  discussed  above.  It  can  be  shown  that  [l8] 

npi  = T Y T“^  (30a) 

= Y"5/y  Z Y (30b) 

/Tz  = Yj^T  Y“1  , (30c) 

The  relations  in  (30)  may  be  easily  shown  [18]  by  forming  Z)  ^'y  Z)  = 

(T  Y T-1)  (T  Y T-1)  = TY®  T“1  = YZ  and  {/^)  {/^)  = (Y"VyZ  Y)  (Y“Vy^Y) 
Z y.  Note  tha t /y  Z 4j"L  Y and  the  order  of  multiplication  of  the  matrices 
cannot  be  interchanged  since  Z and  Y do  nothin  general,  commute. 

If  the  mode  currents,  l^(x),  are  defined  as  in  (23)  and  the  mode  vol- 
tages are  defined  from  (25)  as  ^x)  = Zq  TVjj^(x),  then  it  is  clear  that  the 
mode  quantities  consist  of  n uncoupled  waves  and  each  mode  has  the  propa- 
gation constant  Yj^.  The  velocities,  V£,  and  attenuation  constants,  r|£,  asso- 
ciated with  each  mode  are  found  b^>-  writing  = r|£  + j(u)/v£)  where  t|£  and  v^ 
are  real  scalars.  Thus  one  might  think  of  these  "mode”  quantities  as  being 
somewhat  basic  quantities  in  the  overall  propagation  of  the  waves  since  the 
total  voltages  and  currents  are  linear  combinations  of  the  mode  voltages  and 
mode  currents^ respectively.  This  concept,  however,  is  not  particularly 
useful  in  obtaining  numerical  solutions  to  a given  problem  via  machine  com- 
putation and  is  only  offered  as  a link  to  the  more  familiar  two- conductor 
case  discussed  above*  There  are,  however,  instances  where  this  concept, 
when  related  to  matrix  scattering  parameters,  can  prove  useful  in  certain 
synthesis  problems  [19]. 


The  state  transition  matrix  or  chain  parameter  matrix,  4(x,x-.),  in  (8) 

U 

which  relates  voltages  and  currents  at  the  two  ends  of  a section  of  the  line 
extending  from  Xq  to  x can  also  be  obtained  by  eliminating  and_^"  from 
(23)  and  (25)  as  Fl8,26,47] 


y(x)" 

i(x)^ 


= 1(x,Xq) 


^(xq)! 

(31) 


(32a) 

(32b) 


(32c) 


where  the  nyn  subnaatrices,  $^,.(x,  Xq),  i,  j = 1,  2 are  given  by  [26,47] 

= 1/2  Y-1  T +*  -V(x-xo)^  ^.1^ 

=-1/2  Y-1  T V T-l/TleJ:'*-**)’  - e-X<*-*0>)  t'*  ] 

*„(x.x.)  = .1/2  T (J  Y-‘  T-1  Y 

= -l/2/T(e~^*"*0^  T"^Jt  y“^  T"^  Y 

I ^ ^ 'Ni'  \ fS^  ^ 

}„(x,x.)  = 1/2  T(eX<'‘-=S))  + ,-^(*-*0')  T-‘  . 

From  (21),  one  can  obtain  Y"^Tv  = ZTY'^and  therefore  Y"^Ty  in  (32)  can 
be  written  in  terms  of  Z,  ' 

The  state  transition  matrix  can  also  be  obtained  as  an  absolutely  con- 
vergent matrix  infinite  series  [38,42] 

t(x,X(,)  = jW(*-*o)  ^ ^ + — (33a) 


(32d) 


where  from  (4c) 


1! 

r 


M = 


0 

n~n 


2! 
-Z 


3! 


-Y  0 

~ n~nj 


(33b) 


After  obtaining  the  indicated  products  of  ^ one  can  obtain  [18] 
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+ _^t(Zp=  +...  ,54„ 

= Y-‘fj  ^YzilJ^%(YZ).ir3>\...lY 

~ I n 2!  4!  / i 

♦.,(».*(,>  =-2  -—  <34b) 

- -S(j^)-‘{/yz  (x-x„)  + 

-..j 

= -fyix  <-«o>  ^ 

4 (y^,»  4 -.-]  ^ Y-‘ 

jl,,(x.x„)  =-Y(x.X(,)-  YZY  (YZ)»  Y<!1^  ....  (34c 

= (x.x„,  *(757)3 

^ (yiX)'  ^ 4 ...j 

=-fyiI<>‘-^>4(yYZ)3 

4<yvi.3<i^*-}(/YZ,-lY 

= In  4 X Z + (Y  Z)>  + ...  (34d) 


(34d) 


= yfi  + ZY  + (ZY)»  ^*'*o)* +.__■)  v-i 

~(~n ZT  — 51 — J ~ . 


Matrix  hyperbolic  functions  Cosh  and  Sinh  ixiay  logically  be  defined  in 
the  foUowing  manner.  The  matrix  exponential.  logically 


bs  defined  as  the  absolutely  convergent  matrix  infinite  series  [l8,  38,  42] 


+ (yy|)=  + (yYZ)'>il^%  — j35> 

The  matrix  exponential,  e^^*  can  similarly  be  defined  as  an  absolutely 

convergent  nnatrix  infinite  series  [16,  38,  42] 

^v(x-xo)  ^ ^ ^ , ...  (3„ 

M X ! 2 J ^ 3 f 0 

Since  Y Z is  assumed  to  be  diagonalized  by  T as  in  (21),  then  the  square 

root  of  Y Z nnay  be  defined  as  /Y  Z = T y T"^  as  shown  in  (30).  Therefore, 

(35)  may  be  written  as  [18,  38,  42] 

^fTz  (x-xo),TeX‘*-*o'T-l 


(37) 


whereVz  Y is  defined  in  (30b)  and  (30c),  Thus  the  nnatrix  hyperbolic  func- 
tions  Cosh  and  Sinh  may  be  defined  from  (35),  (36)  and  (37)  as 

C<.,hfyYl  (x.Xo)}  = w[i!n.  '*-*0’  + e-'/n  (38x) 

= 1 + {/yz)^  + (JTz)* 

~n  2!  — 5T 

= l/2TfeX<"-"0)  + ,-X(--0)]T-l 


Sinh{yYZ  (x-Xq)}  = l/2[e' 


^(x-Xq)^ 


(x-Xo)J 


(38b) 


= /7z  + (/YZ)= 

+ ( /y~Z)^  (x-Xq)^^.  ___ 
~ ~ 5 ! 

= l/2T[e^^^-*0).  , 
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Therefore,  it  should  be  clear  by  utilizing  the  relations  in  (35),  (36),  (37) 
and  (38)  that  the  expressions  for  the  chain  parameter  submatrices  In  (32)  and 


(34)  are  equivalent  and  may  be  written  symbolically  as  [18] 

J,,(x,Xo>  = Cosh|'yz7  (x-Xg)}  = y-1  Cosh  JTz  (x-Xg)j  Y (39i) 
Sinh|'yY^(x.xo)^  = -Sinh  (x-x^jj/zT^^Y-J 

= -ZcSinh[yYZ(x-Xo)}  = -Sinh  [/iT  (x-x„ ) } Z^. 

^.(x,Xo)  = -Y(yfY)-‘sinh^yZY(x-Xo)|  = -Sinh|'yYZ  (x-x„)]  (/yZ)-'  Y 
= -Z'‘  Sinh  [,/^  (x-x^)]  = -Sinh  [ (x-x^)}  2.'^ 

(x,  Xq)  = Cosh  fyn  (x-Xq)}  = Y Cosh  {m  <39d) 

where  the  characteristic-impedance  matrix,  Z^,  is  defined  in  (29)  and  (30). 
(Note  that  these  reduce  to  the  scalar  elements  for  two- conductor  (n=l)  lines 
in  (17),)  For  numerical  machine  computation,  however,  one  would  use  the 
forms  of  the  submatrices  given  in  (32)  since  the  equivalent  expressions  in 
(34)  and  (39)  would  be  of  little  practical  value  in  obtaining  numerical  results. 
Also  one  can  show  certain  fundamental  matrix  identities  involving  the 


submatrices  of  ine  chain  parameter  matrix  [18]: 

Identity  1;  - «is»ai=i^n 

Identity  2;  i 

Identity  3:  fgg  $~3  = ^ 


''V/ 


Identity  4;  $g  $3  = 

^ ^ r>^ 

Identib/  5:  = $33 


(40a) 

(40b) 


(40c) 


(40  d) 
(40e) 


I 


ll 


where  f • refers  to  $;i(x,x^).  Identities  land  2 reduce,  in  the  case  of  two- 
conductor  lines  (n=l)  where  ’the  submatrices  become  scalars,  to  familiar 
results,  described  above  i.  e.,  i»a  ^ai  1 thi®  determinant  of 

the  chain  parameter  matrix  is  equal  to  unity.  Similarly,  Identities  3,  4 and 
5 also  reduce,  in  the  case  of  two- conductor  lines,  to  familiar  results,  i.  e. , 
These  identities  maybe  proven  by  substituting  the  forms  of  the 
submatrices  given  in  (32)  and  utilizing  the  fact  that  Y»  ®^and  e 

are  diagonal  matrices  whose  products  may  therefore  be  interchanged.  The 
identities  may  be  more  directly  shown,  even  when  Y Z is  not  diagonalizable 
by  a similarity  transformation,  by  recalling  that  the  inverse  of  the  chain 
parameter  matrix  or  state  transition  matrix  is  given  by  5“^(x,Xq)  = $(xq,x) 
[42],  Forming  this  relation  as  jf(x,XQ)  $(Xq,x)  yields  in  partitioned 

form  [18] 

= _^n  n5nl 

$8i(Xo»x)  $gg(xo,x)  ^0^ 

Multiplying  this  result  out  and  observing  from  (34)  that  ^h(x,Xq)  = ^h(xq,x), 
$ja(x,x^)  = ,(Xq,x),  £2j(x,Xq)  = -$3^(Xq,x),  $aa(x,XQ)  = $sa(xQ,x)  yields 

Identities  1,  2,  3,  4 directly  [18].  Identity  5 is  easily  shown  from  (34a)  and 
(34d)  since  Y and  Z are  symmetric,  i.e.,  Z = Z*"  and  Y = Y^,  and  the  trans- 
pose of  a sum  of  matrices  is  equal  to  the  sum  of  their  transposes  [18,  38,42], 
This  also  shows  from  (34b)  and  (34c)  that  and  are  symmetric,  i.e., 

^ X3  ~ s ^81  ~ lai  [18], 


' \ 
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Thus  the  result  conforms  (symbolically)  to  the  two-conductor  case  in 
which  Y and  Z are  complex  scalars  instead  of  matrices.  This  use  of  sym- 


bolic notation  for  the  square  root  of  a matrix  and  the  matrix  hyperbolic  func- 
tions Cosh  and  Sinh  of  course  rnakes  sense  because  it  was  assumed  that  the 
matrix  product  Y Z was  diagonalizable  by  the  similarity  transformation,  T, 

It  is  not  necessarily  true  that  the  matrix  product  Y Z (and  also  Z Y)  will  be 
diagonalizable  by  a similarity  transformation  [41,  4.2],  If  the  product  is  not 
diagonalizable,  then  a similarity  transformation  may  be  found  to  place  Y Z in 
the  Jordan  Canonical  form  and  this  result  is  found  in  [47]  although  numerical 
results  become  more  complicated  to  obtain. 

Thus  one  of  the  important  simplifying  assumptions  is  that  Y Z is  diago- 
nalizable by  a similarity  transformation  as  in  (21).  It  is  often  assumed  that 
YZ  can  be  diagonalized  by  a similarity  transformation  regardless  of  the 
numerical  entries  in  Y and  Z and  this  is, of  course, not  necessarily  true 


[41,  42],  To  more  completely  investigate  the  problem,  determine  the  eigen- 
values of  Y Z as  roots  of  the  n-th  order  complex  polynomial  in  V®  [18,  41,  42] 

det  - YZ  ) =0  (42) 


where  det  denotes  the  determinant.  If  the  resulting  eigenvalues,  Y?,  are 

distinct,  then  diagonalization  of  Y Z is  assured  and  the  n y 1 columns  of  T = 

[T,,  Ty, , T ],  T.,  are  eigenvectors  of  Y Z sati  fying 

(y?  1 -Yz)  T = 0,  (43) 

~n  ~ ~ — i n-1  ' 

for  i = 1,  , n ri8,  41,  42],  But^-vf  course, one  does  not  generally  know  3^ 

priori  if  the  eigenvalues  will  be  distinct  and  considerable  computation  rnay 
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be  required  to  determine  this.  If  there  exist  repeated  eigenvalues,  then  it 
may  or  may  not  be  possible  to  determine  n linearly  independent  eigenvectors 
via  (43).  If  n linearly  independent  eigenvectors  can  be  found,  then  diagona- 
lization  is  assured  [l8,  41,  42].  It  can  be  shown  that  the  eigenvalues  of  Y Z 
are  the  same  as  the  eigenvalues  of  Z Y (see  [42],  pp.  101-102).  When  either 
Y or  Z are  nonsingular,  this  can  be  easily  shown  by  forming  [l8] 
det  (y2  1 -Yz)  = det  (y  (y®  1 - Z y}  Y-l)  = det  (z'^  (v®  1 -Zy}z)  = 

det  1^  -Z  Y^  since  the  determinant  of  a product  of  square  matrices  is 

equal  to  the  product  of  their  determinants  and  det  (Y)  det^"^^  = 

det  (z~^  det  (Z)  = 1.  Also  one  can  form  (43)  as  Y(y?  1 - Z Y)  T.  = 

0i  and  Z"^  (y?  1 -Z  Y)(ZT.)  = 0,  so  that  if  Y Is  nonsingular  then  each  of 

1V--1  i ^ n ^ ^ ^ • I,  ^ 

the  eigenvectors  of  Z Yis  equal  to  the  product  of  Y“^  and  each  of  the  eigen- 
vectors Cl  Y Z (within  a scalar  constant),  and  if  Z is  nonsingular,  then  each 
of  the  eigenvectors  of  Z Y is  equal  to  the  product  of  Z and  each  of  the  eigen- 
vectors of.  Y Z (within  a scalar  constant)  (18].  These  facts  can  be  used  to 
form  the  relations  in  (23),  (25)  and  (32)  in  terms  of/ZY  and  its  eigenvectors. 

When  discussing  the  question  of  distinct  eigenvalues  in  numerical  com- 
putation, it  is  important  to  consider  the  question  of  "how  distinct".  For 
example,  if  two  of  the  eigenvalues  are  distinct  only  after  the  16-th  digit, 
then  although  they  are  technically  distinct,  the  two  eigenvectors  from  (43) 
associated  with  these  two  "almost-distinct"  eigenvalues  may  be  very  nearly 
collinear  causing  T to  be  an  ill-conditioned  matrix  with  a very  small  deter- 
minant, i.  e,  , T will  be  "almost  singular".  Thus  numerical  instabilities 
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and  other  associated  errors  can  occur  when,  for  example,  computing  the 


inverse  of  T,  T"^,  since  T may  be  an  ill-conditioned  matrix  having  a very 
small  determinant  [49]. 

This  is  one  of  the  reasons  why  determining  numerically  stable  similarity 
transforrnations  such  as  orthogonal  or  unitary  transformations  are  important 
in  numerical  machine  computations[49].  For  example,  a real,  orthogonal 
similarity  transformation,  T,  can  always  be  found  which  will  diagonalize  a 
real,  symmetric  matrix  and  T"^  = where  the  transpose  of  a nntrix  M is 
denoted  by  f4i,49].  Also,  complex,unitary  transformations,  T,  can 
always  be  found  which  diagonalize  complex  matrices  which  are  either  hermi- 
tian  or  norrnal  and  T“^  = where  the  complex  conjugate  transpose  of  a 


matrix  M is  denoted  by  M*  [41,49].  Hermitian  matrices  satisfy  M = and 
normal  matrices  satisfy  (M)(M’i‘)  = [41], 

Machine  computation  of  the  eigenvalues  and  eigenvectors  of  Y Z is  not 
generally  periormed  by  a direct  application  of  (4Z)  and  (43).  Instead  of 
directly  applying  (42)  and  (43),  a more  efficient  method  would  be  to  tranform 
Y Z with  a similarity  transformation  to  some  other  more  convenient  form 


whose  eigenvectors  and  eigenvalues  are  related  to  these  of  Y Z.  For  example 
it  is  known  that  it  is  always  possible  to  obtain  an  n\n  complex, similarity  tran- 
formation,  U,  which  is  unitary  that  will  reduce  any  nyn  complex  matrix  (in 
particular  Y Z)  to  upper  triangular  form,  i.  e.  , YZU  = M and  U*  = U“^, 
where  M has  zeros  below  the  main  diagonal  [41],  Then  since  M is  similar  to 
YZ  (in  the  mathematical  sense  of  similarity),  the  eigenvalues  of  M which  are 


-47- 


the  elements  on  the  main  diagonal  of  Mare  the  same  as  the  eigenvalues  of 

Y Z ’^41,49].  A commonly-used  algorithm  the  QR  transformation  [49].  The 
eigenvectors  of  Y Z,  T-,  are  related  to  the  eigenvectors  of  M,  S.,  by  T-  = 

U where  is  an  n\l  eigenvector  of  M associated  with  the  eigenvalue  and 
corresponds  to  the  eigenvector  T-  associated  with  eigenvalue  v?[41, 49].  The 
transformation  to  Hessenberg  form  is  also  commonly  employed  [49]. 

In  addition  to  the  question  of  the  existence  of  a numerically  stable  simi- 
larity traiipformation  which  diagonalizes  the  matrix  product  Y Z,  there  is 
the  problem  of  recomputing  the  eigenvalues  and  eigenvectors  at  each  fre- 
quency being  considered.  Since  the  matrix  product  Y Z is  a function  of 
frequency,  then  one  is,  in  general,  required  to  repeat  the  determination  of 
the  eigenvalues  and  eigenvectors  of  this  complex- valued  matrix  product, 

Y Z,  at  each  frequency  and  this  can  be  a very  time-consuming  task  when  the 
response  at  a large  number  of  frequencies  is  desired.  There  are,  however, 
certain  practical  cases  where  YZ  can  be  diagonalized  by  a numerically 
stable  transformation  and,  moreover,  for  these  cases,  T is  independent  of 
frequency  and  need  only  be  computed  once.  These  important  cases  will  now 
be  discussed. 

3. 1 Transmission  Lines  in  a Homogeneous  Medium 

This  section  will  consider  the  {n+l)-conductor  lines  in  a homogeneous 
medium  represented  in  Fig.  2.  Although  the  lines  in  Fig.  2a  and  Fig.  2b 
can  only  logically  be  considered  immersed  in  free  space  which  is  considered 
lossless,  the  formulation  which  will  be  investigated  will  assume  losses  in 
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the  medium  in  order  that  the  situation  in  Fig.  2c  may  be  considered.  The 
following  important  relations  which  are  well  known  in  the  case  of  two- con- 
ductor lines  in  a homogeneous  medium  are  shown  in  Appendix  A for  the 
case  of  (n+l)-conductor  lines  in  a homogeneous  medium  which  is  assumed  to 
be  characterized  by  u»  Oi 

LC  = CL  = Uei^  (44a) 

LG  = GL  = Ual  . (44b) 

When  the  dielectric  medium  is  lossy  as  in  Fig.  2c,  the  conductivity  in  (44b) 
refers  to  the  effective  conductivity  as  a=  + me"  = uue^  e^,  tan  6 and  includes 
the  combined  losses  due  to  ohmic  conductivity,  a^,  and  dipole  relaxation 
effects.  The  loss  tangent  of  the  medium  is  denoted  by  tan  5,  e^.  is  ihe  per- 
nruttivity  of  free  space  and  is  the  relative  dielectric  constant.  The  per- 
mittivity, c , refers  to  the  real  part  of  the  complex  effective  permittivity, 
i.  e.,  €=  Cy  Cy,  and  the  permeability,  U,  will  typically  be  that  of  free  space, 
Uy.  In  addition,  since  the  medium  is  homogeneous  it  can  al'^o  be  shown  [54] 
that  C = e K and  from  (44)  it  follows  that 


c = cK 

(45a) 

L = UK-^ 

(45b) 

G = OK 

(45c) 

where  K is  an  n\n  real,  symmetric,  positive  definite  matrix  independent  of 
c(and  therefore  frequency)  and  is  dependent  only  upon  the  cross-sectional 
structure  of  the  line  (conductor  separations  and  wire  radii).  The  matrix 
product  Y Z with  the  relations  in  (44)  and  (45)  becomes 
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(46) 


Y Z = (a  + j(D  e)  K (R^  + + (j'juucr  - uu^uc) 


a*.d  if  perfect  conductors  are  assumed,  then  all  n mode  velocities  and 
attenuation  constants  degenerate  into  one  set,  which  represents  the  true 
TEM  mode  of  propagation. 

If  perfect  conductors  are  assumed,  i.  e. , R = L - 0 , then  from  (46) 

~c  ~c  n-n  ' 

T = 1^  and  y?  = (jit'ua-  (jj®  uc)  in  (?1)  where  ^ is  the  nyn  identity  matrix 
with  ones  on  the  main  diagonal  and  zeros  elsewhere.  Thus,  the  matrix 
chaixA  parameters  for  the  homogeneous-medium  case  with  (n+1)  perfect  con- 
ductors become  from  (32) 


$ii(x,Xo)  = cosh  {y(x-Xq)] 

[in 

(47a) 

= - sinh  jy(x-XQ)j-  [(juj/y)  L] 

(47b) 

=-  sinh  {v(x-xq)}  [(juu/ y) 

(47c) 

«23(x,Xq)  = cosh  {y(x-Xq)j 

[in 

(47d) 

where  Y =yj(DU (cj  + juue)  and  the  characteristic  impedance  matrix  becomes 
from  (29a) 

joiu 

(a  + j(j)c) 


K 


-1 


(48) 


= (juu/  v)  L 


For  a lossless  medium,  a = 0,  y-  and  (47)  becomes  [26] 

«ji(x,Xq)  = cos  j 

(49a) 

$jg(x,Xo)  = - j sin  {3(x-Xq)}[uL] 

(49b) 

$3i(x,Xq)  = - j sin  {g(x-XQ)]-[u  L]-^ 

(49c) 

^aa(x,  Xq)  = cos  ■ 

[6(x-Xo)}j.n 

(49d) 
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\y:- 

I 


p:;' 


p.  I.  Ill) 1 1. 


v^ere  the  wave  number,  fi,  is  given  by  3 = 2tt/X,  \=  u/f,  u=  l/J^  and  the 
characteristic-impedance  matrix  is  real  and  becomes  Zf,  ~ n L. 

If  perfect  conductors  cannot  be  assumed,  then  from  (46)  it  is  sufficient 
to  find  a T which  diagonalizes  K (R^  + jmL  ),  i,  e. , 

^ eo  /w  C 

(50' 

where  A®((ju)  is  an  nyn  diagonal  matrix  with  [a®(u))].  . = A®  ({«)  and  [A^(m)J..  = 0 

^ ^ /w  II  1 ^ 1.J 

for  i ^ j.  The  eigenvalues  can  then  be  found  from  (46)  and  (50)  as 

y\  = (a  + ja>e)  a|  (ot)  +(juiMa-  o)®ue)  . (51) 

In  general,  diagonalization  as  in  (50)  is  not  assured  since  K (R  + jmL  ) is  a 
complex  matrix  with  no  particular  structural  properties  which  would  be 
useful  in  determining  a-priori  whether  the  matrix  is  diagonalizable,  i«e., 
hermitian  or  normal. 

If  one  neglects  the  internal  inductance  of  the  conductors,  i.e, , L.  = 

„0_,  or  neglects  the  resistance  of  the  conductors,  i,  e. , R„  = 0 , then 

n~n  •—  -^c  n~n 

numerically  stable  transformations  can  be  found  which  diagonalize  each  of 
these  cases  but  not  both,  i.e.,  there  exists  a T such  that  [l3j 

T-1  K R T = a®(u))-^{l^  = _0  J (52) 

or  there  exists  a T such  that 

j.i,  T- 1 K L^  T = a®  (u))  — {r,  = „ 0^}  (53) 

but  the  same  T will  not  necessarily  simultaneously  diagonalize  both.  That 


this  can  be  done  relies  only  on  the  fact  that  K is  real,  symmetric,  positive 
definite  and  that  and  L^  are  real  symmetric  [13,  41,  42],  The  construc- 
tion of  a numerically  stable  transformation,  T,  which  will  diagonalize  the 
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product  of  a real,  symmetric,  positive  definite  matrix  and  a real,  symmet- 
ric matrix  will  be  shown  in  Section  3,2  and  may  be  corr.'puted  very  effi- 
ciently with  the  subroutine  NROOT  in  the  IBM  Scientific  Subroutine  Package 
(SSP)  [50].  Generally  for  high  frequencies,  the  entries  in  are  much  less 
than  the  corresponding  entries  in  Land  the  approximation  in  (52)  would  be 
relatively  accurate  [l3].  However,  in  either  case,  since  both  and  L^ 
are  functions  of  frequency,  the  transformation  matrix,  T,  and  the  eigen- 
values  must  be  recomputed  at  each  frequency  under  consideration  and  this 
increases  the  overall  computation  time. 

There  are  cases  where  one  can  include  both  resistance  and  internal 
inductance  of  the  conductors  and  obtain  a numerically  stable,  frequency- 
independent  transformation.  For  example,  consider  Fig,  2a  in  which  all 
(n+1)  wires  are  assumed  to  be  identical.  In  this  case,  (50)  becomes  (see  (6) 
and  (7)) 

(rc  + T-‘  K {i,„  + uj  T = (M) 

where  the  (n+1)  conductors  (including  the  reference  wire)  have  resistance, 
r^,  and  internal  inductance,  and  Uj^  is  the  n^n  unit  matrix  with  one's  in 

every  position,  i.e.,  j = » Note  that  even  though  K and 

"(in  Un}  each  are  symmetric,  it  is  not  necessarily  true  that  their  product 
will  be  symmetric.  Since  K is  real,  symmetric  and  positive  definite  and 
{~n  ^ IJn^  real,  symmetric  then,  as  discussed  before,  the  product  can  be 
diagonalized  and  NROOT  in  SSP  can  be  used  to  perform  the  reduction  [50]. 
Furthermore,  T will  be  independent  of  frequency  and  need  be  computed 


only  once  in  the  frequency  response  solution  and  the  eigenvalues  can  be  re- 
computed very  simply  at  each  frequency  from  (51),  Assuming  that  the  n 
wires  in  Fig.  2b  and  Fig.  2c  are  identical,  then  this  technique  applies  since 
U does  not  appear  in  (54)  because  the  ground  plane  and  circular  shield  are 
assumed  to  be  perfect  conductors.  In  this  case  one  only  needs  to  diagonsl- 
ize  K which  can  be  accomplished  with  the  subroutine  EIGEN  in  SSP  f50]  since 
K is  real,  symmetric. 

3.  2 Transmission  Lines  in  Inhomogeneous  Media 

One  of  the  main  problems  under  consideration  in  this  report  is  the  case 
of  circular  wires  with  circular. dielectric  insulation  as  shown  in  Fig.  3 which 
appear  in  the  form  of  bundles  of  closely  coupled. dielectric- coated  wires. 
These  commonly  occur  in  electronic  systems  in  the  form  of  densely  packed 
cable  bundles  and  flat  pack  or  woven  cables  [51].  The  inhomogeneity  in  the 
surrounding  medium  (free  space  and  insulation  dielectric)  makes  the  identi- 
ties in  (44)  no  longer  true.  However,  it  is  always  possible  to  diagonalize 
the  matrix  product  Y Z with  a numerically  stable  transformation,  T,  when 
perfect  conductors  and  dielectrics  are  considered  regardless  of  the  entries 
in  C and  L. 

First  consider  the  case  where  losses  are  neglected,  i.  e. , G = R_  = L = 
The  matrix  product  becomes 

Y Z = - u)®  C L . (55) 

Recall  that  L and  C will  be  real,  symmetric  and  C will  be  positive  definite 
even  for  this  inhomogeneous  medium  case  [39].  Since  C is  real,  symmetric; 
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then  there  always  exists  an  nyn  real,  orthogonal  transformation  U such 

that  . 

U"‘  C U = D (56) 

where  D is  an  nvn  real,  diagonal  matrix  and  U"^  = f41,  49].  Further- 

more,  since  C is  positive  definite,  the  eigenvalues  of  C which  are  the  ele- 

ments  of  the  diagonal  matrix  D are  all  positive,  real  and  nonzero.  Thus 

one  can  quite  easily  (and  meaningfully)  form  the  square  root  of  the  matrix 

1/2 

D,  D , and  write 

p-1/2  y-1  ^ jj  d»1/2  pl/2  y-1  y j^l/2  ^ ^1/2  L u (57) 


which  is  real,  symmetric.  Thus  (57)  may  be  diagonalized  again  by  an  n\n 
real,  orthogonal  transformation,  S,  such  that 

gt  d1/2  u‘  L U s = a*  (58) 


and  one  can  identify  the  transformation  matrix  T in  (21)  as 

T = U S 


(59) 


and  propagation  matrix  in  (21)  becomes 

y2=_U)2A®.  (60) 

The  propagation  constants  become  from  (60),  = ja>A^  where  [a®]^  = A^^, 


[a®]..  =0,  i j and  it  is  a simple  matter  to  verify  that 
~ ij 

T"  ^ = t‘  C"  ^ 


(61) 


The  matrix  chain  jjarameters  for  this  case  are  given  in  (32)  and  [26]  and  the 
subroutine  NROOT  in  SSP  will  again  perform  this  type  of  reduction  [50].  If 
the  real  parts  of  the  permittivities  of  the  insulations  are  independent  of  fre- 
quency (or  assumed  to  be)  then  this  reduction  need  be  performed  only  once 
and  if  the  real  parts  of  the  permittivities  vary  significantly  with  frequency. 
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one  must  recompute  T and  v (as  well  as  C)  at  each  frequency.  In  either 
case,  T will  be  real- valued  and  numerically  stable. 

In  the  general  case,  the  naatrix  product  Y Z becomes 

YZ=(G+jtt)C}  (R.  + + (G  + jujC)  (jo)L)  . (62) 

Even  if  perfect  conductors  are  assumed,  i.e.,  = L = 0 , diagonaliza- 

C C A 

tion  of  Y Z would  require  the  diagonalization  of  the  complex  matrix  jujGL  - 
tjj®  CL.  However,  G in  general  bears  no  simple  relationship  to  L or  C such 
as  in  (45)  since  the  fields  associated  with  conduction  current  or  dipole 
relaxation  losses  will  be  confined  to  the  insulation  dielectrics  whereas  the 
fields  associated  with  the  real  parts  of  the  corr^lex,effective  permittivities 
of  the  dielectrics  can  fringe  into  the  surrounding  free  space  medium.  Thus 
the  diagonaliza  tion  of  Y Z is  not  assured  a priori.  If  diagonaliza  tion  is 
possible,  T would  in  general  be  complex  and  a function  of  frequency. 

If  the  dielectrics  are  assumed  to  be  perfect  (no  ohmic  conductivity  or 
dipole  relaxation  effects),  then  assuming  all  n conductors  are  identical 
(including  the  reference  conductor  in  Fig.  3a)  Y Z becomes  for  Fig.  3a 

YZ  = j,i,(r_  + joof^)  C (1„  + U_)  - L . (63) 

For  a real,  frequency  independent  transformation,  T,  which  diagonalizes 
Y Z to  exist,  it  would  be  required  in  general  that  the  same  T diagonalize 
both  C (1„  + U_)  and  C L.  This  is,  in  general,  not  possible.  Even  if  the 
reference  conductor  is  assumed  lossless,  i.e.,  U_  = _0_  in  (63)  for  Fig.  3b 
and  Fig.  3c,  the  existence  of  a real,  frequency-independent  transformation 
which  diagonalizes  Y Z would  imply 
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(64a) 


T-1  C T = A? 

T-^  C L T = Af  T-1  L T = A®  (64b) 

where  A®  and  ^ are  nyn  diagonal  matriceR.  This  would  therefore  imply 
that  the  same  T would  diagonalize  both  C and  L and  this  is  generally  not 
possible. 

Therefore,  the  inclusion  of  losses  generally  requires  that  a complex 
transformation  T be  obtained.  The  existence  of  a numerically  stable  trans- 
formation is  not  guaranteed,  in  general,  when  losses  are  included.  T is 
also  a function  of  frequency  which  requires  that  it  be  recomputed  at  each 
frequency  which  increases  the  overall  computation  time. 

3,  3 Cyclic-Symmetric  Matrices 

If  the  n conductors  and  dielectric  insulations  are  identical  and  if  the 
cross-sectional  structure  of  the  line  exhibits  certain  physical  symmetry 
with  respect  to  each  of  the  n conductors  and  the  reference  conductor,  then 


the  matrix  product  Y Z can  be  diagonalized  ^ priori  with  a transformation 
matrix,  T,  which  although  complex,  is  independent  of  frequency  even  when 
lossy  conductors  and  lossy,  inhomogeneous  media  are  considered.  For 
example,  if  the  n conductors  are  identical  with  identical  dielectrics  all  of 
the  same  thickness,  and  are  equally  spaced  with  respect  to  each  other,  on 
a ring  symmetrical  about  the  reference  wire  or  are  equally  spaced  with 
respect  to  each  other  on  a ring  concentric  with  the  circular-shield  reference 
conductor  as  shown  in  F'ig.  8,  then  Y Z is  always  diagonalizable  by  a fre- 
quency independent  tranuformatio'.  T. 
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For  each  of  the  lines  in  Fig.  8,  Z and  Y will  quite  obviously  be  of  the 


form 


Z = 


Zi  Z,  Z,  • • • • Z3  z; 

Yi  YaYg*  * • • Y3Y3' 

^2  ^3 

Ya  Y,  Yg  Y3  Y3 

• , 

* • 

• 

• * . 

Zg  Zj  • . " 

Y3  Ya  . . . 

• * • • . : 

Y = 

• • • * * • 

«> 

• • • • #7 

. • . - 

• . • ‘ . Y3 

« « 

* • « * 

• • 

« 

Y3  • Y,  Y, 

• 

Z2Z9 ZgZaZ^ 

Y2Y3 Y*3Y,Y,_ 

(65) 


where  fZ]..  A 2,,  fY]..  A Y, , [Z].  [Y]..  are  defined  in  (6)  and  (7). 

~ u = ^ ~ u = ^ ~ ij 

Matrices  with  this  special  structure  are  cyciic-symmeti  ic  matrices  and  the 

general  nyn  cyclic-symmetric  matrix,  M,  is  defined  by  [M]..  = M. 

~ ~ [i- j 1 +1 

v/here  M.  . = M.,  M . = M.  and  indices  greater  than  n and  less  than  1 
j+n  j’  n+2-j  j ^ 

are  defined  by  the;  convention  n+j  = j and  n+i  = i [52,  6’  ' Because  of  the 
special  structure  of  the  matrices,  there  al’ways  exists  a transformation,  T, 
which  is  independent  of  frequency  and  the  numerical  entries  in  Z and  Y 

rw' 

which  will  diagonalize  both  Z and  Y,  i.  e. , T~^  Y T = y®  » T"  ^ Z T = y®  and 
y^  = y%  y^  where  y%  and  y^  are  nyn  diagonal  matrices  [l,  b,  52,  65],  The 
elements  of  T which  diagonalize  any  cyclic-symmetric  matrix  of  the  form  in 
(65)  are  F5,  52,  65] 


[T]..= 

~ u /IT 


(66) 


where  a complex  number  c with  magnitude  c^^  and  angle  0^  is  written  as 
T is  unitary  such  that  = T*  an  -'yclic-symmetric  matrices 

i Tl  / r>* 
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can  be  shov/n  to  be  normal  matrices,  i.e. , since  it  can  be 

shown  that  the  product  of  any  two  cyclic-symmetric  matrices  of  the  same 
order  commute  under  multiplication  [52],  The  eigenvalues,  of  the  product 
of  two  cyclic-symmetric  matrices  of  the  form  in  (65),  Y Z,  can  be  shown  to 
be  [5,  52] 

'I  = »-»}  <«) 

where  [2]jp  and  are  the  elements  in  tb®  first  row  and  p-th  column  of 

Z and  Y in  (65)  respectively,  p=l, , n. 

Thus  if  the  line  consists  of  n identical  conductors  with  identical  insula- 
tions and  thicknesses  and  exhibits  certain  cross-sectional  symmetry,  then 
the  matrix  product  Y Z can  always  be  diagonalized  regardless  of  the  numeri- 
cal entries  in  Y and  Z and  the  transformation  matrix  is  independent  of  fre- 
quency. Neither  the  transformation  matrix  T,  T"^  nor  the  eigenvalues  need 
be  computed  since  they  are  known  d priori  through  (66),  (67)  and  T“^  = T’!«. 

Cyclic-symmetric  matrices  are  obviously  quite  desirable  fruiii  a com- 
putational standpoint  and  have  been  used  in  modeling  cable  bundles  under 
the  assumption  that  the  conductors  are  arranged  symmetrically  about  the 
axis  of  an  overall  shield  or  occupy  all  possible  positions  within  the  shield 
randomly  [52],  Special  cases  of  cyclic-symmetric  matrices  are  encountered 
throughout  the  power  transmission  literature  under  the  assumption  that  the 
power  line  is  balanced  or  completely  transposed  and  the  transformation 
matrix  is  often  referred  to  as  a symmetrical-component  transformation 
[6,  13]. 
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This  molhofi  cannot  Kcncsrally  ho  applied  when  the  reference  conductor 
is  an  infinite  ground  plane  since  the  special  structures  of  Y and  Z in  (65) 
will  not  result  unless  the  system  of  n wires  and  the  n images  used  to  re- 
place the  grovind  plane  possess  the  required  symmetry.  However,  for  the 
case  of  a three- conductor  line  (n=2)  consisting  of  two  identical  wires  both 
at  the  same  height  above  a ground  plane,  i.e. , r^j  = ^^2’ 

Oj  = Og,  h^  = hg  in  Fig,  2b  and  Fig.  3b,  then  Y and  Z ^^^.ll  be  cyclic-symmet- 
ric regardless  of  the  wire  spacing,  dj  or  any  form  of  transposition.  In 
this  case,  the  elements  of  the  eigenvectors  become  real  as  = n/^ 

Ta  1 = V/T,  Tja  = 1//T,  Tgg  = -1//Tand  the  eigenvalues  are  easily  shown 
to  be  vf  = (Z  + Z^KY  + Y^),  vS  = (Z  - Z^)(Y  . Y^^)  where  Z = 'Zj.,  = 

= tZ]„  = IZ]„  and  Y = fY]„  = [Y)„,  Y^  = [Y],,  = Iy),,. 


IV.  INCORPORATING  THE  TERMINATION  NETWORKS 


Note  that  in  (23)  and  (25),  the  only  unknowns  are  the  2n  undeterndned 
constants  in  I ^ and  I These  will  be  determined  by  the  boundary  conditions 
(termination  networks)  at  x = 0 and  x = X for  a line  of  total  lont^th  J (sec 
Figure  9).  The  incorporation  of  the  termination-networks  can  consume  con- 
siderable computation  time  for  large  numbers  of  mutually  coupled  conductors 
and  this  necessary  step  in  the  total  problem  solution  is  generally  dismissed 
as  a trivial,  straightforward  problem.  It  is  straightforward  (conceptually) 
but  is  certainly  not  trivial  when  a large  number  of  mutually  coupled  con- 
ductors are  involved. 

For  two- conductor  lines,  the  terminations  (which  are  assumed  to  be 
linear)  are  represented  by  Thevenin  equivalents  as  shown  in  Fig.  6,  The 
terminal  equations  become 

V(0)  = Vq  - Z^I(O)  (68a) 

VU)  = V^+  Z^liZ)  (68b) 

where  Vq  and  Vj,  are  equivalent  open-circuit  port  voltages  with  respect  to 
the  reference  conductor. 

For  multiconductor  lines,  the  termination-networks  are  similarly  con- 
sidered to  be  linear  n- ports  and  are  cha racte rizable  by  "Generalized 
Thevenin  Equ’  /alents"  as 

V(0)  = V - Z 1(0)  (69a) 

— —0  ~0~ 

= V^  + Z^±IZ)  (6Qb) 
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where  ^ Xz  complex- valued  vectors  of  equivalent  open-circuit 

port  excitation  voltages  with  respect  to  the  reference  conductor  and  and 
^ are  n«yn  complex- valued,  symmetric  matrices.  The  linear  n ports  can 
quite  obviously  be  characterized  by  (69)  and  shown 

in  Fig,  10,  For  arbitrary  termination-networks,  the  entries  in  (69)  can  be 
quite  easily  obtained  by  treating  V(0)  and  V(Z)  as  independent  voltage  sources 
and  writing  the  loop  current  equations  of  the  networks.  The  currents  1(0)  and 
I(£)  will  comprise  subsets  of  the  loop  currents  for  the  networks  and  the 
remaining  loop  currents  can  be  eliminated  to  yield  (69).  If  the  i- th  conduc- 
tor is  connected  to  the  reference  conductor  only  through  impedances  and 

Z„.,  then  the  entries  in  Z_  and  Z_  are  easily  obtained  as  [Z„]..  = Z^.,  [Z 
£i  ~0  ~uu  0i~UiJ 

■ f^X^ii  ■ ^jji*  ^5s^ij  ,nand  i?^j. 

Combining  (23),  (25)  and  (69)  one  can  obtain  straightforwardly  [26,  48] 


’■■feoi-r'ii} 

feo  I ► r * I 

X" 

lo' 

|z  T + Y"^  T y}e~^ 

- {z.  T - Y"^  T yle'X^ 

L ' J 

m ^ 

Since  T"^  Y Z T = y®,  then  T y in  (70)  can  be  replaced  by  Z T y"^. 


Once  this  set  of  2n  equations  in  the  2n  unknowns,  I ^ and  I are  solved  (by 
Gfiussian  elimination  and  back  substitution,  for  example,  [49])  then  the 
response,  V(x)  and  I(x),  at  any  point  on  the  line  can  be  determined  from  (23) 
and  (25).  For  two-conductor  lines,  the  matrices  and  vectors  in  (70)  become 
scalars  and  T becomes  1 (see  (13)  and  (14)). 

It  is  also  possible  to  indirectly  solve  for  the  response  via  the  matrix 
chain  parameters.  With  x = X and  Xq  = 0 in  (31)  and  (32)  and  using  (69)  one 
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can  straightforwardly  obtain  [26,  48] 

i(0=&x(0  +[«,,(«- 5,,  U)ZjIi(O)  (71b) 

where  0)  ^ V(x)  and_I(x)  can  be  obtained  for  any  x from  (8)  with 

1(0)  from  the  solution  of  (71a)  and  V(0)  determined  from  (69a),  Here  one 
need  only  solve  n equations  in  n unknowns,  equation  (71a),  as  opposed  to  2n 
equations  in  2n  unknowns  in  (70).  However,  certain  matrix  multiplications 


are  required  in  forming  both  (70)  and  (71). 

Using  the  rnatrix  chain  parameter  identities  in  (40),  it  can  be  shown 
that  (71)  may  be  written  in  an  alternate  form  [18] 


in  1 

i(0) 

in  5l- 

I(£) 

Lr*  J 

(72) 


which  has  a highly  sparse  (large  number  of  isero  elements)  coefficient 
matrix  with  2(n®  - n)  of  the  4n*  elements  identically  zero.  Equation  (72) 
can  also  be  solved  explicitly  for_I(0)  and  1(1)  as  [18] 

i(X)--{«,,(a^.  »„(X)}i(o)+i.(x)^  . 

The  advantage  of  the  formulation  in  (73)  as  opposed  to  (71)  is  that  only  two 
of  the  rnatrix  chain  parameters,  and  need  be  determined  in  solving 

for_I(0)  and_I(£)  via  (73)  (see  (32)). 

The  most  efficient  method  of  solving  n linear,  algebraic  equations  in  n 
unknowns  is  Gaussiain  elimination  with  back  substitution  (LU  decomposition) 
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\R^ich  requires  (n^/3  + - n/3)  operations  (multiplications  and  divisions) 

or  on  the  order  of  n®/3  for  large  n [49].  Thus  solving  (70)  instead  of  (71a) 
or  (73a)  requires  on  the  order  of  (2n)®/3  or  8 times  the  number  of  operations 
and  the  complete  solution  of  the  problem  requires,  at  a minimum,  the  solu- 
tion of  n complex  equations  in  n unknowns.  The  impact  of  this  requirement 
on  the  overall  solution  times  for  large  cable  bundles  can  be  illustrated  as 
follows.  The  time  required  to  solve  50  complex  equations  with  a standard 
Gaussian  elimination  subroutine  with  full  pivoting  (DGELG  in  SSP  [50]  which 
was  converted  for  complex  arithmetic)  was  12.6  seconds  on  an  IBM  360/65 
computer.  So  if  it  is  required  to  solve  for  100  frequencies,  then  the  overall 
computation  time  would  be,  at  a minimum,  on  the  order  of  21  minutes.  It 
is  not  uncommon  to  find  100  conductors  in  cable  bundles  on  modern  avionics 
systems  and  since  the  number  of  operations  required  increases  on  the  order 
of  n^,  then  solutions  for  100  conductors  and  100  frequencies  'would  require, 
at  a minimum,  2.8  hours!  Of  course,  additional  time  will  be  required  for 
matrix  multiplications  (as  well  as  the  computation  of  Y and  Z and  diagonali- 
zation  of  YZ)  as  indicated  in  (70),  (71),  (72)  and  (73).  This  could  be  quite 
substantial  since  n®  multiplications  are  required  to  multiply  two  "full"  nyn 
matrices  which  is  precisely  the  number  of  operations  required  to  invert  an 
nyn  matrix  'which  is  "full"  [49]. 

Using  the  matrix  chain  parameter  formulation  in  (71)  and  (73)  has  an 
additional  advantage  over  (70).  It  allows  a straightforward  incorporation  of 
incident  electromagnetic  fields  into  the  solution.  Consider  (4)  and  (8)  where 
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the  effects  of  an  incident  field  are  included  as  distributed  sources  along  the 


line,  V (x)  and  I (x).  From  (8)  and  (31),  define  equivalent  sources 

Vg(£)  = (£.x)  VgCx)  + iia(X.x)_Ig(A)}dA  (74a! 

ig(X)  = ^^{lsi(£.x)V^(x)  + $aa(£,x)J[g(x)}dx  . (74b! 

A 

One  can  modify  (71)  and  (73)  to  include  these  sources  by  simply  adding  I,(£) 

A 

to  the  right-hand  side  of  (71b)  and  (73b)  and  replacing  V^,  with  V^-  + 

A 

^_Ig(X)  on  the  right-hand  side  of  (71a)  and  (73a).  This  is  quite  obvious 

A 

since  (8)  shows  that  for  x = £ and  x^  i;  0,  I(£)  is  increased  by  I (£)  and  V{Z) 

g *—  — s — 

A 

is  to  be  increased  by^g(£)  over  the  case  without  incident  field  ilJuniina tion. 
From  the  boundary  conditions  (69b),  MX)  - Zjj_I(£),  then  on  the  right 

A A 

hand  side  of  <71a)  and  (73a)  is  to  be  decreased  by  V„(X)  - Z»  I ,(£)  and  I(£)  is 

A 

to  be  increased  by  I„(X)  in  (71b)  and  (73b),  Thus  equations  (71)  and  (73)  are 

•“S 

quite  easily  modified  to  consider  incident  fields  and  the  final  equations 
become  [26] 

$„(X)  - ^ *s,(I)  Zq  - $,,(£)  + $,,(£)  Zq]_I(0)  = (75a 

- Zj.  5,,(x)] 

UX)  = £p,(X)  + [«aa(X)  - |a^(X)  ^^^(O)  +ig(X)  (75b 

or  [18] 

Ifn  - 5t-  5o  - U0>  = - Is 

- :S) 


* 2 is(X) 


_I(X)  = la^(X)  + [«aa(X)  - $ai(£) 


(76b) 


A more  detailed  discussion  of  efficient  incorporation  of  the  boundary  condi- 
tions is  given  in  [Z6]. 

The  equations  in  (75)  become  particularly  simple  for  the  multiconductor 
line  in  Fig.  2a  consisting  of  (n+1)  perfectly-conducting  wires  in  free  space 
illuminated  by  an  incident  electromagnetic  field.  It  is  shown  in  Appendix 
C that  (75)  reduces  for  the  case  of  Fig,  2a  with  incident  field  illumination 


[co5(8X)  {z„  + Z^  + j»m(Si){Zc  +Z^ZcV3-o})  W)  = 

- Vj.  + [jsin(6£)  Zj.  + cos(p£)^j25j 


(inc) 


(77a) 


+ \ {[cos  (0(£-x))^^  + jsin(0(r-x))  Zjj  Z^^]  _E^(x)]-dx 
0 

(inc)  , (inc). 

- E^(X)  + {[cos(e£)^^  + jsin(0D 

I(X  ) = - j sin(0i:)  Zq  Yq  + [cos(0X)  + j sin(BX)  Z^^  Zq]_1(0/  (7Vb) 

-jZ~c  ^ {sin(0(£- x))  (x)/dx 

-jZ^  |sin(e£)  E^(O)]- 
(inc)  (inc)  (inc) 

where  E.(x),  E (£)  and  E (0)  are  ny  1 column  vectors  with  the  entries  in  the 


i-th  rows  given  by 


(inc)  (inc)  (inc) 

[E^(x)].  = E (u.^,  x)  - E.  (0,x) 


(inc) 


i £.  iO' 


(inc)  d.Q  (inc) 

[E  (£)].  = C E (?  , £)  d^ 

C I \ t;  I I 


(77c) 


(77d) 
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(77e) 


i'^nc)  (inc) 

fE(0)]  = V E (?  . 0)  d? 

-t  I tj  t 

for  i=l, ,n.  1^/*  components  of  the  incident 

M lO  JL£ 

electric  field  intensity  vector  in  the  longitudinal  (x)  direction  along  the 

(inc)  Jtnc) 

axis  of  wire  i and  wire  0 respectively.  The  terms  E,  (5.,  0)  and  E (%.,  D 

I i I ti  1 

are  the  transverse  components  (lying  in  the  y,z  plane)  of  the  incident  elec- 
tric field  intensity  vector  at  x = 0 and  x = X respectively  along  the  contour 
between  wire  i and  wire  0.  The  contour  5.  is  a straight-line  path  between 
wire  i and  wire  0 and  perpendicular  to  these  wires.  The  entries  0 and  Zq 
are  the  wave  number  and  characteristic-impedance  matrix  respectively 
with  3 = 2-^/X,  X = v/1',  V = 1/ and  = vL>.  The  corresponding  solu- 
tion for  Fig,  2b  is  also  discussed  in  Appendix  C. 

4, 1 Lumped-Circuit  Iterative  Approximations 

In  deriving  (4),  "electrically  short"  A*  sections  of  the  line  were  con- 
sidered and  since  the  line  was  assumed  to  be  uniform,  all  sections  will 
be  identical.  Requiring  that  the  &x  sections  be  "electrically  short"  for  all 
frequencies,  the  transmission  line  equations  in  (4)  are  obtained  in  the  limit 
as  Ax  **0.  Alternately,  one  can  construct  lumped-circuit  models  for  the 
line  consisting  of  N identical  sections  of  length  T/N  so  that  each  of  the  sec- 
tions would  in  itself  be  "electrically  short",  e.g.,  the  section  would  be  no 
more  than,  for  example,  I/IO  of  a wavelength  long  for  the  frequency  under 
consideration.  Note  that  since  X/N  is  assumed  to  be  "electrically  short  " 
and  the  cross-sectional  dimensions  of  the  line  are  assumed  to  be  "electri. 
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cally  small"  su  that  transmission  line  theory  applies,  then  a lumped*circuit 
representation  of  this  portion  of  the  line  is  valid.  Perhaps  the  m>re  common 
models  of  the  transmission  line  are  the  lumped-circuit  Iterative  approxima- 
tions which  use  this  philosophy.  Models  are  shown  for  (n4-l)-conductor  lines 
in  Figure  11  as  the  lumped  T model,  lumpedTm  :el,  lumped  Pi  (tt)  model 
and  lumped  Tee  model.  The  lumped*"! section  is  similar  to  the  circuit  in 
Figure  7 but  with  Lx  replaced  by  jVN.  The  lumped  T section  is  the  opposite, 
i.  e. , the  capacitance  and  conductance  elements  appearing  at  the  end  of  the 
lumped*!  section  appear  at  the  beginning  of  the  lumped  V section.  The  lumped 
Pi  section  is  similar  to  the  lumpcd*!section  but  has  half  the  values  of  the 
capacitance  and  conductance  parameters  placed  at  the  beginning  and  at  the 
end  of  the  section.  The  lumped  Tee  section  is  again  similar  to  the  lumped 
section  but  has  half  the  values  of  resistance  and  inductance  (self  and 
mutual)  elements  placed  at  the  beginning  and  at  the  end  of  the  section.  The 
lumped  r model  has  been  used  in  the  program  STRAP  [36];  the  lumped  Pi 
model  has  been  used  in  the  program  IVEMCAP  [34],  and  lumped  Pi  and  Tee 
models  have  been  used  in  power  transmission  line  studies  [9]. 

The  2ivy  2n  chain  parameter  matrix  of  a section  of  line  of  length  X/N 
characterized  by  any  the  lumped  iterative  approximations  can  easily  be 
shown  in  terms  of  Z and  Y to  be  (see  Appendix  D) 

- * ew 

{in}  f-Zli/N)} 

!{.Y(I/N)}  {i„  + YZ(I/N)=}J 


Ik' 


( Lumped*!) 


(78a) 
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(b)  The  lumped  j equivalent  circuit. 


Figure  11.  ijum-Ded-circui.t  iterative  models  of  multiconductor 
transmission  lines  (cont.  ). 


{l^  + Z Y(£/N)2|  {- Z(X/N)} 

[ {-YU/N)}  {!„} 

+ 1/2  Z Y(£/N)=} 

{-Y(£/N)  - 1/4  YZ  Y(£/N)®}  {l^  + 1/2  Y Z(X/N)=} 

{ijj  + 1/2  Z Y(£/N)2}  |-Z(£/N)  - 1/4  Z YZ(£/N)3} 


(Lumped  F ) 


(78b) 


( Lumped  Pi)  (78c) 


{-Y(£/N)j. 


{in  ' X5(^/N)2} 


(Lumped  Tee)  (78d) 


These  models  are  referred  to  as  lumped-circuit  iterative  approximations 

since  the  overall  matrix  chain  pai*ameters  for  the  line  of  length  £ and  N 

N 

sections  (all  of  the  same  type)  is  quite  obviously  ($j^)  , i.e.,  multiply 
the  chain  paraiiieter  matrices  in  (78)  together  N times,  since  the  chain  para- 
meters for  each  section  only  relate  the  voltages  and  currents  at  the  two  ends 
of  each  section  as 


where  each  submatrix  i ^ ^ corresponds  to 

submatrices  in  (78)  and  k = 1,  2, , N. 

Generally  an  N section  lumped-circuit  iterative  approximation  is  solved 
(the  boundary  conditions  or  termination-networks  are  incorporated)  as 
strictly  a lumped,  elec trical  circuits  problem  with  circuit  analysis  ,)rograms 
such  as  ECAP,  SCEPTRE,  TRAFFIC,  etc.  [53],  The  node-voltage  equations 
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\ 


or  loop-current  equations  of  the  network  of  N sections  with  termination - 


networks  are  written  and  solved  via  these  programs.  It  is  possible  to  write 
for  example,  the  node-voltage  equations  directly  from  the  matrix  chain 
p>arameters  without  employing  a circuit  diagram.  Using  (69)  and  (79)»  one 
can  straightforwardly  obtain  (see  Appendix  O) 


where 


X ~ ~ lik'.  1 tkia  tkaa  Ikia) 

^N+1  " {iki2  " tkia  tkaa  tkiaj 


In  deriving  (80)  fiom  (69)  and  (79),  use  is  made  of  the  relation 

ikaa  ikii  - iksi  = Ikls 

which  one  can  readily  verify  from  the  partitioned  forms  of  the  chain 
parameters  in  (78).  In  fact,  this  relationship  can  be  shown  to  be  true  in 


general  for  any  lumped,  linear,  reciprocal  2n-port  characterized  by  the 


chain  parameter  matrix  in  (79)  by  writing  (79)  from  the  nodal-admittance - 
rnatrix  characterizing  the  2n-port  and  invoking  symmetry  of  this  nodal 
admittance  matrix  via  reciprocity. 

In  (78a),  (78b)  and  (78c),  = - (N/I)  Z"^  and  in  (78d)  + 

1/4  Y Z(£/N)^ •!- (N/£)Z“ Therefore  when  writing  the  node-voltage 
equations  in  this  fashion,  it  appears  that  the  inverse  of  the  per-unit-length 
impedance  matrix,  Z,  is  required.  The  inverse  of  Z is  needed  when 
writing  node-voltage  equations  strictly  from  a circuit  diagram  when  mutvial 
inductances  are  present  [53],  However,  one  can  show  from  (78)  that 

Ikia  ik?a  iki?  " Ikaa 

A 

which  can  be  used  in  forming  Y and  Y., , , in  (80)  as 

~ '"Nrl 

In+1=  ik.»}  . <8'"'' 

Thus  the  inverse  of  Z is  not  needed  when  writing  the  node-voltage  equations 
in  this  fashion  as  in  (80). 

It  should  be  noted  that  the  formulation  in  (80)  provides  an  additional 
method  of  obtaining  the  exact,  distributed-parameter  solution  for  V(0)  and 
V(£),  For  example,  taking  N=1  and  using  the  distributed  matrix  chain 
parametere  $ig»  $33  from  (31)  and  (32)  one  obtains  a set  of 

2n  equations  in  the  2n  unknowns,  V(0)  and  V(£),  which  can  be  solved  for  iho 
exact,  distributed-parameter  solution  instead  of  using  (70),  (71),  (72)  or  (73). 
For  the  distributed  case,  one  can  also  show  $5ls  ^11  “^21  “ p (see  (40)). 
In  addition,  one  can  show  that  $13^33  ~ ^11  (see  (40))  and  this  can  be 
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used  in  furminj'  in  (HI).  Thus,  for  the  distributer!  case,  (80)  becomes 

{*!«'  5o‘  - 4<»}  in 

in  {i<«  sl‘  - 1;«} 

and  (85)  may  be  explicitly  solved  for  V(0)  and  V(£)  as  was  done  in  obtaining 
(73)  from  (72),  Note  the  similarity  of  (85)  to  (72), 

An  important  consideration  in  using  the  lumped-circuit  iterative  approxi- 
mations is  that  to  obtain  correlation  with  the  distributed-parameter  formula- 
tion described  by  (4)  (which  these  models  are  intended  to  approximate),  each 
section  must  be  electrically  short  and  therefore  the  number  of  sections  used 
to  represent  the  line  must  be  increased  for  increasing  frequency.  For  an 
(n+l)-conductor  line  with  N sections,  (N+l)n  simultaneous,  complex  equations 
in  terms  of  the  node  voltages,  V(0),  V(£/N),---,  V(r)  must  be  solved  at 
each  frequency  as  is  evident  from  (80).  For  N>1,  i.e.,  using  more  than  one 
section  to  represent  the  entire  iine,  the  equations  become  sparse  (large 
number  of  zero  entries)  which  is  clear  from  (80)  and  the  fact  that  each 
section  interacts  directly  with  only  its  two  neighboring  sections.  This  high 
degree  of  sparsity  can  be  used  to  drastically  reduce  the  storage  and  compu- 
tation times  over  that  which  would  be  required  if  the  nodal-admittance 
matrix  were  treated  as  "full"  and  no  advantage  taken  of  the  zero  entries. 
These  considerations  are  implemented  in  the  program,  TRAFFIC  l£3]. 

Even  if  only  one  section  were  used  to  represent  the  entire  line,  i.e.,  N=l, 
(80)  shows  that  2n  complex  equations  in  2n  unknowns  must  be  solved  at  each 
frequency.  However,  these  equations  can  be  solved  explicitly  for  V(0)  or 


V(0) 

y(D 

♦(£) 

~0  — 0 

«(X) 


(85) 
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V(£)  so  that  for  N=1  a minimum  of  n simultaneous  equations  in  n unknowns, 
V(0)  or  V(X)  need  be  solved.  In  obtaining  the  solution  of  the  distributed- 
parameter,  transmission  line  equations  directly  rather  than  approximating 
the  line  with  lumped-circuit  iterative  models,  one  is  also  required  to  solve 
2n  simultaneous  equations  in  2n  unknowns  through  a solution  of  (70),  (72)  or 
(85)  or  n equations  in  n unknowns  through  a solution  of  (71a)  or  (73a)  and  the 
number  of  simultaneous  equations  which  must  be  solved  need  not  be 
increased  with  increasing  frequency  as  is  required  with  the  lumped-circuit 
iterative  approximations. 

Incident  fields  may  also  be  incorporated  into  the  solution  via  the  lumped - 


circuit  iterative  approximations. 

Define  from  (8) 

^(kX/N) 

_I(kX/N) 

, (k-l)X 

= 4(kX/N,  /N) 

^((k-l)X^N) 

+ 

kX/N 

^ i(^^/N,x) 

(k-l)X/N 

Is(x) 

J 

(86) 


= 


r 

-,kX/N 

V (x) 
— s 

1 

(k-l)X/N 

dx  J 

(k- 1)£ 

for  k = 1,  2, , N since  it  is  assumed  that  £{kX/N,  ' /N)  :»k  for  elec- 

trically short  sections.  Thus,  for  electrically  short  sections,  the  lumped - 
circuit  iterative  models  can  logically  be  modified  to  include  incident  fields 
by  adding  appropriate  voltage  and  current  sources  to  the  beginning  of  each 
section  as  indicated  by  (86).  The  node-voltage  equations  in  (80)  will  be 
modified  by  adding  appropriate  additional  forcing  functions  to  the  right-hand 
side  vector  and  the  coefficient  matrix  will  remain  unchanged. 


The  question  of  convergence  of  the  lumped-circuit  iterative  approxi- 
mations to  the  distributed-parameter  solution  is  difficult  to  answer  quanti- 


tatively. A preliminary  indication  can  be  obtained  by  observing  the  con- 
vergence of  the  overall  chain  parameter  matrix  for  an  N section  represen- 
N 

tation,  i,  , to  the  distributed-parameter  chain  matrix  (or  state  transition 

~K 

matrix),  4.  The  state  transition  matrix  f(X)  (equations  (31),  (32)  and  (33) 
with  X = X and  x^  = 0)  can  be  expanded  into  an  absolutely  convergent 
infinite  series  as  shown  in  (34)  as 


i(X)  = 


“l  0 “ 
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Expanding  the  chain  parameter  matrices  for  the  lumped  iterative  approxi- 
mations in  (78)  one  obtains  for  the  lumped  T model: 
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and  for  the  lumped  T model: 
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and  for  the  lumped  Pi  model: 
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and  for  the  lumped  Tee  model: 
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Note  Shat  for  N - 1,  i.e.,  one  sottion  is  u«ed  U>  ropresont  the  cnliro  line, 
the  lumped  Pi  and  lumped  Tee  models  appear  to  be  better  approximations  to 
the  distributed  solution,  In  (87)  than  do  the  lumped  F and  lumped 

models  in  tl  sense  that  the  first  three  terms  of  (88c)  and  (88d)  are  Identi- 
cal with  tile  first  three  terms  of  $(D  and  the  fourth  term  is  only  partially 
the  fourth  term  in  $(£).  In  the  expansions  (88a)  and  (88b),  only  the  first 
two  terms  agree  with  the  first  two  terms  of  $(X)  and  the  third  term  partially 
agrees  with  the  third  term  in  $(£). 

There  are  c<»rtain  other  lumped  approximations  which  at  first  glance 
seem  to  be  not  included  in  this  discussion  but  are,  in  reality,  versions  of 
luinped-circuit  iterative  models  using  only  one  section  to  represent  the 
entire  line  and  with  certain  circuit  elements  neglected  In  addition, 

with  these  approximations  the  raparitivp  and  inductive  coupling  are  corru 
puted  independently  of  each  other  and  added  together  which  is  generally  only 
valid  for  weakly-coupled  lines  (see  reference  f30],  pp.  287-291). 
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V.  THE  PER-UNIT- LENGTH  PARAMETERS 


Dfrrivations  of  tho  per-unit-length  parameters  of  internal  resistance  and 
internal  self  inductance,  i.e.,  the  entries  in  R^,  and  L^,  which  account  for 
skin  effect  associated  with  imperfect  conductors  are  well  known  for  solid, 
round-wire  conductors  and  are  found  in  numerous  texts  ^Z,  3,  30j.  These 
internal  parameters  are  derived  by  assuming  that  the  currents  internal  to 
the  conductors  are  symmetric  with  respect  to  the  centers  of  the  conductors. 
However,  for  closely-spaced  conductors,  this  assumption  may  not  be  valid 
since  proximity  effect  can  alter  the  internal  current  distributions  (see 
reference  *^40],  Chapter  9). 

The  de  'vations  of  the  per-unit-length  external  parameters,  i.e.,  the 
entries  in  G,  L and  C,  assume  all  (n+1)  conductors  are  perfect  conductors 
and  are  more  involved  especially  for  close  conductor  spacings.  These 
parameters  generally  only  exist  in  closed  form  for  the  simple  cases  of  two- 
conductor  lines  in  homogeneous  media  in  Fig.  4 consisting  of  two  bare  wires 
in  an  infinite,  homogeneous  medium  in  Fig.  4a  (reference  [55],  pp.  133-136), 
one  bare  wire  in  an  infinite  homogeneous  medium  above  an  infinite  ground 
plane  in  Fig.  4b  (reference  [55],  pp.  183-185),  and  one  wire  within  a circu- 
lar shield  which  is  homogene ou j. ly  filled  with  a dielectric  in  Fig.  4c  (ref- 
erence [55],  pp.  125-133), 

Having  accurate  values  for  the  entries  in  the  per-unit-length  parameter 
matrices,  especially  the  external  inductance  and  capacitance  matrices  L and 
C,  is  obviously  important  in  obtaining  accurate  solutions  and  the  per-unit - 
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lenn*h  pa  r.t  motors  must  l>o  t^l)^aino<i  ovon  whon  usinj;  lmnpo<l-c  iriaiil  ilor.i- 
tivo  approximations  discussod  in  So c {.ion  4.1.  It  is  important  lij  romombor 
that  with  the  assumption  of  TEM  mode  propagation  on  the  line,  the  trans- 
verse fields  at  each  x along  the  line  satisfy  a static  distribution  and  there- 
fore the  per-unit-length  external  parameters,  i.e.,  e'ntries  in  G,  Land  C, 
are  obtained  as  the  solution  to  a two-dimensional  static  fields  proble'.ii  f39j. 
This  also  is  implied  in  the  inhomogeneous  medium  case  under  the  ''quasi- 
TEM  mode"  assumption. 

5.1  The  Pe r-Unit-Length  External  Parameters  for  Lines  in  a Homogeneous 

Medium 

The  per-unit-length  parameter  matrices,  C,  Land  G for  lines  irnmers- 
ed  in  a homogeneous  medium  possess  the  important  properties  given  in  (44), 
LC  = CL  = uci  and  LG=GL=uCJl„»  It  <^an  be  shown  I 54j  that,  for  a 
homogeneous  medium,  each  of  these  matrices  is  related  to  an  n\n  matrix, 
K,  which  is  independent  of  the  parameters  of  the  medium  and  dependent 
only  ''n  the  cross-sectional  structure  of  the  line  as 


C = € K 

(89a) 

G = aK  = 

(a/e)  C 

(89b) 

L = uK'^ 

= u e C 

• 

(B9c) 

For  two- conductor  lines  in  a homogeneous  medium,  the  per-unit-lenglh 

parameters  in  the  transmission  line  equations  in  (1)  are  obtainable,  exactly 

in  closed  form  even  for  close  conductor  spacings  where  proximity  effect 

produces  a nonuni'brm  charge  distribution  around  the  conductor  peripheries 

15  5],  The  rratrix  t'C  in  (89)  becomes  a scalar  K and  the  parameters  in  (1) 

•>/  * 
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become  c = e K,  g = a K,  i=  uK"^. 

For  two  wires  in  a homogeneous  medium  in  Fig.  4a,  the  per-unit-length 
capacitance  becomes  [56] 


2n  € 


c = 


-08h“'  [(d»-  rJ,j-  r^)/(2  r^  )] 

For  widely-separated  conductors,  (90)  can  be  approximated  by  [3] 


(90) 


c 


2tt  e 




(91) 


■where  in  is  the  natural  logarithm.  For  identical  wires  with  r , = r , , = r , 

wl  wO  w 

(91)  yields  less  than  5%  error  for  (d/r^  > 5 [55,  56].  For  the  case  of  one 
wire  in  a homogeneous  medium  above  an  infinite  ground  plane  as  in  Fig,  4b, 
the  per-unit-length  capacitance  becomes  [55] 


c = 2ne 

coeh“^ 

and  for  (2h/r^)  > 5,  (92)  can  be  approximated  by  [55] 


(92) 


2 n e 
Jin 


(93) 


For  the  case  of  one  wire  within  and  centered  on  the  axis  of  a circular  shield 
which  is  homogeneously  filled  with  a dielectric  as  in  Fig.  4c,  the  per-unit- 
length  capacitance  becomes  [55] 


c = 


2n  G 


?n  (r^/r^) 


(94) 


In  all  of  the  above  cases,  c = G K and  K is  easily  identified.  In  addition, 
^ = Uy  K and  g = a K.  For  Fig.  4a  and  Fig.  4b,  G = g and  a = 0. 
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The  p'irameters  for  lines  consisting  of  more  than  two  conductors  are, 
in  general,  not  obtainable  in  closed  form  for  closely-spaced  conductors  and 
numerical  approximations  must  be  used.  These  techniques  generally  fall 
into  two  classes;  the  meth«d-of-moments  techniques  [56-62,  22,  66  -68]  or 
imaging  techniques  [63,  64,  65],  A particularly  successful  technique  is  the 
use  of  harmonic  expamUon  functions  to  describe  the  charge  distribution 
around  the  conductor  peripheries  [56,  57], 

Consider  the  system  of  (n+1)  bare  wires  in  Fig.  2a.  With  the  moment 
method  using  harmonic  expansion  functions,  the  froe-charge  distribution, 
p. , around  the  cross-sectional  perimeter  of  the  i-th  conductor  per  unit  of 
length  in  the  x direction  is  described  as  a l''ouvioi‘  scries  with  respect  to  a 
cylindrical  coor<linate  system  at  the  center  of  the  i-th  conductor  as  shown  in 
F ig.  12a,  i.  c. , 


Ai 

Pi  ^^0  = ^iO  ^ ^im  Km 

i = 0,1, , n. 

The  absolute  potential^  0 (r  , 6 ),  at  an  arbitrary  point  P located  at  a 
radius  tp  and  angle  0p  shown  in  Fig.  12a  due  to  this  charge  distribution  over 
the  i-th  isolated  conductor  is  [56,  57) 

2tt 


*p'*V  «p’  = -7^  fio  '(Si*  'iSi 


Ai 


02n  p2n 

\ cos  (m0=)  I (Oj)  dO:  + r b;,^,  \ sin(m0i)  l(0|i  cie 

'^0  III  uri  1 


+ F a.  \ cos 

m-l  do 


1.  The  term  "absolute  potential"  reters  to  the  potential  with  respect  to  inl'in- 
Uy  and  the  reference  potential  terms  are  omitted  in  (96).  This  is  \-alid 
for  a system  with  xero  net  charge  and  is  demonstrated  in  Appendix  F, 
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The  parameters  for  lines  consisting  of  more  than  two  conductors  are, 
in  general,  not  obtainable  in  closed  form  for  closely-spaced  con  luc tors  and 
numerical  approximations  must  be  used.  These  techniques  generally  fall 
into  two  classes;  the  method-of-momencs  techniques  [56-62,  22,  66-68]  or 
imaging  techniques  [63,  64,  65],  A particulai  ;uccessful  technique  is  the 
use  of  harmonic  expansion  functions  to  describe  the  charge  distribution 
around  the  conductor  peripheries  [56,  57]. 

Consider  the  system  of  (n+1)  bare  wires  in  Fig.  2a.  With  the  moment 
method  using  harmonic  expansion  functions,  the  free-charge  distribution, 
p.v  around  the  cross-sectional  perimeter  of  the  i-th  conductor  per  unit  i.if 
length  in  the  x direction  is  described  as  a Fourier  series  with  respect  to  a 
cylindrical  coordinate  system  at  the  center  of  the  i-th  conductor  as  shown  in 
Fig.  12a,  i.  e.  , 


A.  Bi 

Pi  ^^0  = ^iO  + ^im  "^®i  + Jli  ^im  "^«i  ( > 

i = 0,1, , n . 

The  absolute  potential^  <t>  (r  , 9 ),  at  an  arbitrary  point  P located  at  a 

P P p 

radius  rp  and  angle  0p  shown  in  Fig.  12a  due  to  this  charge  distribution  o\'or 
the  i-th  isolated  conductor  is  [56,  57] 


■‘p'T'  V ' ■ 


1(6:)  d0, 


f96' 


^i  rZn  ^i  p2n  •, 

^ ^im  \ (tn0:)  I (0-)  d0.  + F , 6-  \ sin(m0i)  I(6T  d0.  i 

m=l  ^ Jq  III  '0  ‘ 


1.  The  term  "absolute  potential"  refers  to  the  potential  with  respect  to  infin- 
ity and  the  reference  potential  terms  are  omitted  in  (96).  This  is  valid 
for  a system  with  zero  net,  charge  and  is  demonstrated  in  Appendix  E. 


-85- 


potentiol  ^ 


whor<>  1(6;)  /'I  f 


- r^j  c-t).s  (6j  - 6pjjr^..  Those  inlourals 

can  be  evaluated  in  closed  form  yielding  1^6,  57]  (See  Apnendix  E ) 


®p>  = -‘iO  (■ 


/wi  (r  ) Ai 

( ^ 4-  V a 1 

cu.(me  ) 
f ’’  'l 

\ e / m=l 

\ 2mcr^^^  ' 

+ ^ 
m 


(m+1) 

r..^  sLn(m0p) 


b. 

=1  2mc  r”'  ^ 


Associated  with  each  conductor  in  Fig.  2a  there  are  + 1 un- 

knowns; the  expansion  coefficients  a_,  a.  , b.  in  (95)^  These  unknowns 
will  be  determined  by  enforcing  the  boundary  conditions  that  each  conductor 

(including  the  reference)  is  at  an  absolute  potential  i=0,  l,---n.  A total 
n 

of  Z (A.  + a.  +1)  match  points  will  be  chosen  on  the  (n+1)  conductors  at 
i=0  I 1 

which  the  potential  due  to  all  charge  distributions  in  the  system  (including 

the  charge  distribution  on  the  conductor  associated  with  the  match  point) 

n 

will  be  enforced.  This  results  in  a set  of  r (A.  + B.  +1)  equations  in 

i=0  ^ ^ 

the  same  number  of  unknowns  and  can  be  written  as  [56,  57] 

P p=  (i  (98) 

where  p is  a vector  of  length  .Fq(A^  + + 1)  containing  a i of  the  unknown 

expansion  coefficients  b..^  and  is  a vector  of  the  assumed  con- 


ductor potentials  as 


; j 


where  is  the  potential  of  the  j- th  conductor  at  each  of  the  chosen  match 
points  on  that  conductor  surface.  Inverting  P in  (98),  one  may  write 


P = P"^  0 
^ ^ • 


(100) 


From  the  solution  in  (100)  for  tlie  expansion  coer^icients,  the  total  free 
charge  on  the  i-th  conductor  is  [56,  57] 


^ Pi(8i)  r^i  d8. 
0 


and  an  (n-*-l)  y (n+1)  "generaliaed"  capacitance  matrix  may  be  written  a» 

’ol  r^oo$)r'C]pol 

• = vio  ^11.  : • (102) 

^nO  • * ' • ^nn  J ^n 

In  (102),  note  that  ^ C. . = q.  if  the  excitations  are  chosen  as  0.  = 1 

~ ij  - ij  I J 

with  all  other  potentials  chosen  zero,  i,  e, , 0p  = 0,p  ^j,  and  i,  j=0,  1, , in. 
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However,  from  (101),  qj^  = 2tt  ^ ^iO*  Therefore,  to  find  simply  add 
all  elements  of  P"^  in  (100)  which  are  in  the  row  associated  with  and  the 

~ lO 

columns  associated  with  4-  and  multiply  the  result  by  2tt  r . [56,  57],  From 

J w I 

energy  considerations,  one  can  show  that  ^ is  symmetric,  L.e., 

[58], 

The  n-y  n external  capacitance  matrix,  C,  used  in  the  transmission  line 


equations  in  (4),  (6)  and  (7)  where  potentials  Vj^,  i=l, ,n  are  defined  with 

respect  to  the  zero-th  conductor  chosen  as  reference  instead  of  absolute 
potentials,  Mj,  can  then  be  obtained  directly  from^  . To  !o  this,  note 

that  choosing  the  zero-th  conductor  as  reference,  = ((J^  - (Jq),  i = l, ,n 

and  the  transmission  line  capi'citance  naatrix  C becomes 


where  [C]..  | C--  for  i,j=l, , n. 

~ ij  ‘J 

Since  the  systerri  is  electrically  neutral,  we  have 


(103) 


‘>0  = - J,  ’l 

and  the  potentials  with  respect  to  the  reference  conductor  become 


(104) 


V;  = (®;  - »„)  . 


Thus  (102)  can  be  written  as 


(105) 


- 1 <?'ol  +<^02  ''2  ^ ^ 


Om  0 


^1  ^11  ^12  ^^2  "V  ^^In  ^n  Jn^O  ^lm'^0 


(106) 


= *nl  '"l  ^ nZ 


+ V„  + ( V C ) 
m=0  nm 


Adding  all  equations  in  (106)  together  we  obtain 


= ( Z C )v,  + ,)V2  + --  + f )v  (107) 

'^m=0  ml-'  1 'm=0^  ni2  ^ ^ m=0  mn' 


or 


*^0  " 


+ ( 

/ ^ Ja  n 

, s c + T 

t + — 

n 

+ r 

'm=0  Om  m=0 

Im 

m=0 

n 

V 

[^m=0^mk^  ^k] 

k=l 

ilo  [m^O  ^?,JnJ 


(108) 


Substituting  (108)  into  the  last  n equations  in  (106)  and  arranging  these  last  n 
equations  in  the  form  of  (103)  yields  the  entries  in  the  per- unit-length  trans- 
mission line  capacitance  matrix  C as 


[C] 


( S w ^ ^ 

-1  - C “dr  ^ni=0  Vra=0  mj/ 

1=0  \m=0  ^ IT\J 


(109) 


r ” ( 

t T 

^ - { "n 

t.  ) 

( ■) 

ij  U=0 

^m=0 

im/. 

J \m=0 

im/ 

\m=0  mj/ 

Vf.=0  m=0  ^ imJ 


Note  that  the  denominator  of  (’/^9)  is  simply  the  sum  of  all  the  elements  of 
d.  The  numerator  of  (109)  can  be  written  as 


dr"  T t i-drr"dvd,.(”d. 

ij  L/=0  m=0  ^mj  ^ij  \m=0  mj/  ij  \m=0  im/ 

mj«fi  m^j 


m^j  m^i 


U 


(110) 
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r " T & - " " dr  ^ 

ij  L /,=0  rn-0  ^ /m  m~0  »nj  m=0  im  ^ij-1 

m^i  m^j 


sum  of  all  terms  in^ex'-ept  those  in  the  L-th 
row  and  j-th  volumn. 


\m=0  im/  \m=0  mj/  • 


rn^i 


sum  of  all  sum  of  all 

terms  in  terms  in  the 

the  i-th  j-th  column  of  it 

row  of  except 

except  . 


Therefore,  can  be  written  as 


Clj  = 


(?ij  (mq  - jM-  Mj-g,j)  ■ jM  M) 

where 

Mg  = Sum  of  all  terms  in 

,M  = Sum  of  all  term.s  in  the  i-th  row  cf^  except 

1 - ij 

Mj  = Sum  of  all  terms  in  the  j-th  column  of  ^except 
For  two  conductors  (n=l)  (ill)  becomes 


(111) 


qi/Vi 


Cl  d c 

00  ~ 01  10 

^00  ^ ^01  ^ ^10  ^ ^11 


c - 


(112) 


This  rosull  for  n~l  is  also  oblaine<l  in  roferonce  (58),  pp.  211-213. 

This  numerical  technique  can  also  be  applied  to  the  systems  in  Fig.  2b 

by  replacing  the  ground  plane  with  the  conductor  images  having  equal  but 

opposite  charge  distributions.  As  an  example,  consider  the  system  of  n 

bare  wires  above  an  infinite  ground  plane  in  Fig.  2b.  Replacing  the  ground 

plane  with  images  results  in  Fig.  13.  Note  that  the  orientation  for  6^  on  the 

image  of  the  i- th  conductor  is  the  same  as  the  orientation  for  0.  on  the  i-th 

conductor  however,  the  cliar^ie  distribution  on  the  i-th  conductor  image  is 

-p^  (-0^).  Also,  note  that  the  potential  of  the  i-th  image  conductor  is  taken 

to  be  By  symmetry  and  the  use  of  image  distributions,  the  voltage  of 

the  i-th  conductor  with  respect  to  the  ground  plane,  V^,  is  equal  to  since 

the  potential  diffexence  between  the  i-th  conductor  and  its  image  is  (/i-  - 

2(i..  Therefore,  we  only  need  to  enforce  the  potential  at  match  points  on 

the  n conductors  above  t.ie  ground  plane  due  to  all  charge  distributions  m the 

system  (those  on  the  n wires  and  on  the  n image  wires).  Thus  a set  of 
n 


E (A.  + B.  +1)  equations  can  be  written  as  in  (98).  However,  these  equa^ 
i=l  ^ 

tions  will  differ  from  those  in  (98)  in  that  the  expansion  coefficients  for  the 

zero-th  conductor  (the  ground  plane),  ^Om’  ^Om’  included 

in  the  vector  £ and  the  potentials  will  not  be  included  in  the  vector  <t> . 

n n 

These  vectors  will  be  of  length  ^ (A.  + B.  +1)  instead  of  T (A.  + B.  +1). 

i=l  * ‘ i=0  * 1 

Furthermore,  we  may  replace  in  d by  Inverting  P we  then  may  obtain 

the  entries  in  the  transmission  line  capacitance  matrix  directly  without  the 

need  for  (ID)  since  d.  = V,.  Therefore,  C. , is  simply  the  sum  of  all  elements 

II  ij 
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in  P"^  which  are  in  the  row  associated  with  a^Q  and  the  columns  associated 
with  cf)j  = Vj  multiplied  by  2rrr^.  for  i=l, ,n  and  j=l, , n. 

Consider  the  case  of  n wires  within  a circular  shield  shown  in  Fig^  2c. 
Rather  than  imaging  the  conductor  charge  distributions  across  the  shield 
boundary,  one  may  expand  the  per-unit- length  free-charge  distributions 
around  the  n conductor  peripheries  in  a Fourier  series  and  also  expand 
the  per-unit-length  free-charge  distribution  around  the  interior  periphery  of 
the  shield.  Note  that  in  this  case,  the  voltages  of  the  n wires  with  respect 
to  the  shield,  V^,  will  be  Thus  (111)  can  be  used  to  obtain  C. 

Note  that  for  all  these  cases,  once  Qis  determined,  and  Q,  are 
obtained  through  (89)  as  L S " C.  For  Fig.  2a  and 

Fig.  2b,  c = €v  a^nd  a = 0. 

It  is  also  possible  to  obtain  closed-form  approximations  for  Q,  L and 
under  the  assumption  that  the  conductors  are  widely  spaced  and  these 
formulas  represent  the  predominant  method  of  computing  the  entries  in 
these  per-unit-length  matrices  [55,  57j.  If  it  is  assumed  that  the  wires  in 
Fig.  2a  are  sufficiently  separated  so  that  the  charge  distributions  around 
the  wire  peripheries  are  constant,  i.  e. , proximity  effect  [55,  56,  57,  58] 
is  not  a factor,  then  only  one  expansion  function  is  needed  in  (95),  a^Q.  In 
this  case  there  are  only  (n+1)  unknowns  , aj^Q,  i=0,  1, , n.  Further- 

more, since  the  wires  are  assumed  to  be  widely  spaced, then  fp  in  (97) 
can  be  taken  to  be  d^Q,  djQ  or  dy,  whichever  is  appropriate,  when  compu- 
ting the  contribution  to  the  potential  of  a conductor  due  to  the  charge  on 
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another  conductor.  When  computing  the  contribution  to  the  potential  of  a 
'onductor  due  to  the  charge  distribution  on  its  boundary,  the  match  point 
is  taken  on  the  conductor  surface,  i.e,,  r =r^j^.  This  assumption  of 

tr 

widely-spaced  conductors  has  been  consistently  used  in  the  power  trans- 
mission area  [3,  4,  13,  55,  57] and  is  generally  valid  if  the  smallest 
ratio  of  wire  separation  to  wire  radius  is  on  the  order  of  5 or  greater 
[ 55,  56,  57]. 

Consider  Fig.  2a.  Assume  as  an  approximation  that  the  (n+1)  wires 
in  Fig.  2a  are  sufficiently  separated  such  that  the  per-unit-length  charge 
distributions  around  the  wire  peripheries  are  constant  with  each  wire 
bearing  a per-iaiit-length  total  free  charge 

^i  ~ ‘^''^^wi  ^iQ  (113) 

and 

Pi  (6i)  = ^iO  (114) 

for  i = 0,  1, , n.  Because  of  the  assumed  large  separations,  this  is 

equivalent  to  replacing  the  wires  with  filamentary  line  charges  [55]. 

Since  the  wires  are  assumed  to  be  widely  separated  so  that  the  assumption 
of  a constant  charge  distribution  around  each  periphery  is  valid,  i.e., 
proximity  effect  [55]  is  not  a factor,  then  we  may  choose  match  points  for 
the  potential  at  the  centers  of  the  conductors  rather  than  at  some  point  on 
the  periphery.  However,  for  a match  point  on  the  conductor  surface 
bearing  the  charge  under  consideration,  we  take  the  match  point  on  the 
conductor  surface,  i.e.,  rp  = in  (97). 
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From  the  previous  results  and  (113),  we  may  then  write  (utilizing  the 
first  term  in  (97)) 


^0 

f’ii(rwo)  -^n(dio)*  • •»£n(dnQ) 

- ” 

^0 

^1 

• 

= - 1 

;.n(diQ)  in(rwj  ) in(d^2)  * 

* • 

• 

^1 

g 

• 

• 

<6n 

2rre 

* • 

« • 

f • • 

0 

0 

^n 

L _ 

Applying  (105)  to  (115)  results  in  a typical  equation 

= c6i  - (f>Q 


+ + 


4 In  AwiN  qi 


for  i,  j - 1, , n.  Applying  (104)  to  (116)  yields 

4- 


Vi  =- 


(115) 


/---  4-  jgn  Awi^'wp  \ q; 

2TTe  ^ V^  io'^iO  ) 

q^  4-- -4-  Zn  /^in^wn  1 q..  \ 

\ J ’ ” J 


--  4-  £n 


Comparing  (117)  to  (10  3)  shows  that 

]ii  ""  77f?  (p^r  ) 

i^’j  2ire  \rwodij/ 

for  i,  j=l, , n. 


(116) 


(117) 


(118a) 

(118b) 
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(Il9ai 


Thus  the  per-unit-leng th  inductance  matrix  is  from  (89) 
[iilii  ~ 


■ 11 a / \ 

2rr  UwFwO y 


[fe]lj  =u€fe-']n 

iij  ijj  2" 

for  i,  j:=l, , n.  Note  that  for  Fig.  2a,  e = ey,  u = Uy  o = o* 


(119h) 


Similarly,  large-separaftion  approximations  can  be  obtained  for  the 
case  of  R wires  above  an  infinite  ground  plane  in  Fig.  2b.  Consider 
Fig,  13  and  assume  that  the  wires  are  separated  sufficiently  from  each 
other  and  the  ground  plane  so  that  the  per- unit-length  charge  distributions 
around  the  peripte  .''.s  of  the  wires  and  their  images  are  constant  and  given 
by  (113)  and  (114)  [ji  i-i, , n.  Again,  due  to  the  large-separation  approx- 

imations^ we  n\ay  take  the  match  points  at  the  centers  of  the  conductors 
when  computing  the  contribution  to  the  potential  due  to  the  charge  distri- 
bution on  another  conductor  and  take  the  match  point  on  the  conductor 
periphery  when  computing  the  contribution  to  the  potential  of  this  conduc- 
tor due  to  its  own  charge  distribution.  From  Fig.  13,  a typical  equation 
for  the  i- th  conductor  may  be  written  as  (again  using  only  the  first  term 
in  (97)) 


+ ^iO 


wi-^n(r-yyi)  ^ ryy.j;,n(2h 


(120) 


f-  (<^ij)  ^ ^wj^^(djj)  \ ^ 

JO  1 € e J "" 

=z4-{---  + ■’j  ^ -} 

for  i,  j=l, , n where  d^'j  is  the  center- to-center  separa  tion  between  the 


i-thwire  (j- th  wire)  and  the  image  of  the  j- th  wire  (i-th  wire)  given  by 
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di'j  = >J  (hi  + hj)®  + _ (hj-hi)3 


(121) 


= J d.?  + 4 hihj 

Comparing  (120)  co  (103)  and  using  (89),  we  obtain 

=4ffW4fi-) 


for  i,  j=l, , n.  Note  that  for  Fig.  2b,  e=e^,  u=  Uv  ^^nd  a=0. 

Large-separation  approximations  may  also  be  obtained  for  the  system 
of  n wires  within  a circular  shield  in  Fig.  2c.  We  assume  that  all  n wires 
are  sufficiently  separated  from  each  other  and  the  shield  so  that  the  wires 
may  be  replaced  by  filamentary  line  charges.  The  circular  shield 
(assumed  to  be  perfectly  conducting),  may  then  be  replaced  by  filamentary 
line  charge  images.  Each  filamentary  line  charge  has  an  image  on  a 
line  joining  the  line  charge  and  the  center  of  the  shield  and  is  at  a dis- 
tance of  Ts/ri  from  the  shield  center  where  rg  is  the  shield  radius  and 
ri  is  the  radial  distance  of  the  i-th  wire  from  the  shield  center  [58].  One 
can  then  straightforwardly  derive,  by  superposition 

^]ii=—  I (I23a) 

9i  I , , qn-  0 


Vi 


'j+l 


— . q„=o 


(123b) 


where  6ij  is  given  in  Fig.  2c. 


i rg*  - 2r[r 


^rs^  cose^i 
Co's  0i 
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To  derive  (123),  consider  an  infinitesimal  line  charge  of  radius  r 


w 


bearing  a total  per~unit-length  free  charge  q.  The  potential  at  a point 
r respect  to  the  line  charge  surface  as  reference  is  (see  reference 

[58],  page  92) 


V = - 


jfn 


2ne, 


(124) 


Equation  (I23a)  can  then  be  derived  from  Fig.  14a  with  (124)  and  q - 0 as 

■•■A  I - (^t^-  • ■■  l%^)j 


q. 

1 

2tt  < 


Equation  (123b)  can  be  derived  from  Fig,  14b  by  making  use  of  the  result  in 
(124)  and  the  law  of  cosines.  With  q^^  = 0,  we  obtain 

fdi 


Vi  = 


2tt£, 


In 


rwj  / Vwj  / \^*‘wj  ) V^wj 


2TTe, 


£ n 


(126) 


Utilizing  the  law  of  cosines,  one  may  obtain 

+ (r^  /r,)^  - (Zr^/rj)  cos  9ij 

4 = - (2r.r^^/rj)  cos  Q-- 


e 


U 


d,  ^ = r ^ + r?  - ^2r  r;)  cos 
1 s j ' s j/ 


(I27a) 

(127b) 

(127c) 

(127d) 


Substituting  (127)  into  (126)  yields  (I23b). 
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Figure  14.  The  replacement  of  the  wires  ard  the  shield  with 
infinitesimal  line  charges  for  shielded,  multicon- 
ductor lines  with  large  separations  in  homogeneous 
media. 


5,  2 The  P»»r-Uni t-Length  External  Parameters  for  Lines  in  an  Inhomo- 
geneous Medium 

Derivations  of  the  entries  in  G,  !.>  and  Q for  the  lines  in  Fig.  3 are 
complicated  by  the  inhomogeneity  of  the  surrounding  medium  introduced 
by  the  interface  between  the  dielectric  insulations  and  free  space.  The 
inhomogeneity  is  introduced  through  the  permittivities  and  ccinductivities 
of  the  insulations  since  the  dielectric  insulations  are  characterized  by 
free  space  permeability,  Uv  Therefore,  the  per-unit-length  external 
inductance  matrix,  1^,  can  be  found  as  where  Cq’-  is  determined 

as  in  Section  5.1  with  the  dielectric  insulations  removed  and  may  be  ob- 
tained accurately  with  moment  methods  and  harmonic  expansion  func- 
tions or  may  be  approximated  for  large  conductor  spacings  by  (119),  (122), 
or  (123), 

The  computation  of  the  entries  in  G introduces  some  conceptual  dif- 
ficulties for  the  inhomogeneous  medium  case.  Consider  the  computation  of 
Q for  static  excitation  (reference  [39],  Chapter  6).  The  transverse  con- 
duction  current  density  p j,  is  related  to  the  transverse  electric  field, 

C'I'i  each  dielectric  as  The  boundary  conditions  on  the 

potential  function  where  P-p  = -grad  is  tha t j,^^mu s t vanish  over  the 
conductor  surfaces  (perfect  conductors  are  assumed  for  this  computation) 
and  the  derivative  of  ^ normal  to  the  dielectric-free  space  boundaries 
must  vanish  at  the  boundary  (see  reference  [39^  chapter  6),  The  latter 

•4 

requirement  insures  that  the  normal  component  of  is  zero  at  <.Iiese 
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boundaries,  i.  e.  no  current  flow  into  the  free  space  medium.  The  con- 
sequence of  this  is  that  if  none  of  the  dielectric  insulations  touch  each 
other  or  the  ground  plane  or  circular  shield,  then  £|computed  for  static 
conditions  with  a straightforward  application  of  the  above  boundary  con- 
ditions would  be  identically  zero  since  no  conductive  path  between  the 
conductors  would  exist  for  nontouching  dielectrics.  However,  dipole 
relaxation  effects  will  nevertheless  produce  certain  losses  even  for  a 
transverse  field  distribution  and  nonstatic  excitation  since  the  trans- 
verse displacement  current  will  have  a portion  in  phase  with  the  transverse 
electric  field.  Therefore,  equivalent  shunt  conductances  should  be  de- 
termined to  represent  these  non- static  losses. 

Assuming  perfect  dielectrics,  however,  one  can  compute  the  en- 
tries in  C in  a straightforward  fashion.  A moment  method  of  solution 
with  harmonic  expansion  functions  as  in  Section  5..'  can  be  used  for  this 

problem  [b6,  57],  Consider  the  system  of  (n+1)  dielectric-insulated 
wires  in  Fig.  3a.  Represent  the  bound-charge  distributions  at  the  di- 
electric-free space  boundaries  with  Fourier  series  as  [56,  57] 

A A 

A . * 

A (6.)  = ^;n  + 5:  a im  cos  mO^  + T b.  sin  md-.  {128  ) 

1 I n-i=i  rn-1  i 

i-u,  i.,  ---,  n. 

Represent  the  charge  distributions  at  the  conductor- dielectric  boundaries 
(which  is  total  charge,  bound  plus  free  for  this  case)  with  Fourier  series 
as  in  (95).  The  contributions  ^o  the  potential  and  electric  field  at  a point 
P in  Fig.  12b  due  to  each  of  the  components  of  the  charge  distributions 

are  given  in  Table  I with  respect  to  Fig.  15  [56,57]. 
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TABLE  I 

Contributions  to  the  Potential  and  Electric  Field  at  a Point  P in  Fig.  15 
due  to  Harmonic  Expansion  Functions  on  a Circular  Boundary. 


Expansion 

Contribution  to  the 

Contribution  to  the 

Potential  at  P 

Electric  Field  at  P 

1 

. rb  /n(rn) 

cos  m6^ 

(m+1) 

r^^  cos  m©p 

2ev  ^ 

,na+l 

/cos  me  r 

+sin  me^0| 

sin  m 

(m+1) 

sin  me^ 

m+1 

(rb/^^  m0  v» 

•y  ^ ill 

^ 

2ev  1 

-cos  m0pe> 

(a)  rp>r^ 


Expansion 

Function 


Contribution  to  the 
Potential  P 


Contribution  to  the 
Electric  Field  at  P 


A A 

Thus  there  are  2 + A.  + B.  + Aj  + B.  unknowns  associated  with  each 

I I ‘ i 

A A A 

wire,  a^Q,  a^m»  ^im'  ®^i0»  *im»  The  boundary  conditions  will  be 

enforced  by  requiring  that  the  potential  on  the  i- th  conductor  due  to  all 
source  distributions  be  and  the  normal  component  of  the  displacement 
vector  due  all  source  distributions  be  continuous  at  the  dielectric-frec 
space  boundaries.  Generally  14-A^  + rnatch  points  are  selected  on  the 

A A 

i-th  conductor  and  1 + match  points  are  selected  on  the  interface 

between  the  i-th  dielectric  and  free  space.  The  component  of  the  total 

electric  field  (from  all  source  distributions)  normal  to  and  just  I- side 

the  '^ielectric-free  space  surface  at  each  match  point  on  this  surface  is 

multiplied  by  a^nd  set  equal  to  the  product  of  and  the  component  of 

the  total  electric  field  (due  to  all  source  distributions)  nornral  to  and  just 

n A A 

outside  the  boundary  at  this  match  point.  A set  of  (2  + Aj  + B^  + Ai  + Bj) 
simultaneous  equations  can  be  written  to  enforce  these  conditions  as 


[56,  57] 


[•  1 [i  1 ■ w 


(129) 


where  and^are  defined  in  (99)  and^is  a column  vector  of  the  expansion 

coefficients  in  (12  8)  arranged  as  in  p and  0 is  a column  vector  of  zeros  of 
n A A 


length  V (1  + Aj  + E;). 

i“0  ' 

Inverting  P in  ^129),  one  can  obtain  [56,  57] 


r £ 1 

r ' 

" 0 ' 



= I 

i ~ 

0 

L — « 

L 

(130) 
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The  total  free  charge  on  the  i-th  conductor  (which  defines  the  general- 
ized capacitance  matrix)  is  given  by  [56,57] 

qf!  = qi  + 4;  = \ PI  (®l)  rwi  ^ ^ ^i(ei)  (r^i  + ti)  d6i  (131) 

0 0 

= a^O  + 2n  (r^^  + tj  )a^Q 

since  q^^  is  the  total  charge  at  the  conductor,  dielec  trie  boundary  which  is 
the  sum  of  the  total  free  charge  and  bound  charge  with  the  bound  charge 
being  identical  in  magnitude  but  opposite  in  sign  to  the  bound  charge  on  the 
dielectric-free  space  surface,  4i*  e. , qj^  = q£^  - and  t^  is  the  thick- 
ness of  the  i-th  dielectric.  The  generalized  capacitance  niatrix  can  be 
written  a?  in  (102)  vdiere  q^  in  (102)  is  replaced  by  from  (131).  By 
using  the  excitation  (4j  = 1,  = 0,p  [$]ij  equals  the  sum  of  two  terms. 

One  term  is  2tt  r^^  multiplied  by  the  sum  of  all  elements  in  P"^  which  are 
in  the  row  associated  with  aj^Q  and  columns  associated  with  (Aj  and  the 
other  term  is  2tt  (r^j^  + tj^)  multiplied  by  the  sum  of  all  elements  in 
which  are  in  the  row  associated  with  and  the  columns  associated  with 
02  [56,  57 ]i  The  entries  in  the  n x*'  per-unit-length  capacitance  matrix, 
used  in  the  transmission  line  equations  can  then  be  straightforwardly  ob- 
tained from  the  (n+1)  y (n+1)  generalized  capacitance  mairUv,  ^ , as  in 
(111). 

Typical  computed  results  for  two  dielectric-insulated  wires  ax  : 
shown  in  Fig.  16.  By  symmetry,  the  coefficients  of  the  terms  in  (95) 

(128)  are  zero,  i.  e.  = 0 and  = 0.  Therefore  the  expansion  func- 
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Figure  16.  Computed  results  for  two  diele  ctric- insulated  wires 


tions  on  each  boundary  consist  of  the  constant  term  and  only  cosine  terms. 
Similar  results  for  a 5-wire  flat  pack  cable  are  shown  in  Fig,  17, 

Selected  entries  in  the  first  tow  of  the  generalized  capacitance  matrix, 

^00"  ^01*  ^02 » ^03'*  ^04  shown.  Again  by  symmetry,  the  sine  ex- 
pansion fxinctions  are  not  included  since  the  coefficients  of  these  terms  will 
be  zero  for  wires  in  a linear  array  such  as  flat  pack  cables. 

These  results  can  be  extended  to  include  the  case  of  n wires  above 
an  infinite  ground  plane  in  Fig.  3b  in  the  following  manner.  Consider 
the  set  of  n dielectric-insulated  wires  above  an  infinite  ground  plane  shown 
in  Fig.  18.  To  treat  these  cases,  we  replace  the  ground  plane  with  a cor- 
responding set  of  image  wires.  The  i-th  wire  image  is  at  a distance  of 
h^  below  the  ground  plane.  Each  image  wire  is  identical  to  its  correspond- 
ing wire  above  the  ground  plane  and  the  potential  of  the  i-th  conductor 
image  is  -(^^.  The  charge  distributions  around  the  i-th  conductor  and 
i-th  dielectric- free  space  boundary  are  denoted  by  0^(0^)  and  q(0p^respec- 
lively.  The  charge  distributions  on  the  corresponding  boundaries  of  the 
image  wires  are  identical  in  magnitude  but  opposite  in  sign  to  those  of  the 
corresponding  wires  above  the  ground  plane,  i.  e.  -p^(-0p  and 
The  charge  distributions  pj^  and  2^  are  again  expanded  into  Fourier  series. 

A set  of  simultaneous  equations  in  terms  of  the  unknown  expansion 
coefficients  in  p-  and  can  again  be  formulated  to  enforce  the  boundary 
conditions  on  the  potential  of  the  i-th  conductor  and  the  continuity  of  the 
normal  components  of  the  displacement  vector  at  the  dielectric-free  space 
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Figure  17.  Comouted  results  for  a five- wire  ribbon  cable 


i' 


I 


boundaries  of  each  wire  due  to  all  charge  distributions  in  the  system 

(the  distributions  on  the  wires  and  their  imagesh  By  symmetry  and  the 

n A 

use  of  image  distributions, we  only  need  to  write  T (2  + A + B. ) 

i=  1 ^ 

equations  to  enforce  these  boundary  conditions  on  only  the  original  n wires 
above  the  ground  plane.  Once  these  equations  are  solved,  the  per-unit- 
length  transmission  line  capacitance  matrix,  C , can  be  directly  obtained 
as  before  since,  for  this  case,  where  are  the  transmission  line 

voltages  with  respect  to  the  ground  plane  as  shown  in  Fig.  18.  Thus  for 
this  case  as  for  Fig.  2b,  there  is  no  need  to  reduce  the  generalized  capa- 
citance matrix  to  the  transmission  line  capacitance  matrix  via  (111). 

The  per-unit-length  transmission  J'  ; inductance  matrix,  L,  can  be 
obtained  accurately  by  repeating  this  solution  with  the  insulation  dielectrics 
removed  and  using  (89)  as  indicated  in  Section  5.1  or  using  the  large-separa- 
tion approximation  in  (122). 

The  solution  for  n wires  in  a circular  shield  in  Fig,  3c  can  be  ob- 
tained in  the  same  fashion  as  discussed  in  Section  5.  1 for  the  case  of 
Fig,  2c. 

The  above  discussion  of  the  solution  for  C assuming  perfect  dielec- 
trics indicates  a method  for  incorporating  dielectric  loss  and  therefore 
obtaining  an  equivalent  per-unit-length  conductance  matrix,  G,  to  repre- 
sent these  losses.  If  each  dielectric  permittivity  is  considered  to  be 
complex,  i,  e.  , ~ " - j(adi/iy)  then  (129)  can  be  formulated  as 

above  with  the  only  difference  being  that  P will  now  be  complex.  In 


mi,  '>  .ii 
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A A 

particular,  the  last  1 + Aj^  + rows  of  P will  be  complex.  Thus  the  n\  n 
complex, capacitance  matrix  can  be  obtained  as  above  a.s  Q,  = + jCj  so 

that  j-jQ  = - 0)^  the  imaginary  part  of  C can  be  identified  as 

G = - (li®  Cj , The  real  part,  Cj^,  is  identified  as  the  usual  capacitance  ma- 
trix and  will  of  course  not  be  the  same  as  the  matrix  which  would  be  com- 
puted assuming  a perfect  dielectric. 

Large-separation  approximations  can  also  be  obtained  in  a fashion 
similar  to  Section  5.  1 by  requiring  that  the  separation  between  all  wires 
be  large  enough  so  that  the  charge  distributions  around  the  dielectric-free 
space  boundaries  and  the  conductor- '^'electric  boundaries  are  essentially 
constant.  However,  this  is  generally  not  the  problem  of  interest  since 
wires  are  closely-coupled  in  densely-packed  cable  bundles  and  flat  pack 
and  ribbon  cables  and  it  is  to  be  expected  that  the  charge  distributions 
around  the  boundaries  will  exhibit  large  variations. 

However,  to  illustrate  the  application  of  the  technique,  the  per-unit- 
length  capacitance  between  two  dielectric-insulated  wires  which  are  widely 
separated  will  be  derived^  Consider  the  case  of  tv/o  dielectric-insulated 
wires  shown  in  Fig.  19.  Because  of  the  assumption  of  large  separations, 
we  may  assume  that  the  charge  distributions  around  the  conductor- dielec- 
tric boundaries  and  the  dielectric-free  space  boundaries  are  constant. 
Therefore,  only  four  expansion  terms  are  needed  in  (95)  and  (128),  ^qO' 
^00’  ^’0’  ^10'  facilitate  the  derivation  and  provide  an  upper  bound  on 
the  pe r-unit-length  capacitance,  we  will  take  the  match  points  at  the 
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boundaries  on  a line  joining  the  centers  of  the  wires  as  shown  in  Fig.  19. 


From  the  previous  results  and  Table  I,  one  may  obtain 

fwO  .tn(r^)  ^ -ain  f 


“ ^00 


Xn(r^)  ^ -ajQ  ^ ^wl  jtn(d-r^)  ^ 


(132a) 


f (J^w0+ k)  ) An  (r^  to)\  * 

«»  V “ 


^ (*^wl  4-  tl)  (d-r^)  \ 


- - “00  ) -“10  ( ) (132b) 

~^00^  + tp)  An  (d-rwl)^  ^ *~1^  An  (r^  + t|^)  ^ 

°°  \ (rwO  + tO)  J \ i<^-rwO-k)  ) / 

^“00  (-‘)  -^0  ( ) 

V (d-TwO  - Iq)  / 

*00  ^ 


(132c) 


^wO  ^ *^rl  ~^) 
(d-Twl  -fcl) 


] .aio  ( ] 

I \ (i-wl  + tl)  > 


(132d) 


= 0 


where  ®r0  = ^0  / ®v  and  fifi  = ®1  / ®v 

These  equations  may  be  solved  for  numerical  values  of  the  parameters 


as  outlined  previously.  However,  in  literal  notation,  the  solution  is  quite 
complicated.  Therefore,  to  simplify  the  solution  and  obtain  a closed  lorm 
expression  for  the  per-unit-length  capacitance,  we  will  assume  that  the 


wires  are  identical,  i.  e. , 


(133a) 


tQ  = ti  = t (133b) 

€q  = 6i  = e (133c) 

= ®/ev  • (133d) 

With  this  assumption  we  may  take  the  charge  distributions  on  wire  0 to  be 
identical  to  the  charge  distributions  on  wire  1 by  symmetry.  Furthermore, 
= 2rrr^  a^^  (134a) 

= q 


qo  = 2i-rr^oO 
= -q 

'll  = 2tt  (r^  + t)  a^Q 

= q 

= 2n  (j  ^ ^ t)  ^00 


(134b) 
(134  c) 

(134d) 


= -q 


Substituting  (133)  and  (134)  into  il32)  yields 


in/ 

\ 

q + Jtoi 

\ 

. a-  x*w  1 

H -r  ^ 

Ld-r  ! 
* w 

m 

/d-r  w \ 

i q + Jtn  j 

/d-r 

r ^ 

\ i 

f 1 

(r^t 

'■) 

)'i 


/(€r-l)  + 

(Gr-i)  \ 

q + 

f (®r-l) 

1 ^ 

\ (^w+t) 

(d-r^t)  J 

(r^  +t)  i 

f (S-l)  (‘^r-D  \ 

Y(d-rw-t)  (’’w 

Note  that  (135a)  and  (135b)  show  that 


1) 


t) 


_L_)  4 = 0 

+«  I 


(135a) 

(135  b) 

(I35c) 

(I35d) 


(t)Q  - - 

and  (135  c)  and  (135 d) are  identical.  Therefore 

V = (J|  - (6q 
= 2^1  . 
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(136) 


(137) 


Note  also  from  (131)  that  = q + q and  therefore 

q = qf  - q • 

Substituting  (138)  into  (135b)  and  (135c)  with  (137)  yields 
V=  — f ) q£  + \ * 1 

I \ / Uw+t  / J 


(138) 


(139^ 


a = (er-l)d 


q = Of 

er(d-r^-t)  ^ 


(140) 


Substituting  (HO)  into  (139)  yields  the  per-unit-length  capacitance  as 

c=^ 

V 


Sf2LzJl±.\tn(  • (HD 

^ V i \er(d-r^t);  [ r^  / 

As  a check  on  this  result,  note  that  for  €j.=l  and  t=0,  (141)  reduces  to 
(91)  for  identical  wires  and  large  separations  such  that  d-r^  s d* 

As  a final  illustration  of  these  methods,  we  will  compute  the  entries  in 
the  per-unit-length  transmission  line  capacitance  matrix  for  the  case  of  two 
dielectric-insulated  wires  above  an  infinite  ground  plane.  In  order  to  sim- 
plify the  procedure  and  to  obtain  closed-form  expressions  for  these  quanti- 
ties, we  will  assuxTie  that  (1)  both  wires  are  identical  and  are  at  the  same 
height  above  the  ground  plane,  and  (2)  the  wires  are  separated  sufficiently 
from  each  other  and  the  ground  plane  so  that  we  may  assume  constant  charge 
distributions  around  the  conductor  and  dielectric  boundaries,  i.e. , 
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(142a) 


i 


q,  = Znr  a 

* I 


w“lO 


^1  = 

^2  = 2rrr^ja2o 


("d  = 


(r  , = r + t) 
' d w ’ 


(142b) 

(142c) 

(142d) 


The  ground  plane  will  be  replaced  by  images  as  shown  in  Fig.  20.  We* 

f A 

choose  the  match  pointsl,  I,  2,  and  2 on  a line  joining  the  two  wires 
(other  choices  are  of  course  possible). 

Applying  the  results  of  this  section  and  utilizing  symmetry  we  only 

A 

need  to  write  constraint  equations  at  match  points  1 and  1.  The  equation  for 
the  potential  at  1 becomes  (see  Table  I) 

*l  = ’1  + “ <‘'d> 

+ fn(d-r^.)q2  + fn(d-r^)q2 

+ jen(dj)(.qj)  + in(djH-qj) 

+ fn(d2)(-q2)  + fn(d2)(-q2)| 

A A 

The  distances  d^,  dj,  d2  and  d2  in  Fig.  20  are  given  by 


(143) 


di  = /4h2  + rj- 

J 


A 

d. 


4h‘‘  + r. 


=/4h^  -I-  (d-r^„)^ 

d^  =/4h^  + (d-rd>^ 


(144a) 

(I44b) 

(144c) 

(144d) 


The  equation  for  the  continuity  of  the  nornial  component  of  the  displace- 
ment vector  at  f becomes  (see  Table  I) 


1 n ^ 1 ^ 1 

£v*-d  ^ ^ 


e a-r,/ 
V d 


(145) 
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• ^(t) 


= e. 


■fe 


+ (0)  - 


^2 


+ eTX 


e &2 

V 


e,(d-r^) 

( -57  ) '"’2>  ■ “CT, 


e„(d-r.) 


(-I2  ) 


A 


A 

(-q^) 


V ‘‘  2 ' V 

Collecting  terms,  (145)  can  be  written  as 

. ■•d  . «r  . 

^ *r  ~.  'y  ~ 
d 


r(»  "V) 


(■7  ^ 

(-7^ 


i". 


r^d 


(146) 


(d-rd) 

«r  'r<'^-'d>  ^ 

‘iz 

(SzF 

(d-rd') 

®r 

Cr(d-rd)  ^ 

1 ^2  = 

(d-Td) 

”757)2“  J 

"d) 

Since  the  wires  are  assumed  to  be  sufficiently  separated  from  each 
other  and  the  ground  plane,  we  may  make  the  following  approximations: 

(147a) 
(147  b) 

do  = 41?+d2  (147c) 


dj  ^ 2h 

= 2h 


A 
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(147d) 


d-r  = d 

w 

1 

-Td  = O . 

Utilizing  these  relations  in  (143)  and  (146)  we  obtain  (V^^  = 0^) 

+ i^■)  ^2  * «”  (-f)  ^2} 

{ ('r  - »(^  ■ ^ <ll  + { (6r  - 1)  + ^ 

-{('r-  1)  (,|S)2  - -^)}  I2 

- |-)}  ““  • 


(I47e) 

(M7f) 

(148a) 

(148b) 


The  total  free  charges  on  the  conductors  are  given  by 

q = q + A 
fl  1 1 

q = q + q 
^f2  2 2 

Therefore  (148)  (along  with  similar  equations  at  match  points  2 

A 

and  2 ) may  be  written  as 

A B 


(14ya) 
(149  b) 


2Tre. 


and 


R]  ■ & ”]  fe]  • [S  c]  [y 

cj  [q^2J  b2J 

[:]  ■[•  .1  [::]  • t 


(150) 


G H 
H G 


(151) 
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H = F 


The  variables  and  can  be  easily  eliminated  from  equations 
(151)  since  D=B  and  H=F  and  the  result  is 


Therefore  the  entries  in  the  inverse  of  the  per-unit-length 
mission  line  capacitance  matrix  are  given  by 


C"^]..  = yi—  f C - G I 

~ “ ^rrsv  C (E-q  ) 


-2^  (■ 


I: 


(e  -1)  ''d 

Ih^ 


(152a) 

(152b) 

(152c) 

(152d) 

(152e) 

(152f) 

(152g) 
(152h) 
(150)  and 

(153) 

trans- 


(154a) 
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(154b) 


L_  / D - H 7 

( (E-G)  j 

■.-^{■■(f)‘“e(»(t))(S^)l 


for  i,  j=l,  2 and  d ''  is  given  by 


f = /liJ 


+ d" 


(155) 


For  this  case,  we  assumed  that  (1)  the  two  wires  are  identical, 

(2)  they  are  at  the  same  height  above  the  ground  plane  and  (3)  the  wires 
are  sufficiently  separated  from  each  other  and  the  ground  plane  such  that 
the  assumption  of  constant  charge  distribution  is  valid.  When  these  assump- 
tions are  no  longer  valid,  the  expression  for  the  entries  in  C“^  or  Q cannot 
be  easily  obtained  in  closed  form  and  a digital  computer  must  be  used. 

The  expressions  in  (154)  for  the  entries  in  will  be  used  in  a later 
publication  in  the  analysis  of  certain  experimental  data  for  which  this 
approximation  is  reasonably  accurate.  It  should  be  noted  that  an  approx- 
imation for  cj2  for  this  specific  example  has  been  obtained  in  [33]  although 
the  derivation  is  not  presented  and  evidently  relies  on  certain  empirical 
data. 
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VI.  SUMMARY 


A complete  and  unified  discussion  of  multiconductor  transmission  line 
theory  has  been  presented.  The  general  solution  of  the  problem  under  the 
assumption  of  TEM  mode  propagation  in  the  case  of  a homogeneous  medium 
or  quasi  - TEM  mode  propagation  in  the  case  of  an  inhomogeneous  medium 
has  been  presented  along  with  parameter  derivations  and  lumped-circuit 
approximations.  If  losses  can  be  neglected,  then  it  appears  to  be  as  efficient 
to  solve  the  transmission  line  equations  directly  and  incorporate  the  termina- 
tion networks  through  the  solution  of  (70),  (75)  or  (76)  as  it  is  to  use  lumped- 
circuit  iterative  approximations  described  in  Section  4.1.  The  matrix  chain 
parameters  for  the  distributed-parameter  approach  can  be  easily  obtained  in 
closed  form  suitable  for  numerical  computation  so  that  "abruptly"  nonuniform 
lines  can  be  handled  and  the  per-unit-length  parameters  must  be  obtained  for 
either  the  distributed-parameter  or  the  lumped-circuit  iterative  approach. 
When  solving  the  transmission  line  equations  directly,  one  is  not  required  to 
solve  an  increasingly- large  (although  sparse)  set  of  equations  for  increasing 
frequencies  when  the  line  is  not  electrically  short  as  is  required  with  the 
lumped-circuit  iterative  approximc- tions.  For  fhe  homogeneous- medium  case, 
a lossy  dielectric  can  also  be  included  with  no  additional  computational  dif- 
ficulties. Lossy  conductors  will  also  present  no  additional  computational 
problems  for  the  homogeneous-medium  case  if  the  n conductors  are  assumed 
to  be  identical.  When  losses  cannot  be  neglected  in  the  case  of  an  inhomo- 
geneous medium,  the  question  becomes  more  difficult  to  answer  since 
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diafj'unalization  of  YZ,  which  is  roquirod  in  a rlirnct  solution  of  the  trans- 
mission line  equations,  is,  in  general,  required  to  be  performed  at  each 
frequency  and  is  not  necessarily  guaranteed  except  in  the  case  of  cyclic 
symmetry  matrices  which  assume  certain  structural  symmetries  of  the 
line  as  described  in  Section  3.  3.  For  this  case,  it  may  be  preferable  to  use 
one  of  the  lumped-circuit  iterative  models  for  frequencies  where  the  line  is 
electrically  short,  e.  g.  , S-<  X,  and  approximate  the  line  with  only  one 
section,  i.  e.,  solve  (80)  with  N=l.  For  frequencies  such  that  £>  Vj^qX,  it 
may  be  preferable  to  solve  the  transmission  line  equations  directly  rather 
than  increasing  the  number  of  lumped-circuit  sections  to  approximate  the 
line  since  no  quantitative  criterion  for  determining  the  required  number  of 
sections  for  a given  approximation  accuracy  can  evidently  be  obtained. 
Numerical  techniques  for  obtaining  the  per- unit-length  parameters  for 
bundles  of  closely-coupled, dielectric- coated  wires  as  well  as  large - 
separation  approximations  for  wires  in  a homogeneous  medium  are  also 
given. 


APPENDIX  A 


The  purpose  of  this  appendix  is  to  demonstrate  certain  important  pro- 
perties of  the  TEM  mode  assumption  given  in  Section  II.  An  appropriate 
reference  for  these  results  is  [40]. 

The  first  objective  is  to  show  that,  for  the  TEM  mode  of  propagation, 
the  transverse  electric  field  vector  and  transverse  magnetic  field  vector 
satisfy  the  same  spatial  distributions  as  static  (DC)  fields  at  each  x along 
the  line.  The  electric  field  intensity  vector  and  the  magnetic  field  intensity 

vector  for  the  steady  state  and  sinusoidal  excitation  are  written  as 

t i'l)  t 

P(x,  y,  z,  t)  = E(x,  y,  z)  e-'  and  M(x,  y,  z,  t)  = H(x,  y,  z)  e respectively  where 

E(x,  y,  z)  = Ej^  X + Ey  y + Ej,  2 (A- la) 

H(x,  y,  z)  = X + Hy  y + z (A-  lb) 

andx,y,  zare  unit  vectors  in  the  x,  y and  z directions  respectively.  Assum- 
ing the  TEM  mode  of  propagat’  in,  E^  = = 0,  the  field  vectors  are  entirely 

transverse  to  the  x direction  and  are  denoted  as 

E'j'(x,  y,  z)  = Ey  y + z (A-2al 

H^(x,  y,  z)  = Hy  y + Hg  z . (A- 2b) 

Now  consider  the  general  (n+D-conductor,  uniform  transmission  line  in 
Fig.  la  :onsisting  of  (n+1)  perfect  conductors  in  a linear,  isotropic  and 
homogeneous  medium.  Faraday's  law  and  Ampere’s  law  become  for  the 
TEM  mode  of  propagation  (in  the  source-free  medium) 
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vX  E<p  - - jjuU  H'p 
vX  Hj  = (cf  + ju>e)  E-p 


(A.  3a) 
(A- 3b) 


where  the  medium  may  be  lossy  through  the  effective  conductivity  a which 
includes  ohmic  conductivity  and  dipole  relaxation  losses  and  e refers  to  the 
real  part  of  the  complex  permittivity. 

Separating  the  curl  operator  into  a transverse  and  a longitudinal  com- 
ponent as 


7 = 


7'j' 


(A. 4) 


and  applying  to  (A- 3)  we  obtain 


7-p  X E,p  +-^(x  X E^) 
X % +-~(x  X H^) 


However,  X E,p  and  V,pX  are  vectors 
fore  we  have  by  matching  components 


= - j(i)U  H.J, 

■4 

= (a  + ju)  e)  E,p  • 

in  the  x direction  only. 


(A- 5a) 
(A. 5b) 
There- 


-♦ 

X E-p  = 0 

(x  X E^)  = - ja,U  Hp. 
V-p  X Hp  = 0 

-^{5^  X Hp)  = (a  + j(i)e)  Ep 


(A- 6a) 


(A-6b) 


(A-6c) 


(A-6d) 
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Note  also  from  v*D  = p and  7*  B = 0 so  that 


V 

vH-j.  = V-p*  H-p  = 0 


(A-7a) 

(A-7b) 


(since  E^  = = 0)  where  p is  the  free-charge  density  in  the  surrounding 

medium  (which  will  decay  to  zero  with  time  cons^^ant  e/a).  Therefore, 
equations  (A-6a),  (A-6c),  (A-7a)  and  (A-7b)  show  that  the  transverse  field 
vectors  E-p  and  H<p  satisfy  the  same  spatial  distributions  as  static  fields  in 
any  transverse  plane  (y,  z)  at  each  x along  the  line. 

This  may  be  more  easily  seen  if  we  write  Faraday's  law  and  Ampere's 
Law  in  integral  form  by  applying  Stokes'  theorem  to  V XE  = - ju)UH  and 
VX  H = (a  + jdi  c)  E +•  J as 


\ H • da 

•^s 

da  + (a  + j(ju  e)  C E 

■s 


(A-8a) 
da  (A.  8b) 


where  C is  a closed  contour  enclosing  the  open  surface  S.  Taking  C to  be  a 
contour  in  the  transverse  (y,  z)  plane  denoted  by  S to  be  a flat 

surface  in  the  transverse  plane  denoted  by  which  is  bounded  by  then 
(A- 8)  becomes 

^ (Ey  dy  + Eg  dz  4-  Ejj  dx)  = - jmM  ^ dy  (A- 9a) 

^yz  Syz 

£ (H  d + d^  dx)  = V Jxd  d^  + (a  + joue)  \ E^  d dz  (A-9b) 

p ' 'c  c ^ 

^yz  *^yz  “^yz 
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v/here  is  any  source  current  in  ^he  x direction  penetratir.a  Sy  v.  Hav”c\oi 


undci-  the  TEM  mode  assumption  = 0 and  (A- 9)  becoi.ios 

^ (Ey  f E^  d^)  = ^ E T • (z  d z + y dy)  = 0 

Cy2  ^VZ 

§>  (Hy  d + Hz  dz)  = ^ Ht  • dz  + ^ dyr)  = ^ d,,  d, 

C-vz  ■'  '-'vz  oyz 


(A-  lOa) 


(A-rOb' 


which  are  of  course  indicates  that  and  are  no  longer  coupled  logeMier 
>.  s is  the  case  for  static  electromagnetic  fields. 

Thus  we  may  uniquely  define  the  voltage  of  the  i-th  conductor  and  liu' 


"urrent  associated  with  the  i-th  conductor  as 


V.(x)  = - E^  * d? 

IJ^(X)  = ftrp  • df, 


(A-  llu) 


(A -lib) 


'-'.  here  Cj^  and  are  shown  in  Fig.  lb  as  contours  in  the  trans-ce  r.‘'.e  pi  i.nc-  e, 
c.  particular  x along  the  line. 

The  second  objective  is  to  demonstrate  that  for  (r-^1)  perlecT  conriuc't  ■ ; 
•L  hc.  "nogeneous  nnediurn,  Die  per-uni;-'c,ng:h  indu-.  laK'.  u .".  'v  ' , . ■ 

•'  . '.h  ca ;vi c i la nco  natri'-,  C,  and  • 

scldsly  the  important  relations  in  l i-i).  T-'n-.'v  f,\  . c...  > ' 

, obtain  by  taking  Ihe  pariiai  c!  > \-a  I ' v.'  r>' 

>■  and  taking  the  curl  of  x with  each  equation. 


(A- 12a) 


- (JU®  U€) 


9P 

3 E»j«  = (jujUcj-  (I)®  u e)  Ey 
ox 


From  equation  (19)  assuming  perfect  conductors 


(A- 12b) 


I(x)  = YZ  I(x)  = (jujGL-  .u^cy  I(x) 


d» 


^ V(x)  = Z Y V(x)  = (jujLG  - tr^LC)  V(x) 


(A- 13a) 
(A- 13b) 


Performing  the  operations  indicated  in  (A- 11)  on  (A-12)  one  obtains  two  sets 
of  n equations 


d8 


li  (x)  = (jujucr  - (w’^ue)  I^(x) 


d2 


V.  (x)  = (jouUCT  - tt)®u  e)  v^(x) 


(A- 14a) 
(A- 14  b) 


Arranging  these  for  i = 1, ,n  as  in  (A-13)  shows,  by  matching  real  and 

imaginary  parts,  that 

LG  = GL  = ucrl„  (A-15a) 


LC  = C L = u € 1. 


(A- 15  b) 
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APPENDIX  B 

The  purpose  of  this  appendix  is  to  demonstrate  the  derivation  of  the 
multiconductor  transmission  line  e<itiations  in  (4),  (6)  and  (7)  from  the  per- 
unit>leugth  lumped  equivalent  circuit  in  Fig.  7. 

Utilizing  Kirchoff's  voltage  law  counter-clockwise  around  the  loop  con- 
sisting of  the  i-th  conductor  and  the  zero-th  conduct^or  in  Fig.  7 yields 

V.(x  + Ax)  - Vg.(x)  Ax  + jai  f i Ax  juun-i^k  ^ 


+ (r^  Ax  + ju)  X Ax)  - Vi(x) 
i i 


Ax  t ju)  Ax)  Iq  - j{«io  Ax  Iq 


Jl  AX  Ik  = 0 


This  equation  can  be  rewritten  as 


V.(x  + Ax)  - Vj(x) 


Ax 


= - joj^i  Ij  - 


kj^i 


(B-2) 


+ jujio  ^0  jiw^'kO  ^k 

k=l 


+ (X) 

1 
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and  the  current  in  the  reference  conductor  satisfies 


I 


0 


n 

E 

k=l 


(B-3) 


Substituting  (B-3)  into  (B-2)  one  obtains 
Vi(x  + ^x)  - Vi(x) 

^ li  - 


(B.4) 


^iO  + «^il  4 


(""cq  + j«’4:o  * h-  ***^  n) 

+ ju) (nijQ  Ij  + + ^n) 


+ Vgjx) 

which  may  be  rewritten  as 
Vi(x  + AX)  - V£(x) 

Ax  “ “ ^*’cq  ^ ^ ■ juD^^lO  “ ju)J^iO  + juum^j)  Ij  (B-5) 


- (^C£  + i'i>^c£  + j«^^i  ^CQ  + juu-^O 


- juun^£Q  - jojnn.Q)  l£  - • . . 


•••  - (^cq  + j«>^CQ 
+ Vs.(x) 

Arranging  these  equations  for  i=l, ,n  and  taking  the  limit  as  Ax  -»0  yields 

(4q)and  the  per-unit-length  impedance  matrix,  Z,  can  be  separated  as  in 


(6a)  with  the  entries  given  in  (7a),  (7b).  (7c). 

The  derivation  of  the  second  transmission  line  equation,  (4b),  proceeds 
similarly.  Utilizing  Kirchoff's  current  law  for  the  i-th  conductor  in  Fig.  7, 
we  may  write 

I|(x  + 6x)  = Ij(x)  - (g^Q  t^)  V.(x  + ^)  (B-6) 

" ^®^ik  ^ A*HV^(x  + 6x)  - Vj^^(x  + t^)) 

1C  "“A 

k^'i 

+ ^ 

which  may  be  rewritten  as 
I;(X  + Ax)  - I:(x) 

— = (gil  + jujc.^)  (X  + Ax)  (B-7) 

n 

+ • • * - (gjQ  + jaiCio  + ^ (gik  + jujc.^^))  V.(x  + Ax) 

+ • • • + ^gin  + ''^n  (^  + ^)  + 

Arranging  for  i=l, ,n  and  taking  the  limit  as  Ax  0 yields  (4b)  and  the 

per-unit-length  admittance  matrix,  Y,  can  be  separated  as  in  (6b)  with  the 
entries  given  in  (7d)  and  (7e). 
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APPENDIX  C 


The  purpose  of  this  appendix  is  to  derive  the  expressions  for  the  per- 
unit-length  equivalent  sources  in  (5)  which  are  induced  by  an  incident  elec- 
tromagnetic field  on  the  uniform  transmission  line  in  Fig.  2a  consisting  of 
(n+1)  perfect  conductors  in  a lossless,  homogeneous  meuium  (see  Fig,  C-1). 
The  expression  for  the  currents  induced  in  termination  networks  given  in 
(77)  will  also  be  obtained.  The  solution  for  Fig.  2b  is  also  discussed. 

The  solution  for  the  special  case  of  a two- conductor  line  (n=l)  was 
obtained  by  Taylor,  Satterwhite  and  Harrison  in  *^20]  and  later  in  a more 
convenient  form  by  Smith  in  *^21].  The  solution  for  the  case  of  a uniform 
plane  wave  incident  on  a three- conductor  line  in  the  transverse  direction 
(perpendicular  to  the  system's  longitudinal  (x)  axis)  with  the  electric  field 
intensity  vector  polarized  parallel  to  the  line  axis  was  obtained  in  fZ4], 
Procedures  for  extending  this  result  to  multiconductor  lines  were  indicated. 

It  is  convenient  to  consider  the  effects  of  the  spectral  components  of 
the  incident  field  as  per-unit -length  distributed  sources  along  the  line  [26], 

The  sources  appear  as  series  voltage  sources  and  shunt  current  sources  as 
indicated  in  Fig.  C-2  for  an  "electrically  small"  6x  section  of  the  line.  The 

multiconductor  transmission  line  equations  may  then  be  derived  for  the 
Ax  subsection  and  are  given  in  equation  (4).  The  termination  networks  are 
given  in  the  form  of  generalized  Thevenin  equivalents  as  in  equation  (69)  and 
the  solution  for  the  termination  currents  is  given  in  equation  (75).  Substiba- 
ting  the  expressions  for  the  matrix  chain  parameters  given  in  (49)  into  (75) 
for  this  case  of 


-133- 


perfect  conductors  in  a lossless,  homogeneous  medium  (free  space) 
described  by  e and  U , we  obtain 

fco.(BX){zo  + Zj,}  +j»in(e«  {zc*2x5c‘'5o}li«»= 

- + ij  «in(8i)  Zj  Zc  ‘ 

+ i,<«  - 5xi.<« 

I(X)  = - j «ln(BX)  Z£.-‘  y,  + fco»(BX)  1^^  t j sm(BX)  gc  ‘ Zq]_I(0)  (C-lb) 
+i,(X) 

where  the  wave  number  is  0 = Zn/\,  \ = y/f,  y = 1//H^  '=  3 x 10®  m/sec. 

The  characteristic-impedance  matrix  Zq  becomes 

Zc  = oL  (C-2) 


A A 

and  Vjj(X)  andj^(£)  are  obtained  by  substituting  (49)  into  (74)  as 


A 

V^(X)  = 

{co8(e(X- :^))  yg(x)  - j sin{g(X-x))  Z<3_Ig(x)  } 

dA 

(C-3a) 

A 

= \ • 

[cos(0  (X- x))_Ig(x)  - j sin(6(X-x)) 

} dx 

. (C-3b) 

Solution  of  (C-la)  for  the  termination  current  vector,  _I(0),  requires  the 
solution  of  n complex  equations  in  n unknowns  (Ij^(O)).  Once  (C-la)  is  solved, 
(C-lb)  yields  the  termination  currents  I(£)  directly. 

Although  the  equations  may  appear  formidable,  they  arc  in  a compact 
form  and  can  be  straightforwardly  programmed  on  a digital  computer. 
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Furthermore,  the  form  is  not  restricted  to  any  particular  value  of  n.  The 

A A 

only  difficulties  are  in  determining  L and  determining  and_Ig  (which 
require  that  we  determine  Vg(x)  and_Ig(x)).  The  determination  of  the  equiva- 
lent sources  V.  (x)  and  I.  (x)  induced  by  the  incident  field  will  be  the  next 

Si  8^ 

objective. 


C,  1,  Determining  the  Equivalent  Induced  Sources 


In  order  to  determine  the  equivalent  induced  sources,  V„  (x)  and  I (x), 

s . s • 

L I 

consider  Figure  C-1,  The  method  used  in  [20]  can  be  adapted  here  in  a 
similar  fashion.  Faraday’s  law  in  integral  form  becomes 


fc  E • di*  = - ^ H • n 


da 


(C-4) 


where  Is  a flat,  rectangular  surface  in  the  x,  y plane  betv/een  wire  i and 
wire  0 and  between  x and  x + (\x  as  shown  in  Figure  C-1.  The  unit  normal  n 
is  n = z where  z is  the  unit  vector  in  the  z dix'ection,  da  = dx  dy  and  is  a 
contour  encircling  in  the  proper  direction  (counter-clockwise  according  to 
the  right-hand  rule).  Equation  (C-4)  becomes  for  the  indicated  integration^ 

,d 

\ 

0 


r iO 

\ [Ej,  (y,  x+^pc)  - Ej  (y,  x)]dy 


(C-5) 


r 


x+^ 


- \ [Ejj  (d^Q,x)  - E,  (0,x)]dx 


X: 


X + iPC  d.Q 

•v  V i 


1.  In  integrating  from  y=0  to  y;:djQ,  we  are  implicity  assuming  that  r^^  and 


wi 


rwo  much  less  than  d^Q,  i.e. , the  wires  are  sufficiently  separated  so 
that  they  may  be  replaced  by  filaments. 
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where  is  the  component  of  the  total  electric  field  (incident  plus  scattered) 
i 

transverse  to  the  line  axis  and  lying  along  a straight  line  joining  the  two 
conductors,  i.e.,  E,  = E ; E,  is  the  component  of  the  total  electric  field 

h y h 

along  the  longitudinal  axis  of  the  line,  i.  e. , S . = E ; and  H io  the  com- 

* ”i 

ponent  of  the  total  magnetic  field  perpendicular  to  the  plane  formed  by  the 
two  wires,  i.e.,  =H  . 

z 

Defining  the  voltage  between  the  two  wires  as 


r iO 

V.(x)=-\  E^(y,x)dy 
0 ^ 


(C-6) 


dVj^(x)  . 

= lirn  \ FEj,  (y,x+Ax)  - E^  (y,  x)]dy  (C-7) 

Ax-»0  ^ ^ i i 


The  total  electric  field  along  the  wire  surfaces  is  zero  since  we  assume 
perfect  conductors.  (One  can  straightfor-wardly  include  finite  conductivity 
conductors  through  a.  surface  impedance  as  wh&  done  in  [20]).  Therefore 
(C.  5)  becomes  in  the  limit  as  Ax-*  0 


dVi(x) 


r^io 

= \ Hj^,(y,  x)dy 


(C-8) 


The  total  magnetic  field  is  the  sum  of  an  incident  and  a scattered  field: 


H (y,,x)  = H (y,x) 
n.  z 

(scat) 


(inc) 


= H^(y,x)  K,j(y,  x) 


scattered  incident 


(C-9) 


/'iTCTJ- 


m 


HUB 


and  the  scattered  field  here  is  produced  by  the  transmission  line  currents. 
The  scattered  flux  passing  between  the  two  conductors  per  unit  of  line  lenj. 
is  directly  related  to  the  scattered  magnetic  field  and  the  per- unit -length 
inductance  matrix,  L,  as 

(scat)  ^dj^Q  (scat) 


' p iO  ' ' 

ai.(x)  - U^H„^(y,x)dy 


f^il’  ^i2'  ' ^n^ 


Ij^(x) 

I?(x) 


In(^) 


(C-10) 


v/here  Substituting  (C-10)  and  (C-9)  into  (C-8)  and  arranging  for 

IJ  ^ Ij 

i = 1,  , n yields 


V(x)  + j(ULI(x)  = 


JO)U, 


d^Q  (inc) 


r iO  ' — ' 

\ H„.(y,x)dy 


(C-11) 


and  the  source  vector  V (x)  in  (4)  is  easily  identified  by  comparing  (4a)  and 
(C-11). 

For  transmission  line  theory  to  apply,  the  cross-sectional  dimensions 
of  the  line  (wire  spacing,  etc.)  must  be  electrically  small,  i.  e.  , 3dj^0  <<  1 
and  pd..  <<  1,  Thus  the  result  indicates  that  the  voltage,  , induced  in 
the  loop  between  the  i- th  conductor  and  the  zero-th  conductor  and  betv'een  x 
and  x+Ax  is  equal  to  the  rate  of  change  of  the  incident  flux  penetrating  this 
"electrically  small"  loop  which,  of  course,  makes  sense. 
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npeve's  law  yields 


(C-12) 


E will  consist  of  scattered  and  incident  field  components  and  is  written  as 


Ej.  (y,  x)  = E (y,  x) 
i ^ 


(scat) 

= E (y,  X)  + 
scattered 


(inc) 

Ey(y.x) 

^ — >■  • 

incident 


(C-13) 


Substituting  (C-12)  into  (C-6)  we  have 


(scat) 


e^iO  1 e io/  ^ 

Mx)  = - V E„(y,x)dy  = \ < — 


(inc) 


(scat) 


x)  9 *) 

— + 

9x 

(inc) 


- 9 ^x^y»  S H 


92 


(C-14) 


9 H (y,  x)  1 
t dy  . 

9z  J ^ 


Utilizing  (C-10)  we  obtain 


V.(x)  = - 


jU)U  € 

V V 


d 

dx 


\z 


, iscat) 

iO  9Hj^(y,x) 

9 z 


dy 


,d.Q  (inc) 


5 

0 ^ 


(y»  x)dy 


(C-15) 


If  we  assume  that  the  currents  on  the  wires  are  directed  only  in  the  x 

direction  i.  e.  , (there  are  no  transverse  components  of  the  currents  on  the 

(scat) 

wire  surfaces),  then  Hj^(y,  x)  = 0 and  (C-15)  becomes 
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V.(x)  = - -1. 

^ Juu^vSr 


{ 


^12' ' ^in 


]l(x)| 


(C-16) 


c^io 

- \ E.(y,x)dy 

'o  ^ 

Arranging  these  equations  for  i = 1, , n we  obtain  the  second  transmission 

line  equation 


j(x)  + L-W(x)  = - jUJU^c^  L-^ 


d.Q  * (inc) 


r io 

E^  (y,  x)dy 


(C-17) 


Utilizing  the  important  relation  for  a homogeneous  medium,  C = L”^  in 

(C-17)  we  obtain  by  comparing  (C-11)  and  (C-17)  to  (4) 


Vs(x)  = j(i)U 


_Is(x)  = -jojC 


iO 

0 . 


(inc) 

H„.(y,  x)dy 

i 


(C-  18a) 


<*10  ■ 


Ej,  (y,x)dy 

• i 


(C-  18b) 


The  shunt  current  sources  in  I (x)  are  therefore  a result  of  the  line  voltage 
induced  by  the  incident  electric  field  being  applied  across  the  per-unit-length 
line- to- line  capacitances  which,  of  course,  satisfies  our  intuition. 
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A A 

C.  2.  Solution  for  V (x)  and  I (x) 

— *8 

A 

The  final  problem  rernaining  is  to  obtain  simplified  versions  of  V and 

A 

in  (C-3)  to  be  directly  used  in  (C-1).  First  consider  the  determination  of 
V (£).  Substituting  (C-18)  into  (C-3a)  yields 

«w3 


V^(X)  = jlDUv  \ t <=0^(0 
0 


I 


4o  * * 

x)<3y 


) 


dx  (C-19) 


- B 


ol 


sin(B(X-  x)) 


d^O  ‘(i«c) 

\ E^(y.x)dy 

’O  . i 


) 


dx 


From  Faraday's  law  we  obtain 


(inc) 


1 


= 

"i  jdiUv 


[ 


(inc) 

BE. 

dy 


(inc) 

aEf  I 

^ J 


(C-20) 


Substituting  this  into  (C-19)  yields 


* cV 

= y 


cos(0  (X-  x)) 


(inc)  . (inc) 

^1.  (0'^) 


} 


dii 


■U 


cos(e  (X-  x)) 


, .(inc) 

'""iO  3Et.(y,Ji) 

\ — dy 

0 Sx 


} 


dx 


(C-21) 
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- 0 V <r 8m(p(X-x)) 


\ Ej  (y.  x)  dy 

•^0  . i 


) 


dx 


Utilizing  Leibnitz's  rule, (C- 21)  is  equivalent  to 


■y 


V„(X)  = \i  eo8(p(X-x)) 


(inc)  . (me) 

Ejj  (d.Q,^)  - (0,5c) 


] 


dji  (C-22) 


pX 

\ i CO8(0(X-X)) 

0 ^ 


-d.Q  -(inc) 


_iQ  

Et.(y,^)  dy 


} 


di 


and  this  may  be  written  as 


Vg(X)  = 


^ jcos(B(X-x)) 

0 ^ 


(inc) 

Ef.(dio.i) 


(inc) 

Ejj^(0,x) 


(C-23) 


• 

>»» 

• 

• 

dio  (inc) 

• 

dio  (inc) 

\ E*  (y,  X)  dy 

+ cos(0  X ) 

\ E^(y,0)  dy 

•^0  . i 

0 . i 

• 
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Similar ily_Ig  (X)  may  be  obtained  as 


is(D  = 


\{ 


am(p(£-  x)) 


(inc)  . (inc) 

*)-  E£.(0,x) 


1 


) 


dx  (C-24) 


\ E^(y,0)  dy 

0 * i 


} 


Th.  important  quantity  in  (C-U)  i.  V.,«  . £ combining  (C-23) 

and  (C>24),  this  becomes 


ls<^)  = ^ f>'co8(B(X-.x))  + j sin(B(X- 

qI  ~n 


X 


*>  5x  (^-25) 


• 

• 

- 

• 

• 

(inc). 

E,.(dio.^)  - 

• 

• 

(inc) 

E^  (0.x) 

) 

• 

• 

d.Q  (inc) 

\ (y,  X)  dy 

0 . i 

• 

• 

-144- 


+ [cos(b£)  1 + j 8in(p£)  Zp  Z^”^] 

Note  that  the  equivalent  forcing  function  on  the  righthand  side  of  (C-la), 

A A 

- 5x  given  in  (C-25)  is  simply  determined  as  a convolution  of 

differences  of  the  incident  electric  field  vector  along  the  wire  axes, 

Cine)  (inc) 

*^nd  a I'near  combination  of  integrals  of  componaits  of 

the  electric  field  vectors  at  the  endpoints  of  the  line  which  are  transverse  to 
(inc)  (inc) 

the  line,  E^^(y,  £)  and  (y,  0),  This  is,  of  course,  precisely  the  result 

obtained  by  Smith  [21]  for  two  conductor  lines.  Substituting  (C-25)  into  (C-la) 

and  setting  Vp=  0 , = 0 , i.  e* , no  independent  sources  in  the  termina- 

* n-1  ^ n-1 

tion^networks,  one  can  verify  that  the  result  reduces  for  two- conductor  lines 
(n=l)  to  the  result  given  by  Smith  [2l]  since  Z^,  Z_,  Z,^  become  scalars  for 
two-conductor  lines  and  (C-la)  becomes  one  equation  in  only  one  unknown 
1(0).  For  uniform  plane  wave  illumination  of  the  lire  (which  is  usually  the 
case  of  interest),  (C-25)  reduces  to  a much  simpler  form  although  the  result 
allows  for  the  more  general  case. 

The  final  equations  for  the  line  currents  then  become  (substituting  (C-25) 
into  (C-D) 

fcos(8r){zQ  + Z^  } +jsin(p£)  {gc  + ^0  ^ = 
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(C-26a) 


I 


- 'j  sin(8X)  Z^  + cos(pX)  1,^] 


-10 


+ (rcos(8(X-  X))^„  + j 8Ln(8(X-:^))  Z^ 


0 

(Lnc) 


- E,(X)  + {rcos(8X)J.^  + j sin(BX)  Zj.  Z^-^]  E^(o/  | 


lU)  - - j sin(8i)  Z^’^  Vq  4 [cos(8£)^^  + j sm(8JC)  ^ Zjj]_I(0)  (C-26b) 

1 r * 1 

- j^c’  {sm(8(X-x))  E^$c)  jd5c 

.1  r 

- J^c  {8in(8£)  Et(0)j 

(inc)  (inc)  (inc) 

where  E^(x),  _E^(X)  andJE^(O)  are  n -y  I column  vectors  with  the  entries  in  the 


i-th  rows  given  by 


(inc) 


(inc) 


(C-26c) 

1 i 

P^iO 

(C-26d) 

/iO 

V Ej  (5i,0)  d5 

>0  1 

(C-26e) 

(inc) 

E^W 

(inc) 

Et(X)j 

(inc) 


for  i = 1, , n which  are  equations(77). 

A word  of  caution  in  the  interpretation  of  the  notation  is  in  order. 

Although  it  should  be  clear  from  the  derivation,  the  reader  should  neverthe- 

(inc) 

less  be  reminded  that  the  integration  path  for  the  component  E^.  is  in  the  y 
direction  when  the  i-th  conductor  is  concerned.  When  other  conductors  are 
concerned,  the  integration  path  is  a straight  line  in  the  y,  z plane  which 
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joins  the  conductor  and  the  zero-th  conductor  and  is  perpendicular  to  these 
two  conductors.  This  is  designated  as  in  (C-26)  and  replaces  the  y vari- 
able for  the  path  associated  with  conductors  i and  0.  The  notation  may  be 
cumbersome  but  the  idea  and  implementation  is  quite  simple. 


C.  3.  The  Per-Unit-Length  Inductance  Matrix,  L 


One  final  calculation  remains;  the  determination  of  the  per-unit-length 
inductance  matrix,  L.  Ordinarily  this  is  a difficult  calculation  as  discussed 
in  Section  V.  However,  if  we  assume  that  the  conductors  are  separated 


sufficiently  such  that  the  charge  distribution  around  the  periphery  of  each 
conductor  is  constant,  then  the  conductors  can  be  replaced  by  filamentary 


lines  of  charge.  Typically,  this  will  be  quite  accurate  if  the  smallest  ratio 
of  conductor  separation  to  wire  radius  is  greater  than  5 fS6].  In  this  case, 


the  entries  in  L are  shown  in  Section  V to  be 


wi  wO 


fL]..  = U e [C-^]  . = ^ /n  ) 

vv~  2TT  . 


(C-27a) 


(C-27b) 


For  closer  conductor  spacings,  proximity  effect  will  alter  the  charge  distri- 
butions from  constant  ones  and  numerical  appx’oximations  must  be  employed 
to  find  L as  was  discussed  in  Section  V. 

The  entries  in  the  per-unit-length  inductance  matrix 
for  large  wire  spacing  given  in  (C-27)  can  be  derived  in  an  alternate  manner 
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which  more  clearly  illustrates  its  relation  to  the  total  scattered  flux  passing 

(scat) 

between  the  wires.  The  matrix  L.  relates  the  scattered  flux  ^ passing 


between  the  wires  to  the  wire  currents  as 


(scat) 

0 


(scat)*" 

0l 

• 

ill  . . . 

• 

• 

^In 

• 

• 

h 

• 

• 

(scat) 

L^n  J 

# • • • 

L^nl 

M 

^nn 

I 

n 

m m 

ii 


(C-28) 


The  respective  entries  are  determined  as 
(scat) 

0- 

je..  = ^ 


“ Ii  'Ii, , --- 


I =0 
’ n 


(scat) 

^ I 

i^j 


Ij.l'  Vl = ° 


(C-29a) 


(C-29b) 


and  Large  wire  sep>aration8  are  assumed  so  that  the  wires  may  be 

replaced  by  filaments  of  current. 

Consider  Fig.  C-3a.  The  magnitude  of  the  magnetic  field  intensity 
vector  due  to  on  wire  i at  a distance  r > away  from  wire  i is 


H 

r 


h 

Zrrr 


(C-30) 


and  the  total  flux  passing  between  wire  i and  wire  0 due  to  both  currents  is 
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(scat) 

<t). 


* 'wL  J 


(C-31) 


= la  (_£iol__) 

'r  • T * 

^wi  wO 


2n 


Thus  £..  is  easily  identified  as  in  (C-27a). 


11 


(scat) 

Consider  Fig.  C-3b.  The  portion  of  the  flux  passing  between 

wire  i and  wire  0 due  to  -Ij  in  the  reference  conductor  is  as  above 


2TT  ^ 


(C-32) 


and  the  portion  of  the  flux  passing  between  wire  i and  wire  0 due  to  I.  in  the 


j>th  conductor  can  be  found  to  be 


(scat) 


(6 


V H dp=-^  I V r 

0 “ 277  3 I \ Q 


‘0  dp  1.  (C.33) 


■} 


= ^x,fjL„  r <1^11 

I " L ?o  " ^^iO  - J ' 


Combining  (C>32)  and  (C-33)  we  obtain 


(scat)  (scat)  (scat) 


= <t> 


iO 


+ <6 


kiv’^ 


u 


2tt 


j-  In 


(C-34) 


Since 


“ij'  = ?o'  " '“io  - 


(C-35a) 
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. 2 » 2 j.  ^ < 

jO  ^0  ^ PQ 


(C-3Bb) 


and  is  easily  identified  as  given  in  {C-27b). 


C.4,  Computed  Results 

To  show  the  simplicity  of  the  result  and  indicate  its  equivalence  to  the 
result  obtained  by  Harrison  in  [24],  Example  1 considered  in  ['24]  will  be 

_3 

computed  by  this  method.  Three  wires  all  of  radius  10  m lie  in  the  x,  y 
plane  as  shown  in  Figure  C- 4.  A uniform  plane  wave  with  an  electric  field 
intensity  magnitude  of  IV/m  is  propagating  in  the  y direction  and  5000 
(purely  resistive)  loads  connect  each  line  to  conunon  nodes.  The  various 
distances  in  Figure  C-1  are  dj^g  = 10”^m,  d2g  = 2 x 10  and  dj2  = 10  ^m. 


Zg  and  can  be  easily  shown  to  be 


1 

5o  " 


1000 

500 

_500 

1000 

] . 


The  characteristic-impedance  matrix,  using  the  values  for  the  per-unit* 
length  inductance  matrix  given  in  (C-27),  becomes 

fin(lOO)  /n(20)  1 


Z_  = uL  = 60 


tn(ZO) 


In(ZO)  ' 
/n(400)J  , 


(inc) 

= 0 in  (C-26)  and  the  electric  field  intensity  of  the  wave  is 

(inc) 

E,(y,x)  = E e'Jfiy 

i X 
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(I,(0)+I,(0)) 


Uniform  Plane 
\Htove  Propagating 
in  + y Direction 


r»i 


m 


Ex-IV/m  e’^^y 


Figure  C-4.  The  geometry  for  the  example. 


From  this,  one  can  determine 


(inc)  (inc)  -oin-Z 

E^^(djQ,x)  - E^^(0,x)  = .1 

Two  frequencies  are  considered  in  [24]  in  terms  of  p£;  p£=  1.5,  pX=  3,0. 
Equations  (C-26)  with  the  above  items  were  programmed  on  an  IBM  370/165 
computer  in  double  precision  arithmetic.  The  execution  time  (cpu  time) 
was  .01  seconds  (1/100  sec)  and  the  results  are 


pX  = 1.5 


’^1Iq(0)I=  1.7662556E  - 5A 
‘ Il^(0)I  = 9.0756083E  - 8A 

]l2(0)l=  1. 7671218E  - 5A 


/Iq  (0)  = 70.77* 
/ll(0)  = -13.9“ 
/l2(0)  = -109.  52“ 


0 X=  3.0 


I 


ll  (0)  I = 5.4543875E  - 5A 
0 

Il^(0)I  = 7.7363155E  - 7A 
Il2(0)l=  5.4608I10E  - 5A 


/Iq(0)  = 9.845“ 
(0)  = -75.8“ 
/l2(0)  = -170.96“ 


The  computed  results  obtained  by  Harrison' 


method  and  given  in  [24]  are 


f [1^(0)  I = 1.766E  - 5A 
llj(0)|  =9.076E  - 8A 
1X2(0)  I = 1.767E  - 5A 
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'TTvfti'Jur 


B£=  3.0 


r IIq(0)1  = 5.454E  - 5A 
" ” (0)1=  7.736E  - 7A 


]l2(0)l  = 5.461E  - 5A 


I 

Sf; 


The  results  computed  by  this  method  are  exactly  those  computed  by 
Harrison's  method  in  [24],  However,  with  this  method  only  2 simultaneous 
equations  in  the  2 unknowns,  1^(0)  and  l2(0)»  ^.re  required  to  be  solved 
(Iq(0)  = -lj(0)  - 12(0)).  Harrison's  method  required  the  solution  of  10  simul- 
taneous equations  in  10  unknowns.  Furthermore,  Harrison's  method  was 
restricted  to  uniform  plane  wave  illumination  of  the  line  with  the  wave  inci- 
dent perpendicular  to  the  line.  Since  = Zq  for  this  example  and  since  the 
uniform  plane  wave  is  propagating  broadside  to  the  line,  1(0)  = I (£). 

C-5.  Extension  of  the  Method  to  Wires  Above  a Ground  Plane. 

Consider  the  system  of  n wires  in  free  space  above  an  infinite 
ground  plane  shown  in  Fig.  2(b).  The  result  for  (n+1)  wires  given  in  (C-26) 
can  be  extended  to  this  case  with  the  following  observations.  Consider  Fig. 
C-5.  Clearly  we  may  apply  Faraday's  law  in  (C-4)  and  the  previous  devel- 
opment to  the  flat, rectangular  surface  in  the  x,  y plane  shown  in  Fig,  C-5b 
between  the  ground  plane  and  the  i-th  wire  and  between  x and  x+^.  This 
flat, rectangular  surface  Sj^  lies  in  the  x,y  plane„  Equations  (C-26a)  and 
(C-26b)  will  again  be  obtained.  Equations  (C-26c),  (C-26d)  and  (C-26e)  be- 


come for  this  case 


(inc)  (inc)  (inc) 

[E^(x)]i  = Ej,.  (h.,x)  - E^.  (0,x) 


(C-36a) 
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J 


(C-36b) 


(inc)  (inc) 

E,.  (?i.  i)  d?i 


(inc)  ^hi  (inc) 

(0)]i  = \^  Et.  (?i,0)d§. 


(C-36c) 


where  is  a straight-line  contour  in  the  x,  y plane  between  the  position 

of  the  ground  plane,  y=0,  and  the  i-th  wire  which  is  perpendicular  to  the 

ground  plane,  i.e.  = y.  This  is  indicated  in  Fig.  C-5a. 

(inc) 

E..  (h;,  x)  is  the  component  of  the  incident  electric  field  parallel 

r,i  l 

(inc) 

to  the  axis  of  the  i-th  wire  at  y=h^  and  E^.  (0,x)  is  the  component  of  the 

incident  field  parallel  to  the  ground  plane  directly  beneath  the  i-th  wire, 
^(inc)  component  of  the  incident  electric  field  parallel  to  5^  and 

directed  in  the  +y  direction. 

The  per- unit-length  inductance  matrix,  L,  can  be  obtained  in  a 
fashion  similar  to  Section  C,  3 by  determining  the  scattered  magnetic  flux 
passing  through  the  surface  between  the  i-th  wire  and  the  position  of  the 
ground  plane^the  ground  plane  is  replaced  by  image  wires)  and  is  given 
in  (122)  as 


for  i,  j=l, , n where 

d.r  = /d.  /'  + 4h;h. 
ij  V U 1 J 


(C-  37a) 
(C-37b) 


(C-38) 
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APPENDIX  D 


The  purpose  of  this  appendix  is  to  derive  the  simplified  form  of  the 
nodal-admittance  matrix  for  the  lumped-circuit  iterative  approximations 
given  in  (80)  and  the  forms  of  the  rnatrix  chain  parameters  given  in  (78). 

To  derive  the  first  equation  in  (80)  consider  the  equations  for  the  termin- 
ation-networks in  (69): 

V(0)  = - ZqJ[(0)  (D-la) 

VU)  = + Z^_I(£)  (D-lb) 

and  the  matrix  chain  parameter  equations  for  the  lumped-circuit  iterative 
models  given  in  (79)  for  a line  of  N sections: 


V 

= $ 
*kll 

^ fkl2  -i 

(¥^) 

(D-2a) 

_I  1 

= i 
^21 

* ikZZ  - 

. (D.2b) 

For  the  first  subsection  rewrite  (D-la)  as 

i«»  = Yfl  ^0  - 

where  Yg  = Zg”^  and  substitute  into  (D-2a)  with  k=l  to  yield 

V(X  /N)  = V(0»  + *^2  (Y„  - Yg  V(0»  (D-4) 

= «kU-iu2  2o>i<*»Mw2  2oIo 

which  is  the  first  equation  in  (80). 
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To  rierivn  th<!  last  equation  In  (80),  rewrite  (IX lb)  as 


U£)  = Yjj  V(£)  - ^ 


(D-5) 


where  Yj.=  Z£  , Substitute  this  into  (D-2b)  for  k = N to  obtain 

2x  :5Cl  = &21  y *)  + &22i  . (D-^l 


Substitute  (D-2a)  for  k = N into  (D-6)  to  obtain 


(D-7) 


Yj  V(X)  - V,  Vf  = v(V*)  + il.22{t;;}2  V(X)  . v(lBx)} 


which  can  be  rewritten  as 


{ fyzi  - ik22  !kU  } { !k22  SI2-  2x  }x(»  = 


However,  one  can  easily  prove  the  identity  in  (82) 


tk21  “ ik22  tkl2  ikll  ^ ~ ikl2 


- Yp  V,  (D-8) 

^ «U  '■  "'JL*  0, 


(D-9) 


associated  with  the  forms  of  the  lumped  iterative  matrix  chain  parameters 
in  (78).  Therefore,  substituting  (D-9)  into  (D-8)  and  multiplying  on  the  left 
by  -1.JJ2  equation  in  (80). 

The  derivation  of  the  intermediate  equations  proceeds  similarly. 
Substituting  (D-  2b)  for  k=m  into  (D-2a)  for  k=m+l  yields 

s)  = + ikl2  {iai  + *k22  i(-^  i)}  p-‘0> 


Writing  (D-2a)  for  k=m  as 


i = iw2  *)  - iila  ikii  1 ^ 


(D-11) 


and  substituting  into  (D-10)  yields 

{ ikl2  ^k21  " ^M2  hiZZ  ^klZ  ^kll  } - ("n"  ^ (D-12) 

~kll  ~kl2  ik22  ~kl2  } — W ^ 

- v(^~£)  =0 

— ' N ' ri-1 

Again  substituting  the  identity  in  (D-9)  into  (D-12)  we  obtain 

z(^^)-{»kU+«W2ik22iil2  <D-«) 

which  is  the  fornn  of  any  of  the  intermediate  equations  in  (80), 

In  addition,  it  can  easily  be  verified  from  the  forms  of  the  naatrix  chain 
parameters  for  the  lumped-circuit  iterative  models  given  in  (78)  that 


or 


£ A ^ $ $ 

~kl2  ~k22  ~k22  ~kl2 


A A & X ft  ^ 

-kl2  ~k22  *kl2  *k22 


(D-14a) 

(D-14b) 


Substituting  (D-14b)  into  (D-8)  (along  with  the  identity  in  (D-9))  and  substi- 
tuting (D-14b)  into  (D-13)  along  with  (D-4)  yield  the  final  nodal-admittance 
matrix  equations  in  (80): 


(ikiz  lo  - «kll)  = lkl2  Yo 


(D-IS) 
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- 'Iku  + iwz'  V i)  + V (i^  i)  = 0 

I **  * 

t 

^ <ikl2  - 4z> 

for  m=l, , N-1. 

The  forms  of  the  matrix  chain  parameters  for  the  lumped-circuit 
iterative  approximations  given  in  (78)  can  be  derived  in  the  following 
manner. 

For  the  lumped  ^ model,  the  terminal  equations  are  (see  Fig.  11(a)) 

iili  ■ i{”lT  ^)  ■ X (■§■)  ~ (w^)  • (D-16b) 

Equation  (D-I6a)  corresponds  to  the  first  matrix  chain  parameter  equation 
in  (78a).  Substituting  (D-l6a)  into  (D-I6b)  yields 

i (1^)  = -X  (i)X  ) 4 in  ^ 

which  is  the  second  matrix  chain  parameter  equation  in  (78a), 

For  the  lumped  F model,  the  terminal  equations  are  (see  Fig.  11(b)) 

-a(f)i(|[i)  (o->8=» 

i . (D-18W 

Equation  (D-18b)  is  the  second  matrix  chain  parameter  equation  in  (78b). 


(D-17) 


Substituting  (D-18b)  into  (D-18a)  yields 


}2(^)  -5('R)i(^") 


which  is  the  first  matrix  chain  parameter  equation  in  (78bj, 


(D- 19) 


For  the  lumped  Pi  model,  the  terminal  equations  are  (ser  Fig.  11(c)) 
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(D-20a) 


(I,")  = ■ <l)ii  - J Ki)  ^ (^*)} 


k-1 


1 „/£ 


^k-1 


= { J-n  H }v(^i)  -Z  (^i(^z) 


(D-20b) 


-ixy^'C— 0 • 

Equation  (D-20a)  is  the  first  matrix  chain  parameter  equation  in  (78c).  Sub- 
stituting (D-20a)  into  (D-20b)  yields 

which  is  the  second  matrix  chain  oarameter  equation  in  (78c), 

For  the  lumped  Tee  model,  the  terminal  equations  are  (see  Fig.  11(d)) 


■r  2(|)i(?,i) 


(D-22a) 


'-{I  z)  y(£)(v(>^x) 


'k-1 


k-1 


(D-22b) 
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Equation  (D-ZZh)  is  the  second  matrix  chain  parameter  equation  in  (78d). 
Substituting  (D-22b)  into  (D-22a)  yields 


{if] 


(D-23) 


which  is  the  first  matrix  chain  parameter  equation  in  (78d) 


APPENDIX  E 


The  purpose  of  this  appendix  is  to  provide  justification  i'or  omitting 

the  reference  potential  terms  in  the  potential  expressions  in  Chapter  V. 

Consider  Fig.  E-l(a)  in  which  infinitesimal  line  charges  lie  on  a cylindrical 

surface  of  radius  r . The  potential  (r  ,0  ) with  respect  to  the  potential 

w p p p 

reference  point  due  to  one  of  the  line  charges  is  {reference  [58],  pp.  91-92) 


Vi>-  8p'  = 


In 


(E-l) 


where  the  distances  from  the  line  charge  to  the  potential  and  reference 


points  are  given  by 

2^2 
d = r + r,„ 

p p w 


2r  r cos  {0-  0 ) 
p w ' p ' 


(E-2a) 


dy  = r^co3  (6-  0y)  . (E-2b) 

If  the  cylindrical  surface  supports  a per- unit-length  charge  distribution  of 
the  form 

A B 

o(6)  = + '^  a cos  m0  + S b sin  mft  . (E-3) 

w m=l  m m=l  m ' 

then  the  potential  0-(r  , 0 ) can  be  obtained  as  the  limiting  case  of  an 

P 

infinite  number  of  infinitesimal  line  charges  with  appropriate  weighting 
given  in  (E-3)  as  [56] 


^(rp  . ep)  = 0 in 


d0 


(E-4) 


Substituting  the  expressions  for  and  given  in  (£-2)  into  (E-4)  yield 

integrals  which  can  be  evaluated  in  closed  form.  The  result  is  [56] 

ip(rp  t Bp)  - &Q  Dq  (rp  , 0p,  r^  ) - aQ  DQ(rj, , , r^)  (E- 5) 

A 

■‘‘rfll  I « ®p»  ‘ ^nPiA^r  ' ®r*  J 

® f 1 

i ^ nr.  ®no(*’r  ^ ®r»  J 

where 


“o'-’p'  «p-  '•w' 


'w  >-1 


■'p>  “■» 


*■  r 
p « w 


(E-6a) 


r m+1 

\ cos  m0o  r > r 


^ 

m p p w 


2cm(rp)”^ 


p w 


2^-)m-l  cosmep  ^ 


(E-6b) 


D (r  , e , r ) 
m p p ’* 


2cm  (rp) 


m sin  me 


% >^w 


(E-6c) 


(r  jm 

— U. — , sin  m0^  r 4;  r 

2cm  P w 


The  third  argument  in  the  expressions  of  (E-6)  will  designate  the  radius 
of  the  boundary  svipportisig  the  charge  distribution.  The  terms  DQ(rj,,  i*^), 
0y,  r^)  and  r^)  in  (E-5)  are  the  reference  potential 

terms  which  were  omitted  from  the  potential  expressions  in  (96)  and  (97). 


Consider  a typical  system  of  (n-tl)  wires  shown  in  Fig.  E-l(b) 
bearing  per* unit- length  charge  distributions  of  the  form 

Ai  Bi 

Pi(fii)  = a.Q  a.j„  cos  m6i  + sin  me^  (E-7) 

for  i = 0,  I, , n,  A typical  expression  for  the  potential  at  a match 

point  on  the  U th  conductor  in  Fig.  E'>'i(b)  is 

*^i  “ ^00^0  ^ ’^pO'  ®p0'  *'w0^  “ ^OG  ^0  ' ®r0,  ^wO  ^ (E-8) 

^rEi  ^p0»^w0^“  ^Om  ^^rO' ®r0' ^vO^} 

Bq  si 

l^Om  ^m  '®p0'*’w0^  "^Om  ^m^’^rO’ ^rO' ^woU 


■*^*i0  ^0  ^**wi»  ®oi  * *wi  ^ “^0  ^0  ^**ri  ' ®ri ' ^wi^ 
r c 

^m=l  ^ m^^^wi  ’ %i  ’ *^wi^  "^im^ m ^’^ri ' ^ri  ’ Svi  ^ 1 

® i r ® 1 

+ S -lb.  D (r  .,0  .,r  .)-b.  D®(r.  ,0  .,r  .)  r 

m=l  ^ m ri  ri  wi  J 


' ^pj'  ^wj'  ”^j0  ^0  ^^rj  ' ®rj»  ^wj^ 

A. 

-f-  S fa.  ( r . , 9 . , r . ) - a.  ( r . , 0 . , r . ) 1 

m=l  ' rn'  oj  ’ ''pj  ' wj  ' jni  rn'  rj  ' '^rj'  wj' J 

B. 

’ ®oj  ' ^wj  ^ '^rj'^^wj^} 
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E 


■P|IPPH|PP|mpHn^^  imui  I j j 1 1 jiipiip. 


If  we  allow  the  reference  noint  for  the  oofcentials  to  move  to  infinity, 
the  reference  potential  terms  for  the  a^^cos  mf).  and  b.^  sin  mS.  terms, 

i.  e.  , (r  . , 0 . , r . ) and  D®  (r  . , 0 r .)  go  to  zero  for  i=0, 1, , n 

as  is  clear  from  (E-6b)  and  (E-6c).  In  addition,  if  the  total  ner-unit-length 
charge  on  the  system  of  (n+1)  conductors  is  zero,  then  the  reference  poten- 
tial terms  due  to  the  constant  expansion  terms  i.  e.  , Dq  (r 

ma  ,•  also  be  removed  for  all  i=0,l, , n.  This  can  be  shown  in  the  following 


manner.  The  total  per-unit-length  charge  on  the  i- th  conductor  is 

'2n 

^0 


\ f " Pi  («ii  >\.-i 

->n 


(E-9) 


= 2rr  r • a, . 
wi  lO 

and  the  portion  of  the  potential  expression  in  (E-8)  consisting  of  the  refer- 
ence potential  terms  due  to  the  constant  charge  expansion  terms  is 

---  - ^0 0 Dq  (rj.,0 , 0i.o»  

“ ■‘iO  ^0  ^^ri  ' ®ri  » ^wi^  “ 

---  -ajo  Eb  (rrj  , Gj-j  , r^j) 

n 

~ ^kO^O^^rk' ®rk' ^wk^ 

K=0 


(E-10) 


Utilizing  the  expression  for  of  the  form  given  in  (E-6a),  equation  (E-10) 


can  je  written  as 


n 

- r,  a 
k=0 


- r^i^  £n  (r^i^  ) 


) . 


(E-11) 


With  the  expression  for  the  total  per-unit-length  charge  on  the  i- th  conductor 
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given  in  (E-9),  equation  (E-11)  can  be  written  as 


n 

Y;  q 

k=0  2tt6 


Requiring  the  system  to  be  electrically  neutral,  i.e.  , 


n 


kit) 

equation  (E-12)  can  be  written  as 


in(r^Q) 

2ttc 


+ ?:  q^  <^rk) 

k=l  ^ 2tt  c 


(E-12) 


(E-13) 


(E-14) 


By  combining  associated  terms,  equation  (E-14)  can  be  written  as 


n 


^k 


Z/  —fL.  An 

1 2rr  e 


(E-15) 

As  the  reference  potential  point  moves  to  infinity,  the  distances  from  the 

centers  of  the  conductors  to  the  reference  point  become  equal,  i.e,, 

r^Q  = r^.^  = ----  = and  (E-15)  approaches  zero.  Therefore,  the 

reference  potential  terms  in  the  potential  expressions  may  be  omitted. 

Implicit  in  this  is  the  fact  that  the  potentials,  are  with  respect  to 

infinity.  This  is  permissible  as  was  shown  in  this  appendix  only  if  the 

net  per- unit-length  charge  on  the  system  is  zero,  i.e,, 

n 


k=0 


(E-16) 
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